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Thermodynamic  properties  have  been  determined  for  compound  semiconductors 
by  solid  state  electrochemical  technique  and  used  successfully  in  making  assessments  to 
the  available  data  on  these  systems.  The  measurement  technique  used  a solid  electrolyte 
and  two  custom  designed  unique  cell  configurations  to  maintain  stable  open  circuit 
electromotive  force  (emf)  values,  and  minimize  volatile  V-component  (-As,  -Sb)  loss  due 
to  vaporization  during  the  experimental  stage. 

The  experimental  data  on  Gibbs  energy  data  for  the  binary  systems  Ga-As,  In- As 
and  Sb-Te  are  the  first  of  their  kind  and  are  very  important  for  a good  understanding  of 
their  thermodynamic  properties  and  their  alloys.  In  addition  to  these  systems,  In-Sb  and 
Ga-Sb  systems  were  also  assessed.  The  assessments  of  the  III-V  binary  systems  were 
used  as  start  points  for  thermodynamic  predictions  into  the  two  ternary  systems,  namely 
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Ga-Sb-As  and  In-Sb-As.  Solid  state  emf  measurements  were  also  conducted  on  the 
ternary  Ga-In-As  to  understand  the  thermodynamic  characteristics  in  the  pseudobinary 
section  of  the  GaAs-InAs  phase  diagram.  The  experiments  confirmed  the  presence  of  a 
miscibility  gap  in  the  GaxIni.xAs  solid  phase. 
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CHAPTER  1 
INTRODUCTION 

Phase  Diagrams  and  Their  Importance  to  Compound  Semiconductors 

Phase  diagrams  are  graphical  representations  of  the  relations  between  state 
parameters  and  the  delineation  of  phase  fields  within  this  parameter  space.  Phase 
diagrams  are  useful  in  selecting  initial  compositions  and  subsequent  processing  conditions 
for  thermal  processes  including  the  separation  and  purification  of  materials  since  phase 
diagrams  are  a graphical  representation  of  the  conditions  of  phase  equilibria.  Extensive 
thermodynamic  information  can  be  extracted  from  them.  A thermodynamic  description  of 
a phase  is  related  to  the  interaction  of  the  species  and  the  structure  of  the  phase  at  the 
molecular  level.  Thus  phase  diagrams  can  be  calculated  from  known  thermodynamic 
phase  functions  and,  conversely,  thermodynamic  properties  and  phase  structure 
characteristics  can  be  obtained  from  phase  equilibria  data. 

Since  temperature,  pressure  and  concentration  are  typically  measurable  and 
controllable  state  parameters  in  many  processes;  the  Gibbs  energy  is  the  important 
thermodynamic  function,  the  value  of  this  function  being  a minimum  at  equilibrium  for 
constant  temperature,  pressure  and  overall  composition.  Graphically,  the  compositions  of 
two  phases  in  equilibrium  at  a certain  temperature  and  pressure  is  determined  by  the 
common  tangent  to  the  two  curves  representing  the  concentration-dependent  Gibbs 
energy  of  each  phase.  This  is  a graphical  interpretation  of  the  equilibrium  condition  that 
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the  chemical  potentials  of  the  species  in  the  coexisting  phases  are  equal.  It  is  noted  that 
the  pressure  dependence  in  solid-liquid  phase  equilibria  problems  is  often  assumed 
negligible  for  pressures  less  than  lO’  MPa. 

Of  interest  to  the  present  study  are  the  class  of  semiconductors  constituted  by  the 
elements  from  the  Groups  III  and  V columns  of  the  periodic  table,  and  termed  III-V 
compound  semiconductors.  These  compounds  exhibit  extended  miscible  substitutional 
solid  solutions  when  mixed  on  the  same  sublattice,  thereby  allowing  the  semiconductor 
device  designer  control  of  the  lattice  parameter  and  electrical  or  optical  characteristics  of 
the  material  to  conform  to  the  required  device  specification.  The  increased  number  of 
elements  in  III-V  compounds  as  opposed  to  the  elemental  nature  of  Si,  however,  make 
their  property  characterization  and  processing  technology  more  difficult. 

The  attractive  characteristics  of  compound  semiconductors  include  the  possibility 
of  direct  bandgap  energy  and  the  concomitant  high  electron  mobility.  As  a result,  these 
materials  are  important  in  high  speed  electronic  and  optoelectronic  device  applications. 
Furthermore,  these  materials  offer  a broad  selection  of  optical  and  electrical  properties 
that  are  often  superior  to  those  of  elemental  Si.  Group  III-V  compounds  have  been  used 
in  a wide  range  of  applications  as  the  host  semiconductor  material,  including  high-speed, 
optoelectronic  and  integrated  optical  devices  [Hei78,  Sew75  and  Yos75].  In  particular, 
(In,  Ga)  Sb  and  (In,  Ga)  As  form  important  constituents  in  the  solid  solutions  used  for 
long  wavelength  photodetectors,  light  emitting  diodes,  diode  lasers  and  heterogeneous 
bipolar  transistors.  Moreover,  the  semiconductor  InSb  has  the  lowest  bandgap  energy  (Eg 
= 0.18  eV)  and  highest  room  temperature  electron  mobility  (p.„  = 10^  cmW-s)  of  any  III- 
V compound  semiconductor  [Wel52,  Wel53  and  Bre54]. 
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The  processes  for  fabrication  of  compound  semiconductor  devices  often  require 
near  equilibrium  contact  between  solid  and  liquid  phases  (e.g.,  bulk  crystal  growth,  liquid 

phase  epitaxy)  or  between  solid  and  vapor  phases  (e.g.,  chemical  vapor  deposition,  rapid 

« 

thermal  annealing).  Phase  diagram  ai\d  thermodynamic  property  functions  are  important 
in  specifying  the  necessary  boundary  conditions  in  the  analysis  of  such  device  processing 
steps.  It  is  observed  that  the  optical  and  electrical  properties  of  compound 
semiconductors  are  sensitive  to  the  processing  conditions.  As  an  example,  variations  in 
the  melt  stoichiometry  during  bulk  crystal  growth  or  in  the  vapor  phase  IIW  input  molar 
ratio  during  CVD  significantly  changes  the  electrical  properties  of  the  grown  material. 
This  implies  that  extremely  accurate  thermodynamic  properties  and  phase  relations  are 
sometimes  necessary  for  their  meaningful  application  to  process  analysis. 

Typical  Phase  Diagrams  of  III- V Systems 

A typical  III- V binary  phase  diagram  is  shown  in  Figure  1 - 1 . It  is  noted  that  pure 
solid  elements  show  negligible  solubility  of  the  other  element  and  each  half  of  the  phase 
diagram  is  termed  simple  eutectic.  In  this  phase  diagram,  the  eutectic  towards  the  Group 
Ill-rich  region  of  the  phase  diagram  is  shown  as  being  nearly  degenerate,  which  means  that 
the  eutectic  composition  is  so  close  to  pure  Group  III  component  that  it  cannot  be 
depicted  in  the  figure.  The  eutectic  composition  is  determined  in  the  first  order  by  the 
differences  between  the  compound  and  pure  component  melting  temperatures,  and 
modified  by  the  liquid  solution  behavior  and  component  thermodynamic  properties.  For 
example,  the  relatively  close  melting  temperatures  of  GaSb  and  Sb  produces  a eutectic 
composition  that  is  relatively  rich  in  Ga.  The  eutectic  compositions  in  the  Group  V-rich 


4 


portion  of  the  III- As  or  III-P  phase  diagrams,  however,  can  be  experimentally  difficult  to 
determine  because  of  the  experimental  problems  encountered  by  high  vapor  pressures. 

The  upper  curve  in  Figure  1 - 1 stretching  across  the  limits  of  the  phase  diagram  is 
the  liquidus  curve  which  is  a boundary  of  the  two-phase  region  (liquid  + AB)  and 
separates  the  pure  liquid  phase  from  the  two-phase  region.  The  vertical  line  drawn  at  the 
composition  0.5  atom  fraction  of  the  phase  diagram  represents  the  position  of  the 
essentially  stoichiometric  compound  AB  (i.e.,  a line  compound).  In  actuality  the 
compound  AB  in  III-V  systems  is  a solid  solution  phase  with  negligible  extent  of 
solubility,  and,  when  comparing  this  boundary  to  other  phase  boundaries,  AB  is  treated  as 
a line  compound.  Although  the  solid  solution  extent  of  the  line  compound  has  negligible 
influence  on  the  other  boundary  positions,  the  extent  is  of  utmost  importance  to  the 
semiconducting  properties.  Any  deviation  in  the  ideal  stoichiometry  of  the  III-V 
compound  A0.5  B0.5  in  this  region  of  solid  solution  (typically  of  the  order  1 0'^  to  1 0'^  mole 
fraction)  is  the  result  of  native  point  defects.  Many  of  these  defects  are  ionic  and  can 
significantly  change  the  impurity  behavior  and  electrical  and  optical  properties  of  these 
materials.  Thus,  by  being  able  to  accurately  describe  the  thermodynamics  of  these  solid 
solutions  it  is  possible  to  predict  bulk  phase  equilibria  as  well  as  the  stability  and  extent  of 
solid  solution. 

Figure  1 -2(a)  shows  a typical  III-III  -V  ternary  phase  diagram,  where  A and  B are 
the  group  III  elements.  The  vertical  plane  through  the  elements  A and  C is  the  binary  III- 
V limit,  while  a similar  plane  through  the  points  B and  C is  the  binary  IIl'-V  limit.  Both 
the  binary  phase  diagrams  are  of  the  type  just  discussed  and  as  shown  in  Figure  1 - 1 . The 
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Figure  1 - 1 . Schematic  III-V  binary  phase  diagram. 
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AxBi_jjC  Pseudo— Binary  Liquidus 


L, 


ISOTHERMAL  LIQUIDUS  (c)  PSEUOOBINARY  CUT 
PROJECTIONS  ^ 

Figure  l-2(a).  Schematic  HI-III’-V  ternary  phase  diagram, 

(b) .  Isothermal  liquidus  projections  LI  and  L2, 

(c) .  Pseudobinary  AC-BC  phase  diagram. 
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vertical  plane  passing  through  the  points  A and  B represents  the  III-III  binary  systems 
which  exhibit  a range  of  phase  diagram  types.  The  surface  enveloping  the  entire  ternary 
diagram  is  the  liquidus  surface,  which  denotes  the  upper  temperature  limit  of  the  two 
phase  (liquid  +solid  phase/pure  element)  region.  The  eutectic  valley  located  in  the  Groups 

I 

III/III  -rich  region  of  the  phase  diagram  is  almost  degenerate,  while  the  region  in  the  phase 
diagram  towards  the  Group  V-rich  (i.e..  As,  Sb,  P)  is  not  well  characterized  owing  to 
problems  associated  with  the  high  vapor  pressures  of  Group  V elements  at  those 
conditions  of  temperature  and  composition. 

As  in  the  case  of  binary  phase  diagrams,  the  vertical  lines  at  the  Xv  = 0.5 
composition  in  the  binary  subsystems  A-C  and  B-C  are  the  locations  of  the  line  compound 
phases,  IIIV  and  IIIV  respectively.  The  vertical  plane  passing  through  the  compounds 
AC  and  BC  includes  a region  of  solid  solution  between  IIIV  and  III  V,  and  is  termed  the 
pseudobinary  cut.  When  the  ternary  phase  diagram  is  cut  by  a vertical  plane  through  AC 
and  BC  and  viewed  from  the  side,  the  typical  pseudobinary  phase  diagram  is  seen  as 
comprising  three  regions  sandwiching  a 'cigar'  shaped  two-phase  region  representing  solid- 
liquid  equilibrium  (isomorphous  phase  diagram).  The  phase  above  the  two-phase  region  is 
the  ternary  liquid  solution  of  fixed  Xv  = 0.5,  while  that  below  the  lower  curve  of  the  two 
phase  region  is  the  pseudobinary  solid  solution.  An  understanding  of  the  correlation 
between  structural  properties,  band  structures,  and  the  phase  diagram  will  aid  in 
evaluating  these  alloys  for  use  in  long  wavelength  IR  material  (LWIR),  lasers, 
photodetectors,  and  other  device  applications. 

The  systems  of  specific  interest  to  this  study  are  the  binary  Group  III  (Ga,  In) 
arsenides  and  antimonides,  and  the  associated  ternaries  (Ga-In-As,  Ga-Sb-As  and  In-Sb- 
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As).  The  binary  sytem  Sb-Te  is  also  addressed  in  this  study  because  of  its  application  as 
an  eraseable  optical  storage  and  recording  media.  The  motivation  of  this  research  is  to 
provide  accurate  thermodynamic  information  for  the  liquid  and  solid  phases  of  these 
technologically  important  semiconductor  systems.  This  information  will  be  helpful  to 
understand  the  interactions  in  the  alloys  at  the  molecular  level,  and  to  identify  the  design 
parameters  (e.g.,  temperature,  pressure,  compositions,  molar  flow  ratios,  growth 
conditions)  at  the  macroscopic  level  for  material  processing.  Accurate  thermodynamic 
properties  are  essential  in  selection  of  starting  materials  (source)  and  in  the  estimation  of 
extents  of  chemical  reactions.  Of  particular  importance  are  the  solid-gas  phase  reactions 
that  occur  during  chemical  vapor  deposition  (CVD)  and  their  manifestation  in  a CVD 
phase  diagram.  These  are  guiding  tools  for  establishing  the  boundaries  of  process 
operation  to  achieve  certain  film  properties. 

Thermodynamic  Information  From  Phase  Diagrams 

Phase  diagrams  can  be  thought  of  as  graphical  solutions  of  the  conditions  of  phase 
equilibria  and  metaphorically  as  storage  media  of  thermodynamic  properties.  It  is  often 
the  case  that  several  experimental  determinations  of  different  thermochemical  properties 
and  phase  diagrams  exist.  The  relationships  of  thermodynamics  provide  a means  to  check 
the  consistency  of  the  experimental  thermochemical  and  phase  diagram  data.  This 
methodology  of  performing  consistency  tests  by  the  process  of  comparing  predictions 
based  on  one  set  of  data  (e.g.,  phase  diagram)  to  experimental  data  based  on  other  sets 
(e.g.,  enthalpy  or  Gibbs  energy  data)  and  vice  versa,  has  been  illustrated  by  Aselage 


[Ase85]. 
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Compound  semiconductors  of  interest  in  this  study  pose  a situation  where  phase 
diagram  data  is  readily  avaliable  in  literature,  however  thermodynamic  properties  such  as 
vapor  pressure  and  enthalpy  data  is  not  so  accurately  available.  For  this  reason,  estimates 
have  been  made  by  numerous  authors  to  calculate  thermodynamic  properties  from  phase 
diagram  data.  Such  estimates  coupled  with  little  experimental  data  on  these  systems 
render  the  understanding  of  the  thermodynamic  properties  difficult.  The  large 
uncertainities  involved  in  such  estimates  make  it  challenging  to  reliably  perform 
experiments  for  measurements  of  thermodynamic  properties. 

The  data  by  Hoshino  et  al.  [Hos65]  on  In-Sb  system  indicate  that  the  In-Sb 
solution  is  not  subregular  as  earlier  assumed,  and  the  concentration  dependence  of  the 
exchange  energy  is  nonlinear.  Regular  solution  behavior  was  assumed  in  the  In-As 
solutions  near  the  melting  temperature  of  InAs  [Kar75].  In  the  Ga-Sb  system,  the 
experimental  values  of  excess  heat  and  entropy  of  solutions  have  not  been  well 
represented.  The  current  predictions  on  the  heat  of  mixing  differed  from  the  data  given  by 
Predel  and  Stein  [Pre71]  and  Gambino  and  Bros  [Gam75]  by  the  order  of  2 kJ/mol.  The 
excess  entropy  differed  from  the  data  of  Bergman  et  al.  [Ber74]  by  about  2.5  J/mol-K. 
The  values  of  heat  and  entropy  of  formation  of  solid  InAs,  GaAs  and  GaSb  have  not  been 
theoreticaUy  well  represented  [Glu72,  Luk77,  Pan72,  Gor75,  Yam88  and  Fol58].  In 
particular,  the  data  on  In-As  and  Ga-As  in  the  As-rich  segment  of  the  phase  diagram  is 
lacking,  and  the  values  that  are  available  are  subject  to  significant  dispersion  when 
compared  with  each  other.  Tmar  et  al.  [Tma84]  addressed  the  discrepancy  in  data  for  In- 
As and  Ga-As  systems  due  to  lack  of  information  on  gaseous  phase  composition  and  the 
species  contribution  to  the  total  pressure  in  the  gaseous  phase.  Both  As2  and  As4  are  the 


10 


main  constituents  of  the  vapor  phase,  depending  upon  the  location  in  the  phase  diagram. 
The  authors  calculated  partial  pressure  of  As2  and  AS4  based  upon  total  pressure 
measurements  of  Rau  [Rau75],  enthalpy  of  fusion  data  of  Kaufman  et  al.  [Kau81],  and  the 
data  of  Drowart  et  al.  [Dro78]  for  the  energy  of  dissociation  of  AS4  into  As2.  Such 
calculations  revealed  that  any  uncertainity  with  the  dissociation  enthalpy  has  negligible 
impact  on  the  phase  diagram  optimization,  while  has  a serious  impact  on  the  absolute 
partial  pressures  of  As2  and  AS4  species  calculated  above  the  In- As  and  Ga-As  systems.  In 
all  such  cases  direct  measurement  of  Gibbs  energy  data  is  desirable.  Such  is  possible  by 
emf  measurements  in  semiconductor  alloys  using  galvanic  cells,  the  cells  being  associated 
with  a functional  error  in  the  order  of  0. 1 to  0.3  %. 

Calculations  have  also  been  performed  for  activity  coefficients  of  Sb  and  Te  in  Sb- 
Te  alloys  [Gla69,  Gho89],  however  no  experimental  data  exist.  Furthermore,  the  Sb-Te 
system  is  associated  with  extremely  small  excess  entropy  of  mixing,  which  indicates  the 
predominance  of  entropy  increase  on  mixing  due  to  chemical  interactions  over  the 
decrease  due  to  the  resulting  ordering.  In  view  of  insufficient  thermochemical  data  on 
these  technologically  important  compound  semiconductor  systems,  an  attempt  is  made  in 
the  present  work  to  perform  experiments  and  assessments  to  obtain  the  most  consistent 
set  of  thermodynamic  data  in  order  to  characterize  these  sytems  fully. 

Solid  State  Activity  Measurements 

The  use  of  solid  electrolytes  in  galvanic  cells  has  been  a common  method  to 
measure  the  equilibrium  properties  of  electrode  materials.  The  critical  issues  in  solid  state 
ionics  are  material  stability,  transport  processes,  and  interfacial  phenomena.  The 
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application  of  solid  electrolytes  to  the  measure  Gibbs  energies  in  compound 
semiconductor  systems  is  discussed  next. 

Solid  Electrolytes 

The  solid  electrolyte  used  in  this  study  is  a solid  solution  between  Y2O3  and  Zr02. 
The  yttria  stabilizes  the  material  in  a fluorite  crystal  structure,  and  thus  is  termed  yttria 
stabilized  zirconia  (YSZ).  It  exhibits  abnormally  high  oxygen  ion  conductivity.  The 
fluorite  structure  is  a face-centered-cubic  arrangement  of  cations  with  the  oxygen  atoms 
occupying  the  tetrahedral  sites.  Such  a structure  has  a large  number  of  octahedral 
interstitial  voids.  In  YSZ,  the  size  of  the  O'^  ion  is  sufficiently  larger  than  the  metal  cation 
to  permit  easy  diffusion  of  O'^  ions  through  the  crystal  structure.  The  presence  of  8 to  10 
wt  % Y2O3  is  essential  to  stabilize  the  flourite  phase  in  zirconia  for  use  at  temperatures 
below  1500°C. 

The  conductivity  mechanism  of  YSZ  is  based  upon  its  defect  structure.  The  most 
common  defect  types  are  the  Schottky  defects  (vacancy  pairs)  and  Frenkel  defects 
(vacancy-interstitial  pairs).  The  total  electrical  conductivity,  a,  through  the  YSZ  is  the 
sum  of  the  effects  due  to  the  concentrations  of  ionic  and  electronic  (electrons  and  holes) 
defects  in  the  solid,  and  is  given  by 

ct  = Iciqj/Xi +nqeMe  + pqhMh  (1-1) 

i 

where  c is  the  ionic  defect  concentration,  q is  the  charge,  fj.  is  the  mobility,  and  n and  p 

are  the  electron  and  hole  concentrations,  respectively.  The  subscripts  i,  e and  h represent 
ions,  electrons  and  holes,  respectively.  Useful  application  of  solid  electrolytes  requires  the 
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electron  and  hole  conductivity  to  be  considerably  less  than  the  total  conductivity.  That  is 
the  case  in  YSZ  that  oxygen  ion  conductivity  is  essentially  the  only  mechanism. 

Principle  of  Electrochemical  Cells 

The  activity  of  component  A in  a binary  A-B  solution  was  determined  in  an 
electrochemical  cell  of  the  type: 

A,  AOy/A"SZ//AyBl.y,  AOy  (I) 

In  the  experimental  measurements  presented  here,  two  types  of  measurements 
were  made.  In  one  configuration  component  activities  in  solution  were  determined  and  in 
the  other  the  Gibbs  energy  of  compound  formation  was  measured. 

The  reference  electrode  consists  of  a two  phase  mixture  of  A and  its  stable  oxide 
AOx,  while  the  working  electrode  contains  the  solution  AyBi.y  and  the  most  stable  oxide 
AOx.  The  half  cell  reactions  and  the  overall  rections  are: 

A (pure)  + X O'^  = AOx  + 2x  e'  (la) 

AOx  + 2x  e'  = A (in  solution  AyBi.y)  + x O'^  (Ib) 

A (pure)  = A (in  solution  AyBi.y)  (Ic) 

The  overall  cell  reaction  is  simply  the  transfer  of  A in  the  pure  (reference)  state  to 
the  solution  phase.  The  Gibbs  energy  change  of  the  reaction  is: 

AGic  = AG°  + RTlnaA  (1.2) 

when  specifying  the  phase  of  the  reference  state  the  same  as  that  of  pure  A at  the 


measurement  temperature. 
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The  value  of  AG  is  obtained  from  the  open  circuit  potential  according  to  the 
following  equation: 

AG  = -nFE  = -2x  FE  (1.3) 

Therefore  the  activity  of  A in  solution  AyBi.y  is  directly  calculated  from  the 
measured  open  circuit  potential  by 

aA  = exp  (-2x  FE/RT)  ( 1 .4) 

Another  view  of  this  process  is  to  consider  the  associated  reaction  occuring  at  both 
electrodes 

X O-'  = Vi  O2  (g)  (Id) 

Since  the  aA  in  the  reference  electrode  (pure  A)  is  higher  than  that  in  the  working 
electrode,  the  partial  pressure  of  O2  that  would  be  in  equilibrium  with  each  electrode  is 
lower  in  the  reference  electrode.  This  chemical  potential  difference  causes  oxygen  ions  to 
transfer  from  the  working  electrode  to  the  reference  one. 

The  transport  process,  however,  results  in  an  electron  accumulation  at  the  low 
oxygen  chemical  potential  electrode  (left  side),  and  electron  depletion  is  produced  at  the 
electrolyte  surface  at  the  high  oxygen  chemical  potential  electrode  (right  side).  This 
produces  an  electric  potential  gradient  across  the  electrolyte.  At  the  point  of  equilibrium, 
the  chemical  potential  gradient  of  the  oxygen  ions  is  balanced  by  the  electric  potential 
gradient  across  the  electrolyte  so  that  O'^  electrochemical  potential  gradient  is  zero.  This 
balance  at  equilibrium  is  mathematically  represented  by  the  following  equation: 

H 

-nFE=  Jti„„d/t 

# ^ 

a 

02 


(1.5) 
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where  tion  is  the  ionic  transference  number  defined  as  the  fraction  of  ionic  conductivity  to 
the  total  electric  conductivity.  When  tion~  i the  above  relation  is  reduced  to 

-nFE  = /io,  - Mo,  ( 1 .6) 


where  E is  the  measured  open  circuit  emf,  F is  the  Faraday  constant,  n is  the  number  of 
electrons  participating  in  the  half  cell  reaction.  Thus,  by  careful  design  of  such 
electrochemical  cells  the  measurement  of  the  open  circuit  emf  is  related  to  the  chemical 
potential  difference  for  the  transported  species  (O'^).  If  the  chemical  potential  is  well 
known  at  one  electrode  (termed  reference  electrode),  then  the  only  unknown  is  the 
chemical  potential  of  the  transported  species  at  the  other  electrode  (termed  working 
electrode).  The  difference  between  reactions  la  and  Id  is 

A + x/2  O2  = AOx 

and  the  sum  of  Ib  and  Id  is 

AOx  = A (in  solution  AyBi.y)  + x/2  O2 

The  overall  cell  reactions  are  identical  and  thus  equation  1 .4  is  the  same  as  equation  1 .6. 

The  Gibbs  energy  of  formation  of  compound  AB  was  measured  in  a cell  of  the 


type: 


A,  AOx//YSZ//AB,  B,  AOx 

(II) 

A + x/2  O'^  = AOx  + 2e' 

(Ila) 

B + AOx  + 2e  = AB  + x/2  0'^ 

(Ilb) 

A(aA)  + B(aB)  = AB(aAB)  (lie) 

Expressing  the  reactions  in  which  all  components  in  their  standard  states  (aj  = 1) 
are  reversibly  changed  from  the  reference  state  (aj  = 1)  to  the  state  at  which  the 
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measurement  is  made(ai),  the  Gibbs  energy  associated  with  each  change  is  expressed  as 
follows: 


A(aA)  = A(aA  = 1) 


AG,.7  = (G°A-GJ  = RTln 


V aA  / 


(1.7) 


B(aB)  = B(aB=  1) 


AG,.8  = (G”-GB)  = RTln 


V aB  / 

AB(aAB  = 1 ) = AB(aAB)  AG, , = (G^^  " G“ab)  = RT  In  a^^ 

Along  with  the  standard  Gibbs  energy  of  formation  of  compound  AB : 

A(aA  = 1)  + B(aB  = 1)  = AB(aAB  = 1)  AG,.,o=  AG”,^ 


(1.8) 

(1.9) 


(1.10) 


Upon  summation  of  reactions  (1.7)  through  (1.10),  the  overall  cell  reaction  lie  is  obtained 
and  the  expression  for  the  standard  Gibbs  energy  of  the  reaction  (1.10)  is  established  as 
follows: 


AG,6  = 0 = RTln 


V J 


+ ^G°,o 


(1.11) 


The  value  of  AG°,q  is  determined  experimentally  by  measurement  of  the  open  circuit 


potential  and  is  related  to  the  activity  values  through  the  equation  (1.11). 


CHAPTER  2 

EXPERIMENTAL  STUDIES 
Description  of  Apparatus  and  Measurement  Techniques 

The  solid-state  electrochemical  cell  enclosure  assembly  used  in  this  study  is 
depicted  in  Figure  2-1.  The  cell  enclosure  was  a closed-one-end  alumina  tube  45.7cm  in 
length  and  5.08  cm  outside  diameter.  It  was  secured  onto  a water-cooled  brass  cell  head 
by  means  of  a Viton  0-ring.  Four  0.318cm  cajon  through-bore  fittings  were  located  on 
the  cell  head  while  another  0.318cm  through-bore  fitting  was  located  at  the  center  of  the 
cell  head.  Alumina  sheaths  were  passed  through  these  fittings.  The  central  fitting  passed 
an  alumina  sheath  (single  bore,  0.318  cm  O.D.)  through  which  an  inconel  wire  was 
inserted  to  suspend  the  cell  assembly  within  the  enclosure.  Inconel  lead  wires  were  passed 
through  a double  bore  alumina  sheath  through  one  of  the  outer  four  fittings  on  the  cell 
head  to  make  connections  to  the  reference  and  working  electrodes  of  the  cell.  A 
thermocouple  (Type  S:  Pt-10%  Rh,  Pt)  was  passed  through  another  outer  fitting.  The 
remaining  two  bores  on  the  brass  cell  head  were  provided  for  the  argon  purge  gas  inlet 
and  outlet.  The  exposed  ends  of  the  alumina  sheaths  were  sealed  by  means  of  epoxy  resin. 
The  thermocouple  tip  was  placed  in  an  alumina  casing  and  sealed  by  means  of  high 
temperature  alumina  cement. 

The  temperature  of  measurement  was  typically  in  the  range  773  to  1473  K.  At 
high  temperature,  there  was  concern  that  volatilization  losses  would  affect  the  working 
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electrode  composition.  For  this  reason  two  different  cell  constructions,  designed  to 
minimize  electrode  material  evaporation,  were  used. 

The  cell  design  designated  Type  I and  termed  'compression  design'  is  a 
modification  of  the  original  design  developed  by  Egan  and  coworkers  [Del79,  Alq83, 
Ega85],  and  is  shown  schematically  in  Figure  2-2.  Alumina  disks,  2.0  cm  diameter  and  0.6 
cm  thick  with  a countersunk  hole  (1.0  cm  diameter  and  0.3  cm  in  depth),  were  used  as 
crucibles  for  the  electrode  material  and  were  custom  fabricated  by  Bolt  Technical 
Ceramics  (Conroe,  Texas  USA).  The  countersunk  holes  of  the  necessary  diameter  were 
created  by  ordering  a special  tool  incorporating  a diamond  cutting  edge  on  the  outer 
surface  of  the  cylindrical  drill.  The  working  electrode  was  comprised  of  a mixture  of  the 
metal  alloy  and  the  stable  oxide,  while  the  reference  electrode  contained  a mixture  of  the 
pure  metal  and  metal  oxide  powder.  A yttria  stabilized  zirconia  (YSZ,  8 to  10  % yttria  in 
zirconia)  disk,  and  2.54  cm  diameter  and  0.16  cm  thick,  was  used  as  the  oxide  solid 
electrolyte,  supplied  by  Zircoa,  Inc.  (Solon,  Ohio).  The  faces  of  the  alumina  crucibles  and 
YSZ  electrode  were  carefully  polished  using  a carbide  coated  grinding  conveyer  belt  to 
roughen  the  contact  surfaces  to  provide  good  seal  between  the  electrolyte  and  the 
crucibles.  The  alumina  crucibles  containing  the  electrode  materials  were  placed  on  either 
side  of  the  YSZ  electrolyte.  The  cell  was,  in  turn,  pressed  between  stainless  steel  plates 
and  the  entire  assembly  was  held  together  by  four  tantalum  rods. 

A tungsten  foil  0.005  cm  thick  with  a hole  ~ 0.6  to  0.8  cm  diameter  was  placed 
between  the  electrolyte  and  each  of  the  crucibles.  The  diameter  of  the  hole  in  the  tungsten 
foil  was  slightly  smaller  than  the  1.0  cm  countersunk  hole  in  the  crucible.  This  was 
necessary  because  after  filling  up  the  crucibles  with  electrode  material,  an  overlap  of 
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inconel  Suspension  Wire 


Figure  2-1.  Schematic  of  the  solid  state  electrochemical  cell  assembly. 
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Alumina  Crucible 
(Working  Electrode) 


Tungsten  Foil 


YS2  Electrolyte 


Tungsten  Foil 


Alumina  Crucible 
(Reference  Electrode) 


Stainless  Steel  Plate 
Molybdenum  Cap 


Stainless  Steel  Spacer 


Stainless  Steel  Ball 
Molybdenum  Plate 
Stainless  Steel  Plate 


Tantalum  Rod 
and  Nut 


Figure  2-2.  Schematic  of  the  compression  design  for  the  electrochemical  cell  (type  I). 
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TOP  VIEW  SjDE  VIEW 


EXPLODED  SIDE  VIEW 


Figure  2-3.  Schematic  of  the  encapsulated  design  for  the  electrochemical  cell  (type  II). 
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about  0.2  to  0.4  cm  of  tungsten  metal  was  maintained  for  ensuring  electrical  contact  with 
the  electrodes.  The  rough  surfaces  of  the  crucibles  and  the  electrolyte  gave  a good  seal 
when  the  tungsten  foil  was  placed  in  between  the  electrode  and  the  electrolyte,  and  the 
tantalum  rods  were  screwed  on  tight.  The  choice  of  tungsten  foil  over  commonly  used 
platinum  was  based  upon  the  fact  that  platinum  when  used  in  such  experiments  was  readily 
involved  in  the  reaction  with  the  vapor  phase  of  the  electrodes  causing  severe  damage  to 
the  point  of  breakage.  Furthermore,  platinum  being  significantly  more  expensive  than 
tungsten  was  not  a good  choice.  A test  cell  was  run  with  to  verily  the  acceptable  use  of 
tungsten  foils.  It  was  observed  that  even  though  the  tungsten  was  readily  oxidized,  the 
measured  voltage  difference  across  the  electrodes  changed  by  less  than  0.5%.  The  reason 
being  that  the  oxidation  reactions  in  the  two  compartments  had  a natural  tendency  to 
cancel  each  other’s  effect.  Tungsten  foils  of  thickness  of  the  order  of  0.005  cm  were 
optimum  in  their  performance  as  electrical  contacts,  especially  because  of  their  greater 
capacity  to  deform  compared  to  thicker  foils  which  being  brittle  are  subject  to  cracks 
under  the  effect  of  mechanical  stress.  Care  was  taken  to  ensure  that  the  centers  of  the 
holes  m the  alumina  crucibles  and  tungsten  foils  were  aligned  with  each  other  while 
assembling  the  working  and  reference  electrodes. 

This  cell  design  relied  on  the  difference  in  thermal  expansion  between  the  Ta  rods 
and  stainless  steel  plugs  to  form  a compression  seal  between  the  electrode  crucible  rim  and 
the  electrolyte  disk.  The  celt  was  held  aligned  with  the  vertical  axis  of  suspension  and  was 
imparted  some  degree  of  rotational  freedom  by  means  of  a 0. 1 cm  diameter  stainless  steel 
ball  seated  at  the  base  of  the  molybdenum  plate.  This  was  necessary  because  of  the  brittle 
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nature  of  the  YSZ  electrolyte.  Any  additional  stress  due  to  uneven  mechanical  tightening 
of  the  tantalum  nuts  could  cause  the  electrolyte  to  break. 

Inconel  lead  wires  were  well  suited  for  spot  weld  contact  with  the  tungsten  foils 
and  was  used  instead  of  tungsten.  Tungsten  wire  posed  a problem  during  spot  welding, 
because  of  its  high  melting  temperature  (~  3683  K)  and  was,  thus,  replaced  by  inconel.  It 
was  determined  by  electrical  measurements  of  cell  resistance  that  the  spot  weld  contact 
and  the  lead  wire  IR  drop  did  not  seriously  affect  the  voltage  readings  and  that  the 
variations  were  less  than  1%  when  compared  to  the  observed  voltages  during  the 
experiment. 

Cell  Type  II,  also  termed  'encapsulated  design',  was  a cell  arrangement  in  which  a 
maximum  of  five  alloys  compositions  could  be  studied  simultaneously.  As  shown  in 
Figure  2-3(a)  and  (b),  the  cell  holder  was  designed  with  five  cylindrical  holes  (15  mm 
diameter  and  20  mm  depth)  along  the  periphery.  These  holes  housed  the  YSZ  electrolyte 
crucibles  that  were  cylindrical,  flat-one-end  closed  tubes  (13  mm  diameter  and  50  mm 
length).  The  annular  space  between  the  YSZ  crucible  and  alumina  cell  holder  contained 
the  reference  electrode  material  comprised  of  the  metal  and  the  metal  oxide.  The  working 
electrode,  which  consisted  of  the  alloy  of  the  desired  composition  and  the  most  stable 
oxide,  was  placed  inside  the  YSZ  crucible.  Tungsten  wires  were  used  as  the  leads  for 
both  electrode  compartments.  The  open  ends  of  the  YSZ  tubes  were  sealed  by  means  of 
tin  caps  pasted  onto  the  tubes  using  alumina  cement. 

Owing  to  the  open  design  of  such  a cell  arrangement  the  problem  arose  in  not 
being  able  to  contain  the  vapors  of  the  more  volatile  component  of  the  alloy  under 
investigation,  especially  at  high  temperature.  This  was  avoided  by  the  incorporation  of 
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B2O3  powder  (typically  0.5  g)  in  both  the  electrode  compartments.  B2O3  melts  at  about 
450°C  and  is  a viscous  liquid  immiscible  with  the  liquid  metal  electrode  materials.  It  is 
believed  that  the  liquid  B2O3  fully  encapsulates  the  electrode  material  and  provides  an 
effective  mass  transfer  barrier  to  As  or  P escape  from  the  melt.  This  approach  to  reducing 
volatilization  losses  has  long  been  used  in  the  bulk  crystal  growth  of  group  III  arsenides 
and  phosphides.  Several  test  cells  were  assembled  to  quantify  the  influence  of  B2O3 
encapsulant  on  the  experimental  results.  The  following  systematic  cell  was  used: 

A,  hOJIYSZIIk,  AOx,  B2O3  (III) 

The  cells  were  run  for  both  indium  and  gallium  (A=In,  Ga)  at  the  typical 
temperature  range  773  to  1 123  K and  open  circuit  emfs  were  recorded.  The  values  were 
of  the  order  of  1 to  4 mV  (equivalent  to  a maximum  value  of  1.16  kJ/mol)  amounting  to 
less  than  5%  of  the  experimental  measurements  in  such  cells.  With  B2O3  in  both  the 
compartments,  the  measured  voltage  difference  was  less  than  1 mV  which  corresponded 
to  an  error  of  less  than  0.5%.  The  effect  of  B2O3  in  both  the  compartments  was 
determined  to  well  offset  each  other  so  that  the  impact  of  B2O3  was  presumed  negligible. 

The  cell  assembly  was  heated  in  a resistance  furnace  that  consisted  of  two  heating 
zones,  the  temperature  of  each  zone  being  controlled  by  a microprocessor-based  digital 
process  controller  (Micristar  Model  828D,  Control  Systems,  Research  Inc.).  The  furnace 
had  the  capacity  to  reach  a temperature  of  1400°C  and  maintain  temperature  uniformity  to 
within  ± 0.5  K.  The  entire  cell  assembly  was  blanketed  in  a flowing  gas  stream  of  Ar. 
This  purge  stream  was  purified  by  passing  the  gas  sequentially  over  hot  Ti  sponge  (1123 
K),  a molecular  sieve  trap,  and  a liquid  nitrogen  cold  trap. 
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The  cell  emf  was  measured  with  a Keithley  Model  1 96  digital  multimeter  of  high 
input  impedance.  A pressure/flow  controller-Type  250B  (MKS  Instruments,  Inc.)  was 
used  to  monitor  the  gas  pressure  and  purge  gas  flowrate.  The  pressure  in  the  cell 
assembly  could,  however,  be  manually  adjusted  with  the  inlet  and  exit  valves  to  the 
system.  Any  pressure  in  excess  of  the  preset  value  was  released  through  the  automatic 
on/off  pneumatic  valves  triggered  by  the  Micristar  process  controller. 

The  experimental  apparatus,  as  outlined  in  the  previous  section,  is  schematically 
shown  in  Figures  2-1  to  2-3.  Electrodes  of  the  cells  were  prepared  from  components  in 
pure  states.  They  were  mixed  together  in  the  electrode  crucibles  and  prepared  to  give  a 
well-packed  mixture  in  the  electrode  compartments.  The  electrode  crucibles  were  then 
physically  contacted  with  the  YSZ  electrolyte  to  achieve  the  desired  cell.  Prior  to  the  cell 
operation  it  was  necessary  to  leak  check  the  assembly.  This  was  achieved  by  a series  of 
successive  evacuation  and  purging  steps  and  monitoring  the  pressure.  Evacuation  was 
performed  by  a high  vacuum  system  to  pressure  levels  in  the  range  lO*^  to  lO'^  Torr.  An 
acceptable  leak  proof  cell  condition  was  determined  when  the  cell  assembly  pressure 
stayed  within  the  vacuum  range  and  exhibited  less  than  0.0001  Torr  increase  in  pressure 
over  a duration  of  120  s.  High  purity  argon  was  used  as  the  inert  purge  gas  with  the 
purge  flow  rate  adjusted  between  1 to  2 standard  cm^/s.  The  cell  temperature  was 
programmed  to  rise  at  the  rate  of  0.02  Ks‘l  to  the  measurement  temperature. 

Measurements  were  initiated  typically  after  2 to  3 days  when  equilibrium  was 
achieved  at  the  first  temperature  setting.  The  cell  was  determined  to  be  at  equilibrium 
when  the  open  circuit  emf  values  were  stable  (within  0. 1 mV  over  a period  of  two  hours). 
Cell  reversibility  was  verified  by  applying  a small  current  to  the  cell  in  both  directions  and 
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verifying  that  the  cell  voltage  returned  to  its  original  value.  The  cell  was  also  observed  to 
be  reversible  for  small  temperature  increments  and  decrements  (order  - 5 K).  At  the 
conclusion  of  the  experiment,  the  cell  compartments  were  again  checked  for  leaks  by 
successive  evacuation  and  purge  steps  and  pressure  monitoring.  The  electrodes  were 
examined  for  evidence  of  oxidation  and  were  analyzed  by  X-ray  diffraction  to  confirm  the 
absence  of  side-reactions. 


Experimental  Results  in  Binary  Systems 

As  discussed  in  Chapter  1,  the  indium  and  gallium  arsenides  and  antimonides  are 
the  systems  of  interest  in  this  study.  These  materials  are  the  building  blocks  of  the 
semiconductor  alloys  used  in  the  optoelectronics  industry.  The  electrical  and  optical 
properties  of  these  alloys  depend  on  the  stoichiometry  of  the  alloys  and  since  they  are 
often  grown  with  near  equilibrium  boundary  conditions  the  alloy  compositions  can  be 
estimated  by  equilibrium  considerations.  Activity  measurements  in  both  solid  and  liquid 
phases  of  these  alloys  have  been  performed  in  this  study  using  solid  state  galvanic  cells. 

Ga-As  System 

Phase  diagram  data.  The  phase  diagram  for  the  Ga-As  system  has  been  fairly  well 
studied  by  numerous  authors  [Sol74,  Kos55,  Hsi74,  Hal63,  Rub66,  Rak70,  Vig73,  PerSO, 
DeW86,  Dut83,  Ceo89,  Yam88,  Ans84],  and  has  been  shown  to  be  characterized  by  the 
existence  of  the  congruently  melting  compound  GaAs.  Because  of  high  pressures,  the 
phase  diagram  region  close  to  pure  As  composition  is  not  yet  been  accurately  determined. 
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Most  of  the  data  are  in  fairly  good  agreement  except  for  the  results  of  Koster  and  Thoma 
[Kos55].  A summary  of  literature  data  is  presented  in  the  next  chapter. 

Calorimetric  and  enthalpy  data.  The  enthalpies  of  formation  and  fusion  as  well  as 
low  temperature  Cp  value  for  GaAs  have  been  reviewed  by  Sirota  [Sir68],  Martosudirjo 
and  Pratt  [Mar74],  Piesbergen  [Pie63],  Holste  [Hol72],  Cox  and  Pool  [Cox67],  Dash  et 
al.  [Das74],  Lichter  and  Sommelet  [Lic69]  and  Richman  and  Hockings  [Ric65]  by 
calorimetric  methods  such  as  calorimetric  precipitation,  drop  calorimetry,  differential 
thermal  analysis,  thermometry  and  oxygen  calorimetric  bomb.  The  mixing  enthalpy  for  the 
liquid  has  been  determined  by  Yamaguchi  et  al.  [Yam88]  by  drop  calorimetry  over  a 
temperature  range  of  800  to  1 550  K.  The  heat  capacity  data  for  GaAs  compound  [Cox67, 
Das74  and  Lic69]  show  large  discrepancies.  Due  to  this  reason,  Chatillon  et  al.  [Cha90] 
performed  a least  square  fit  to  the  heat  capacity  data  of  these  authors  to  obtain  an 
expression  for  optimized  heat  capacity  data  in  the  temperature  range  1 50  to  1511  K.  This 
expression  was  used  in  performing  third  law  analysis  of  the  measured  Gibbs  energy  data 
for  Ga-As  system  in  this  study. 

Vapor  pressure  and  Gibbs  energy  measurements.  The  equilibrium  vapor  pressure 
above  solid  GaAs  and  liquid  solutions  has  been  measured  by  several  researchers.  Above 
1510  K,  GaAs  vaporizes  incongruently  with  an  excess  of  As.  Although  several  polymers 
of  As  exist,  the  predominant  species  below  the  melting  temperature  are  As2  and  As4. 
Unless  care  is  taken  to  maintain  a particular  stoichiometry,  when  the  temperature  is  above 
the  maximum  congruent  vaporization  temperature,  the  preferential  loss  of  As  will  lead  to  a 
two-phase  Ga-rich  melt-GaAs(s)  equilibrium.  Both  partial  pressure  and  total  pressure 
measurements  have  been  reported  in  the  literature.  The  partial  pressures  of  As2  and  As4 
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have  been  measured  [Pup74,  Dro58,  Gut61,  Art67,  DeM70,  Fox74  and  Hir85]  using 
Knudsen  cell  mass  spectrometry.  Total  arsenic  pressure  measurements  have  been 
performed  by  van  den  Boomgaard  and  Schol  [B006I],  Folberth  [Fol58],  Lyons  and 
Silvestri  [Lyo61]  using  the  dew  point  technique,  by  Rakov  et  al.  [Rak70]  using  continuous 
weighing  technique,  by  Vigdorovich  et  al.  [Vig73],  Richman  [Ric63],  Hegewals  et  al. 
[Heg88]  using  Bourdon  gauge  and  by  Khukhryanskii  et  al.  [Khu74],  and  Panish  [Pan74] 
using  dynamic  methods  incorporating  transport  by  H2  gas  flow  or  equilibrium  with  H2  and 
AsHa  gas  flows.  No  experimental  information  is  available  on  the  standard  Gibbs  energy  of 
formation  of  GaAs  (s)  at  a temperature  above  743  K.  The  Gibbs  energy  of  formation  of 
GaAs(s)  was  measured  by  Abbasov  et  al.  [Abb66],  Krestovnikov  et  al.  [Kre71],  and  Sirota 
[Sir68]  by  electrochemical  method  in  the  temperature  range  637  to  743  K.  These  emf 
measurements  were  performed  at  low  temperature  and  were  found  inconsistent  with  the 
phase  diagram  data  [Cha90].  Recently,  Lee  [Lee92]  conducted  emf  measurements  for  Ga 
activity  along  the  Ga-rich  liquidus  curve  (two-phase  equilibrium  GaAs-liquid)  in  the 
temperature  range  923  to  1373  K. 

There  exists  significant  discrepancy  between  the  vapor  pressure  data  of  the 
numerous  authors  [Dro58,  Gut61,  Art67,  DeM70,  Fox74,  Fox73,  Rub66,  Boo57,  Fol58, 
Lyo61,  Rak70,  Ric63,  Hsi74,  Khu74a,  Khu74b  and  Pan74].  A major  problem  has  been 
the  complications  associated  with  the  As  vaporization  and  the  subsequent  formation  of 
As2  and  As4  species.  The  isolation  of  the  gaseous  species  and  their  individual  contribution 
to  the  total  measurement  in  the  mass  spectrometric  measurements  and  vapor  pressure 
studies  is  not  easily  quantified.  Depending  upon  the  temperature,  either  As2  or  AS4  is  the 
main  constituent  of  the  vapor  phase  [Art67,  Ric63,  Pup74,  Dro58,  DeM70,  Pan74,  Lyo61 
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and  Rak70].  For  the  total  pressure  measurements  or  when  only  one  species  partial 
pressure  is  measured,  the  remaining  species  partial  pressures  are  calculated  from  the 
equilibrium  constant  Kp  for  the  dissociation  reaction: 


One  of  the  difficulties  in  comparing  the  data  is  the  large  uncertainity  in  the  value  of 
Kp.  Chatillon  et  al.  [Cha90]  have  recalculated  the  literature  As  vapor  pressure  data 
referred  to  As(s)  and  have  presented  an  equation  for  calculations  based  upon 

transformation  equation  involving  Kp.  Even  after  such  corrections  were  applied  to  the 
vapor  pressure  data,  significant  disagreement  still  existed.  This  being  the  case,  the 
motivation  behind  this  study  was  to  measure  the  standard  Gibbs  energy  of  formation  of 
GaAs  (s)  using  a solid  state  electrochemical  cell.  Cell  type  II  was  used  to  extend  the 
temperature  limit  of  measurement  beyond  that  of  the  previous  molten  salt  studies. 

Standard  Gibbs  Energy  of  Formation  of  GaAs  (s') 

The  cell  that  was  used  in  the  experiment  is  shown  below  along  with  the  half  cell 
reactions: 


As2(g)  = |As^(g) 


(2.1) 


B203(1),  GazOaCs),  Ga  (1,  pure)/AfSZ//GaAs  (s),  As(s),  GaiOjCs),  BjOjCI)  (IV) 


Ga  (1,  pure)  + 3/2  O'^  = 1/2  Ga203  + 3e- 


(2.2) 


1/2  Ga203  + As  (s,  pure)  + 3e-  = GaAs  (s)  + 3/2  O'^  (2.3) 


Ga  (1,  pure)  + As  (s,  pure)  = GaAs  (s) 


(2.4) 
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Applying  Wagner's  approximation  [Wag57]  to  the  chemical  potential  of  the 
essentially  stoichiometric  compound  GaAs  the  overall  cell  reaction  is  transformed  to  the 
following: 

Ga  (1,  pure)  = Ga  (ss)  (2.5) 

where  ss  denotes  the  solid  solution  Ga-As  with  composition  fixed  to  lie  along  the  solidus 
line  towards  the  GaAs(s)  + As(s)  region.  This  is  equivalent  to  measuring  the  chemical 
potential  of  Ga  along  the  As-rich  solidus  of  the  compound  GaAs.  The  standard  Gibbs 
energy  of  formation  of  GaAs  referred  to  Ga  (1)  and  As  (s)  states  is  given  by  the  following 
relation: 

A“Gf(GaAs)=‘>Go,^,3-“G«3, = -3FE  (2.6) 

where  F is  the  Faraday  constant  (96,485  C/eq-mV)  and  E is  the  open  circuit  emf  (mV). 

The  measured  values  of  the  open  circuit  potential  at  the  various  measurement 

temperatures  are  listed  in  Table  2-1.  The  emf  values  and  the  calculated  A“Gj(T)  values 

have  a near  linear  temperature  dependence  as  shown  in  Figure  2-4.  The  A°Gf(T)  values 

are  significantly  lower  (more  negative)  than  the  earlier  reported  emf  measurements  when 
extrapolated  to  this  temperature  range.  One  possible  reason  is  the  uncertainity  in  the 
oxidation  state  of  Ga  in  the  molten  halide  cells  of  the  earlier  work.  Gallium  has  both  a +1 
and  +3  oxidation  state  and  the  Gibbs  energy  difference  between  GaCl  and  GaCh  is 
relatively  small.  This  would  make  the  effective  value  of  n less  than  3 and  thus  make  the 


value  of  A°Gf(T)  less  negative. 
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Two  methods  were  employed  to  check  the  consistency  of  the  obtained  results.  In 
the  first  method,  termed  the  second  law  method,  the  standard  enthalpy  of  formation  is 
related  to  the  standard  Gibbs  energy  of  formation  by  the  following  equation: 

A°Hf=-Rd(-A‘’Gf/RT)/d(lA')  (2.7) 

Thus,  the  slope  of  -A°Gf/RT  vs  1/T  gives  -A"Hf/R.  For  a linear  temperature 

dependence  of  the  Gibbs  energy  function  the  value  of  - A°Hf  is  a constant.  The  slopes 

were  measured  graphically  for  each  measurement  temperature.  The  values  of  standard 
enthalpy  of  formation  of  GaAs  by  the  second  law  method  at  different  temperatures  are 
shown  in  Table  2-2  and  Figure  2-5.  It  is  important  to  note  that  the  second  law  analysis 
involves  slope  measurement,  which  is  subject  to  considerable  uncertainity.  In  the 

temperature  range  of  the  present  work  the  average  value  for  A“Hf  (Tav=  903  K)  is  -28.27 
kcal/mol.  The  value  obtained  from  the  second  law  method  after  extrapolation  of  the  - 
A°Gf /RT  vs  1/T  curve  to  298.15  K is  A°Hf(298.15  K)“  = -25.83  kcal/mol. 

The  second  consistency  test,  termed  the  third  law  method,  includes  high 
temperature  Cp  and  standard  entropies  (hence  the  third  law  method).  The  standard 

enthalpy  of  formation  -A°Hf  (298.15  K)  is  calculated  using  the  following  formalism: 
-A"Hf(298.15  K)  = -A°Gf(298.15  K)- A,  |°H(T)-°H(298.15  K)j  + 

r (2.8) 

TA°Sf(298.15  K)  + TAf{"S(T)-“S(298.15  K)} 

where  A,  indicates  thermal  function  difference  of  the  formation  reaction  (eqn  2-4). 

In  the  above  treatment,  pure  component  specific  heats  and  the  standard  elemental 
entropies  at  298.15  K were  taken  from  the  Janaf  Thermochemical  Tables  [Jan85]  and 
were  consistent  with  those  from  Hultgren  [Hul73].  The  standard  entropy  of  GaAs  at 
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298. 15K  was  taken  from  Piesbergen  [Pie63]  and  Holste  [Hol72].  The  heat  capacity  data 
for  GaAs  (s)  was  taken  from  the  optimized  correlation  by  Chatillon  et  al.  [Cha90]. 

It  is  interesting  to  see  the  temperature  variation  in  the  - A°Hf(298.15  K)  values 

obtained  from  the  third  law  analysis  as  shown  in  Figure  2-6  and  Table  2-2.  The  values 
show  systematic  temperature  dependence,  suggestive  of  systematic  errors.  The  error 
could  be  in  either  the  A°Gf  (T)  values  or  in  the  high  temperature  heat  capacity  difference 

data.  However,  the  temperature  dependence  fluctuates  and  is  acceptable  within 
experimental  errors. 

The  mean  value  of  A“Hf(298.15  K)  in  this  temperature  range  is  -19.94  kcal/mol 

which  agrees  excellently  with  the  measured  value  of  -20.96  kcal/mol  by  Yamaguchi  et  al. 
[Yam94]  using  tin  solution  calorimetry  and  that  of  -19.69  kcal/mol  measured  by 
Martosudirjo  and  Pratt  [Mar74]  using  precipitation  calorimetry.  It  also  compares  well 
with  the  calculated  value  of  -21.24  kcal/mol  reported  by  Chatillon  et  al.  [Cha90].  This  is 
comparable  to  the  experimental  errors  of  +0.12  kcal/mol  reported  by  Yamaguchi  et  al. 
[Yam94]  and  +0.11  kcal/mol  reported  by  Martosudirjo  and  Pratt  [Mar74].  It  also 
compares  well  with  the  calculated  value  of  -21.24  kcal/mol  reported  by  Chatillon  et  al. 
[Cha90].  The  value  dilfers  from  the  data  of  Thurmond  [Thu65],  Drowart  and  Goldfmger 
[Dro58],  and  Gutbier  [Gut61].  The  method  used  by  Sirota  [Sir58],  Abbasov  et  al. 
[Abb66]  and  Krestovnikov  [Kre71]  reported  emf  measurements  at  temperatures  below 
743K,  so  the  calculations  of  - A“Hf(298.15  K)  were  based  upon  extrapolations  of  their 
data.  On  the  other  hand,  the  data  by  Pupp  et  al.  [Pup74]  and  Arthur  [Art67]  were  too 
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low.  Thus,  these  data  were  discarded  for  comparison  purposes.  The  reported  values  from 
literature  have  errors  associated  with  errors  in  the  values  of  specific  heats,  Cp,  of  pure 
components.  The  errors  get  easily  propagated  due  to  the  temperature  dependency  of  the 
Cp  values.  This  is  also  seen  in  the  calculated  values  in  the  present  study.  The  temperature 
dependence  seen  in  the  present  study  also  includes  the  errors  associated  with  the  Cp  data. 
However,  a possible  reason  for  the  temperature  dependence  is  the  increase  in  As  loss  as 
the  temperature  is  increased. 

In  summary,  the  standard  Gibbs  energy  of  formation  of  GaAs  (s)  at  temperatures 
above  743  K was  measured  by  emf  method  and  was  checked  for  thermodynamic 

consistency  using  both  second  and  third  law  analysis  and  comparing  - A°Hf(298.15  K) 

values  with  the  experimental  measurements  of  Martosudirjo  and  Pratt  [Mar74],  and 
Yamaguchi  et  al.  [Yam94].  The  results  from  the  third  law  analysis  were  in  excellent 
agreement  with  the  experimental  data,  however,  the  discrepancy  between  the  measured 
and  calculated  values  indicates  random  errors  inherent  to  experimental  measurement,  and 
systematic  errors  inherent  to  the  choice  of  Cp(T)  data  used  in  the  calculations  by 
numerous  authors.  Considering  all  uncertainities  including  those  in  the  Gibbs  energy 
measurements  it  is  clear  that  the  present  calculations  are  in  excellent  agreement  with 
literature  values. 

In- As  System 

Phase  diagram  data.  Phase  diagram  data  for  In-As  has  been  studied  fairly 
extensively  by  several  authors  [Liu53,  Hal63,  Kar73,  PerSO,  Kno85  and  DeW86].  Most 
of  the  phase  diagram  data  is  available  for  In-rich  segent  of  the  phase  diagram.  Only  the 
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Figure  2-5.  Calculated  second  law  values  of  the  enthalpy  of  formation  of  GaAs  at 

298. 15K. 


35 


- - C^  <N  CS 

(|OUJ/|BO>j)  „|^5^'86Z'^HoV 


o 

00 

o> 


O 

CD 

O) 


O 

CN 

(J> 


O 

O 

O) 


O 

CO 

CO 


O 

CO 

00 


O 

00 


o 

CN 

00 


0 


CO 

0 

Q. 

E 

0 

K 


Uh 


igure  2-6.  Calculated  third  law  values  of  the  enthalpy  of  formation  of  GaAs  at  298. 1 5K. 
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data  by  Liu  and  Peretti  [Liu53]  and  Karataev  et  al.  [Kar73]  extend  to  the  As  rich-segment 
of  the  phase  diagram  (xas  > 0.70).  The  measurements  of  Hall  [Hal63]  and  Perea  and 
Fonstad  [PerSO]  were  made  by  heterogeneous  equilibrium  and  weight  loss  technique.  Liu 
and  Peretti  [Liu53]  used  differential  thermal  analysis  (DTA)  and  X-ray  analysis,  while 
Knobloch  and  Thoma  [Kno85]  and  DeWinter  and  Pollack  [DeW86]  used  source 
dissolution  method.  Karataev  et  al.  [Kar73]  measured  values  by  dew  point  and 
continuous  weighing  technique. 

Calorimetric  and  enthalpy  data.  The  enthalpy  of  formation  of  InAs  has  been 
determined  by  Yamaguchi  et  al.  [Yam88]  over  a wide  range  of  temperature  700  to  1500  K 
using  drop  calorimetry.  Standard  enthalpy  of  formation  of  InAs  at  298  K was  determined 
by  Schottky  and  Bever  [Sch58]  using  dissolution  calorimetry  and  by  Sharifov  et  al. 
[Sha63]  using  a calorimetric  bomb.  Enthalpy  of  fusion  data  is  available  from  Cox  and 
Pool  [Cox67]  and  Yamaguchi  et  al.  [Yam88]  at  their  measured  fusion  temperatures  of 
1213  K and  1221  K respectively  by  drop  calorimetry.  Heat  capacities  have  been 
determined  in  the  temperature  range  298  to  1200  K by  drop  calorimetry  [Cox67,  Lic69]. 
The  measured  heat  capacity  data  was  correlated  by  means  of  a least  square  fit  using  a 
polynomial  expansion  in  terms  of  temperature  in  the  range  160  to  1211  K [Cha90].  This 
correlation  was  used  to  perform  a third  law  analysis  on  the  measured  Gibbs  energy  data  in 
the  present  work.  The  pure  compound  standard  entropy  for  InAs  at  298.15  K was 
selected  as  20.84  cal/mol,  a value  recalculated  by  ChatiUon  et  al.  [Cha90]  from 
thermometric  measurements  of  Piesbergen  [Pie63],  Dash  et  al.  [Das74]  and  Holste 


[Hol72]. 
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Vapor  pressure  and  activity  data.  Partial  pressures  of  As2  and  As4  have  been 
measured  by  Knudsen  cell  mass  spectrometry  [Pup74,  Gut61,  Gut59  and  Gol59],  Total 
pressure  measurements  were  made  using  static  methods  [Boo57,  Fol58  and  Kar73]  such 
as  dew  point  technique,  and  by  dynamic  method  such  as  transport  of  As2  and  AS4  in  a H2 
flow  [Khu74,  Khu77], 

The  vapor  pressure  data  reported  for  In-As  system  show  a wide  distribution. 
Activities  of  In  in  the  two  phase  region  InAs-As  was  measured  by  electrochemical 
methods  [Kre71,  Abb64]  upto  a maximum  temperature  of  784  K.  More  recently  Tmar  et 
al.  [Tma87]  measured  As2  pressure  along  the  liquidus  curve  in  the  In-rich  segment  of  the 
phase  diagram  over  a temperature  range  of  846  to  1086  K.  Direct  measurement  of  In 
activity  has  been  made  in  the  single  liquid  phase  extending  all  the  way  to  the  equilibrium 
liquidus  curve  (liquid+InAs)  using  emf  technique  in  the  temperature  range  of  920  to  1 260 
K [Ito82]. 

In  this  study,  the  activity  of  indium  was  measured  along  the  In-rich  liquidus  in  the 
In-As  system  in  the  temperature  range  833  K to  1153  K using  a solid  state  galvanic  cell. 
Furthermore,  since  no  data  is  available  on  the  standard  Gibbs  energy  of  formation  of  InAs 
(s)  at  temperatures  above  784  K,  such  measurements  were  also  targeted.  An  assessment 
of  the  phase  diagram  and  thermodynamic  data  for  this  system  was  performed  by  Tmar  et 
al.  [Tma84].  In  their  preliminary  assessment  work,  liquid  solution  model  parameters  were 
estimated  from  an  optimization  of  other  data,  primarily  phase  diagram  data  and  arsenic 
partial  pressure  measurements.  As  discussed  in  the  previous  section  for  Ga-As  system, 
arsenic  activities  in  the  liquid  phase  calculated  from  arsenic  vapor  pressure  measurements 
is  subject  to  considerable  uncertainity  because  arsenic  sublimes  to  form  polyatomic  vapor 
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species  (primarily  As2  and  AS4)  at  atmospheric  pressure  which  have  a wide  range  of  scatter 
for  equilibrium  dissociation  constants.  In  general,  the  thermodynamics  of  liquid  As  is  not 
well  defined.  The  situation  is  aggravated  by  the  significant  scatter  in  the  partial  pressure 
measurements  that  are  likely  produced  by  parasitic  evaporation  and  chemical  side 
reactions. 

Owing  to  such  reasons,  the  liquid  solution  thermodynamic  behavior  has  not  been 
accurately  ascertained.  The  more  recent  assessment  [Cha90]  used  recently  measured 
liquid  phase  activity  of  Itoh  et  al.  [Ito82],  enthalpy  data  of  Glazov  et  al.  [Gla75]  and 
Yamaguchi  et  al.  [Yam88],  and  As2  vapor  pressure  measurements  of  Tmar  et  al.  [Tma87]. 
Predictions  have  suggested  moderate  negative  deviations  from  ideal  solution  behavior. 

In  these  experiments,  both  pure  As  (s)  and  As  (1)  can  be  used  as  the  standard 
reference  states  for  the  As2  and  AS4  vapor  pressure  data.  Selected  values  for 

A°Hf(298.15  K) , A"Sq^5(298.15  K)  and  Cp(T)  are  available  from  the  recent  assessment 
work  by  Tmar  et  al.  {Tma85j.  As  in  the  case  for  Ga-As  system,  Cp(T)  and  the  value  for 
°Sqjas(298.15  K)  for  In- As  system  were  recalculated  [Cha90]  using  experimental  Cp(T) 

data  of  Cox  and  Pool  [Cox67],  and  Lichter  and  Sommelet  [Lic69].  The  authors  [Cha90] 
also  recalculated  As  vapor  pressure  data  by  using  an  empirical  equation  for  P^^  and 

transforming  the  data  using  an  equation  for  dissociation  constant  Kp.  These  calculations 
were  referred  to  pure  As  (s)  in  the  standard  state.  However,  even  after  such  corrections 
were  applied  to  the  vapor  pressure  data,  significant  disagreement  still  existed. 
Furthermore,  the  emf  measurements  [Kre71,  Sir67,  Abb66]  in  literature  were  performed 
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Table  2-1.  Experimental  Gibbs  energy  of  formation  of  GaAs  (s). 


Temperature  (K) 

EMF  (mV) 

A“Gf  (kJ/g-at) 

833 

232.5 

-33.65 

853 

224.0 

-32.42 

873 

220.0 

-31.85 

892 

214.0 

-30.97 

913 

207.5 

-30.04 

933 

203.5 

-29.46 

953 

201.0 

-29.10 

973 

198.5 

-28.74 

Table  2-2.  Enthalpy  of  formation  of  GaAs  (s)  by  second  law  analysis  and  the  values  at 

298. 15K  by  the  third  law  method. 


Temperature  (K) 

A°Gf  (kJ/g-at) 

A°Hf(T)“ 

(kcal/mol) 

A”Hf(298.15K)“ 

(kcal/mol) 

833 

-33.65 

-40.93 

-20.80 

853 

-32.42 

-23.05 

-20.35 

873 

-31.85 

-30.60 

-20.23 

892 

-30.97 

-33.78 

-19.95 

913 

-30.04 

-29.01 

-19.66 

933 

-29.46 

-24.24 

-19.53 

953 

-29.10 

-21.86 

-19.50 

973 

-28.74 

-22.65 

-19.47 
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at  low  temperatures  and  were  found  inconsistent  with  the  phase  diagram  data  [Tma85, 
Cha90].  Keeping  in  mind  the  uncertainity  in  the  available  thermodynamic  data,  especially 
the  Gibbs  energy  data,  experiments  in  this  study  were  designed  to  obtain  liquid  phase  In 
activity  and  high  temperature  Gibbs  energy  of  formation  of  InAs  (s). 

Activity  of  In  Along  the  In-rich  Liquidus 

The  experimental  cell  used  for  the  activity  measurements  along  the  In-rich  liquidus 
line  is  schematically  shown  below: 

6203(1),  In203(s),  In  (1,  pure)/AfSZ//InxAsi.x  (1),  InAs  (s),  In203(s),  6203(1)  (V) 

In  (1,  pure)  + 3/2  O'^  = Iu203  -i-  3e'  (2.9) 

Iu203  + 3e’  = In  (in  InxAsi.x  (1))+  3/2  O'^  (2. 10) 


In  (1,  pure)  = In  (in  InxAsi.x  (1)) 


(2.11) 


The  In  activity  was  calculated  from  the  following  expression: 


a^  =exp 


-3FE 

RT 


/ 


(2.12) 


where  F = Faraday  Constant,  96,485  C/eq-mV 
and  E = Open  circuit  emf,  mV 

The  results  for  the  In  activities  along  the  liquidus  are  shown  in  Table  2-3  and 
Figure  2-7.  In  this  experiment  the  working  electrode  contained  the  two  phase  mixture  of 
InAs  solid  and  In-rich  liquid  aUoy.  As  the  temperature  was  increased,  some  of  the  solid 
InAs  dissolved  in  the  In-rich  liquid,  decreasing  the  In  composition  of  the  equilibrium  liquid 
solution.  According  to  the  phase  diagram  data,  the  composition  Xi„,  decreases  in  a nearly 
exponential  manner  with  temperature  for  In-rich  compositions.  The  measured  emf  values 
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are  expected  to  exhibit  an  exponential  dependence  and  such  is  observed  in  the  results. 
The  results  are  in  good  agreement  with  the  liquid  phase  In  activity  measurements 
conducted  by  Itoh  [Ito82]  using  an  electrochemical  technique.  However,  the  present  set 
of  data  represent  a stronger  negative  deviation  from  ideality  than  that  suggested  by  Itoh. 
Standard  Gibbs  Energy  of  Formation  of  InAs  (s) 

The  cell  that  was  used  in  the  experiment  for  the  measurement  of  standard 
Gibbs  energy  of  formation  of  InAs  (s)  alongwith  the  half  cell  reactions  is  shown  below: 
6203(1),  In  (1,  pure),  InaOs  (s)//YSZJ/lnAs  (s).  As  (s,  pure),  Iu203  (s),  B2O3  (1),  (VI) 

In  (1,  pure)  + 3/2  O'^  = 1/2  In203  (s)  + 3e-  (2. 1 3) 

1/2  In203  (s)  + As  (s,  pure)  + 3e‘  = InAs  (s)  + 3/2  O'^  (2. 14) 


In  (1,  pure)  + As  (s,  pure)  = InAs  (s)  (2. 15) 

The  above  equation  is  transformed  into  the  following  overall  cell  equation  by  applying 
Wagner's  approximation  to  the  chemical  potential  of  the  essentially  stoichiometric 
compound  InAs  (s): 

In  (1,  pure)  = In  (ss)  (2. 1 6) 

where  ss  denotes  the  solid  solution  In- As  with  composition  fixed  to  lie  along  line  towards 
the  InAs(s)  + As(s)  region.  This  is  equivalent  to  measuring  the  chemical  potential  of  In 
along  the  As-rich  solidus  of  the  compound  InAs.  The  standard  Gibbs  energy  of  formation 
of  InAs  referred  to  In  (1)  and  As  (s)  states  is  given  by  the  following  relation: 

A°G,(InAs)=°G.„^3,-“G^,-°G^3,,,„  = -3FE  (2. 17) 

The  emf  values,  and  the  calculated  standard  Gibbs  energy  of  formation  of  InAs  are 
tabulated  in  Table  2-4  and  can  be  seen  in  Figure  2-8  to  exhibit  a linear  temperature 
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dependence.  The  A°Gf  values  compare  well  with  the  earlier  emf  measurements  [Abb64 

and  Kre71]  when  extrapolated  to  the  low  temperature  range.  Some  variation  in  the 
temperature  dependence  of  the  present  data  is  indicative  of  possible  As  loss  by 
volatilization  (especially  at  higher  temperatures)  thereby  resulting  in  lowering  of  the  open 
circuit  emf  values.  The  values  in  the  low  temperature  region  have  uncertainties  related  to 
difficulties  in  achieving  equilibrium  owing  to  slow  diffusion  and  non-homogeneous  two- 
phase  solid  mixture.  Steady  values  were  recorded  at  each  temperature  as  being 
representative  of  the  open  circuit  emf  for  the  equilibrium  cell  at  those  temperatures. 

As  in  the  previous  section  for  Ga-As  system,  both  second  and  third  law  methods 
were  employed  to  check  the  consistency  of  the  Gibbs  energy  data  for  In-As  system. 
Using  the  second  law  method,  the  standard  enthalpy  of  formation  was  obtained  by 

calculating  the  slope  of  -A°Gf/RT  vs  1/T  curve  which  gave  the  value  of -A°Hf/R  at  the 
temperature  of  interest  using  the  following  expression: 

A"Hf  = - R d(-  A°Gf  /RT)/d(  1/T)  (same  as  equation  2.7) 

The  values  of  standard  enthalpy  of  formation  of  InAs  obtained  from  the  second 
law  analysis  at  different  temperatures  is  shown  in  Table  2-5  and  Figure  2-9.  In  the 

experimental  range  of  temperature,  the  average  value  of  A°Hf(Tav  = 878  K)  was  -12.83 
kcal/mol.  The  value  obtained  using  the  second  law  method  after  extrapolation  of  the  - 
A°Gf/RT  vs  1/T  curve  to  298.15  K was  A“Hf  (298.15  K)“  = -11.92  kcal/mol.  Owing  to 

extrapolation  of  experimental  data  and  errors  involved  in  slope  measurement,  the  second 
law  value  is  not  as  accurate  as  the  value  obtained  from  the  third  law  analysis  as  described 


next. 
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Using  pure  component  and  solid  compound  Cp  data  and  standard  entropies,  the 
standard  enthalpy  of  formation  -A“Hf  (298.15  K)  is  determined  by  the  third  law  method 
using  the  following  formalism: 

-A°Hf(298.15  K)  = -A‘’Gf(298.15  K)- A^  {°H(TKH(298.15  K)}  + 

TA°Sf  (298.15  K)  + TAf  |‘’S(T)-”S(298. 15  K)|  (same  as  equation  2.8) 

where  Af  indicates  thermal  function  difference  of  the  formation  reaction. 

It  is  interesting  to  see  the  temperature  variation  in  the  A°Hj(298.15  K)  values 

obtained  from  the  third  law  analysis  as  shown  in  Figure  2-10  and  Table  2-5.  The  values 
show  non-systematic  temperature  dependence,  suggestive  of  random  errors  typical  in  the 
experiment.  The  temperature  dependence  varies  and  is  acceptable  within  experimental 

errors.  The  mean  value  of  A°Hf(298.15  K)  in  this  temperature  range  is  -14.53  kcal/mol. 

This  compares  excellently  with  the  calorimetric  measured  values  of  -14.34  kcal/mol  by 
Yamaguchi  et  al.  [Yam94],  -13.80  kcal/mol  by  Sharifov  et  al.  [Sha63],  and  the  calculated 
value  of  -14.05  kcal/mol  reported  by  Chatillon  et  al.  [Cha90].  This  is  also  in  agreement 
with  the  reported  value  of  -14.80  kcal/mol  by  Schottky  and  Bever  [Sch58].  The 

deviations  in  the  calculated  values  of  A"Hf(298.15  K)  are  of  the  order  of  +0.17  kcal/mol 

when  compared  to  the  mean  value  of  -14.53  kcal/mol.  This  is  comparable  to  the 
experimental  errors  of  +0.24  kcal/mol  reported  by  Yamaguchi  et  al.  [Yam94]  and  ±0.19 
kcal/mol  reported  by  Sharifov  et  al.  [Sha63]. 
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Figure  2-10.  Calculated  third  law  values  of  the  enthalpy  of  formation  of  InAs  at  298.15  K. 
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The  experimental  data  exhibits  tight  distribution  when  compared  to  the  deviations 
in  the  calculated  values  by  Tmar  and  Chatillon  [Tma87]  and  the  measurements  of  Schottky 
and  Bever  [Sch58].  The  value  differs  from  the  data  of  Pupp  et  al.  [Pup74],  Gutbier 
[Gut61]  and  Goldfinger  and  Jeunehomme  [Go  159].  The  emf  method  used  by  Abbasov  et 
al.  [Abb66]  and  Krestovnikov  [Kre71]  reported  measurements  at  temperatures  below  784 

K,  so  the  calculations  of  -A°Hf(298.15  K)  were  based  upon  extrapolations  of  their  data. 

Thus,  these  data  were  discarded  for  comparison  purposes.  The  reported  values  from 
literature  have  errors  associated  with  errors  in  the  values  of  specific  heats  of  pure 
components,  Cp(T).  The  errors  get  easily  propagated  due  to  the  temperature  dependency 
of  the  Cp(T)  values.  This  is  also  seen  in  the  calculated  values  in  the  present  study.  The 
temperature  dependence  seen  in  the  present  study  also  includes  the  errors  associated  with 
the  Cp(T)  data.  A possible  cause  exists  in  the  temperature  dependence  of  Gibbs  energy 
data  due  to  As  losses  due  to  volatilization. 

In  the  above  calculations,  pure  component  specific  heats  and  the  standard 
elemental  entropies  at  298.15  K were  taken  from  the  Janaf  Thermochemical  Tables 
[Jan85]  and  were  consistent  with  those  from  Hultgren  [Hul73].  The  standard  entropy  of 
InAs  at  298.15  K was  taken  from  Piesbergen  [Pie63].  The  results  from  analyses  were  in 
excellent  agreement  with  the  experimental  data  however;  the  discrepancy  between  the 
measured  and  calculated  values  indicates  random  errors  inherent  to  experimental 
measurement,  and  systematic  errors  inherent  to  the  choice  of  Cp  data  used  in  the 
calculations  by  numerous  authors.  Considering  all  uncertainities  including  those  in  the 
Gibbs  energy  measurements  it  is  clear  that  the  present  calculations  are  in  excellent 


agreement  with  literature  values. 
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Sb-Te  System 

Recent  investigations  in  the  field  of  erasable  optical  recording  technology  have 
shown  that  materials  which  undergo  a reversible  phase  transformation  between  the 
crystalline  and  glassy  states  upon  irradiation  are  promising  alternatives  to  magneto-optical 
materials.  Many  of  these  alternative  recording  media  are  based  on  chalcogenide  alloys, 
most  notably  tellurium-containing  alloys  [Mar90].  In  particular,  the  antimony-tellurium 
system  is  an  attractive  binary  constituent  for  multicomponent  optical  storage  media  [Uij85, 
Poe85,  Cha90  and  Yam91].  Since  the  glassy  state  is  obtained  by  quenching  from  the 
liquid  state,  insight  into  the  thermodynamic  properties  and  degree  of  ordering  in  the  glassy 
state  can  be  gained  from  study  of  the  thermodynamic  behavior  of  the  equilibrium  liquid 
state. 

A thorough  review  of  the  thermodynamics  of  the  Sb-Te  system  has  been  reported 
[Gho89].  The  solid-liquid  Sb-Te  phase  diagram  is  dominated  by  the  fairly  well 
characterized  congruently  melting  compound  SbaTea.  The  assessment  results  of  Ghosh  et 
al.  [Gho89]  also  included  two  solid  solutions  in  the  Sb-rich  segment  of  the  phase  diagram, 
though  no  experimental  determinations  of  the  solution  thermodynamics  for  these  two 
phases  have  been  reported.  The  only  theirnodynamic  property  measurements  for  the 
liquid  solution  are  two  calorimetric  determinations  of  the  enthalpy  of  mixing  throughout 
the  entire  composition  range  at  a temperature  near  the  melting  point  of  Sb2Te3  [Mae72, 
Feu88].  These  results  showed  an  exothermic  enthalpy  change  upon  mixing  with  a 
minimum  value  reported  at  a composition  close  to  the  stoichiometry  of  the 
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Table  2-3.  Summary  of  emf  measurements  for  In  activity  along  liquidus. 


Temperature 

(K) 

EMF  (mV) 

am 

X*i„ 

r*m 

833 

0.70 

0.971 

0.970 

1.003 

883 

1.53 

0.941 

0.945 

0.996 

918 

3.30 

0.882 

0.927 

0.952 

953 

5.90 

0.806 

0.908 

0.888 

994 

6.80 

0.788 

0.876 

0.899 

1033 

9.20 

0.733 

0.843 

0.870 

1073 

15.40 

0.607 

0.803 

0.756 

1108 

19.50 

0.542 

0.763 

0.710 

1128 

21.0 

0.523 

0.736 

0.711 

1153 

27.5 

0.436 

0.703 

0.620 

* Using  Optimized  Phase  Diagram  Data  from  [Tma84] 
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Table  2-4.  Experimental  Gibbs  energy  of  formation  of  InAs  (s). 


Temperature  (K) 

EMF  (mV) 

A“Gf  (kJ/g-at) 

803 

120.5 

-17.44 

833 

118.5 

-17.15 

853 

117.0 

-16.94 

867 

115.0 

-16.65 

873 

114.3 

-16.55 

891 

114.0 

-16.50 

913 

112.9 

-16.34 

921 

112.6 

-16.30 

953 

108.5 

-15.70 
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Table  2-5.  Standard  enthalpy  of  formation  of  InAs  (s)  by  the  second  law  analysis  and 

values  at  298. 15K  by  the  third  law  method. 


Temperature  (K) 

A°G,  (kJ/g-at) 

A"Hf(T)" 

(kcal/mol) 

A°Hf(298.15  K)“ 
(kcal/mol) 

803 

-17.44 

-12.12 

-13.94 

833 

-17.15 

-12.12 

-14.02 

853 

-16.94 

-13.75 

-14.06 

867 

-16.65 

-15.90 

-14.00 

873 

-16.55 

-11.52 

-14.03 

891 

-16.50 

-9.94 

-14.13 

913 

-16.34 

-10.33 

-14.22 

921 

-16.30 

-13.91 

-14.26 

953 

-15.70 

-15.90 

-14.20 
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compound  Sb2Tc3.  In  this  study,  the  first  experimental  measurement  of  activity  of  the 
antimony  in  Sb-Te  melts  was  performed. 

Activity  of  Sb  in  Sb-Te  Melts 

The  antimony  activity  was  measured  using  a solid-state  galvanic  cell  at  a 
temperature  near  9 1 1 K.  The  oxygen-ion  conducting  electrolyte  yttria  stabilized  zirconia 
(YSZ)  was  selected  in  preference  to  a molten  salt  electrolyte  because  of  the  low  boiling 
temperatures  of  the  stable  halides  (e.g.,  494  K for  SbCb).  With  this  selection  of  the 
electrolyte,  the  electrochemical  ceU  used  to  measure  the  Sb  activity  is  schematically 
represented  as: 

W/Sb,  SbzTej/A^SZ//  SbxTe,.„  SbzOjAV  (VII) 

The  reference  electrode  consisted  of  pure  liquid  antimony  in  equilibrium  with  its 
stable  oxide,  Sb203.  Assuming  the  solubility  of  oxygen  in  the  SbxTei.x  liquid  phase  is 
negligible,  the  activity  of  antimony,  asb,  is  related  to  the  measured  open  circuit  potential, 
E,  by  the  following  expression: 

RTlnasb  = RTln  = -3FE  (2. 1 3) 

where  T is  the  absolute  temperature,  R is  the  gas  constant,  F is  the  Faraday  constant,  and 
Pq  “ is  the  equilibrium  oxygen  partial  pressure  in  the  reference  (I)  or  the  working 
electrode  (II). 

The  temperature  of  measurement  was  about  911  K;  a temperature  between  the 
melting  temperatures  of  pure  Sb  (904  K)  or  Sb2Te3  (891  K)  and  that  of  antimony 
sesquioxide  (930  K).  At  this  temperature,  there  was  concern  that  volatilization  losses 
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would  affect  the  working  electrode  composition.  For  this  reason  both  cell  constructions 
depicted  in  Figures  2-2  and  2-3,  designed  to  minimize  electrode  material  evaporation, 
were  used. 

Measurements  using  the  encapsulated  cell  assembly  design  was  also  used  to 
conduct  experiments  on  selected  alloy  compositions  to  confirm  the  findings  of  the 
compression  cell  assembly.  The  working  electrode  was  comprised  of  a mixture  of  SbxTci. 
X and  SbaOa  Sb203  powder  while  the  reference  electrode  contained  a mixture  of  pure  Sb 
and  Sb203  powder.  The  annular  space  between  the  YSZ  crucible  and  alumina  cell  holder 
contained  the  reference  electrode  material  (Sb  + Sb203),  while  the  YSZ  crucible  contained 
the  working  electrode,  which  consisted  of  the  alloy  (SbxTei.x)  of  the  desired  composition 
and  the  most  stable  oxide  (Sb203). 

In  these  experiments  B2O3  was  not  used.  Careful  design  and  assembly  of  these 
electrodes  were  able  to  yield  stable  emfs,  however,  the  incorporation  of  B2O3  proved 
successful  in  reproducing  these  results  in  later  experiments. 

Nineteen  alloys  were  studied  throughout  the  entire  composition  range  at  a 
temperature  near  9 1 1 K.  The  antimony  activities  were  calculated  from  the  measured  emf 
and  temperature  with  the  use  of  equation  (2.13)  and  are  shown  in  Figure  2-11  as  a 
function  of  composition.  The  results  indicate  that  the  solution  is  nearly  ideal  for  Sb-rich 
compositions  (Nxe  < 0.45).  For  Te-rich  compositions,  however,  moderate  negative 
deviations  from  ideal  solution  behavior  are  evident. 
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Neglecting  the  small  variation  in  activity  in  the  measurement  temperature  range 
(908.6  to  915.7  K),  the  corresponding  activity  coefficient  of  tellurium  in  the  melt  was 
determined  by  integration  of  the  Gibbs-Duhem  relation  using  Darken’s  method  [Dar53]; 

Np 

f Te  ~ ~^Sb  Nre  Ngb  ~ J®Sb*^^Te  (2.14) 

n;=. 

In  y 

where  and  Nj  is  the  mole  fraction  of  species,  i.  A summary  of  the 

(1  ~ Ngb ) 

experimental  results  and  the  calculated  component  activities,  partial  molar  Gibbs  energies 
Agj,  and  Gibbs  energy  of  mixing  (AG^^)  is  given  in  Table  2-6.  Smoothed  values  of 

AGnj,  were  combined  with  the  calorimetric  measurements  of  Feutelais  et  al.  [Feu88]  to 
yield  the  entropy  of  mixing,  AS^^ , as  displayed  in  Table  2-7  and  Figure  2-12. 

Discussion.  There  is  considerable  experimental  evidence  for  the  formation  of  short 
range  order  in  Sb-Te  melt.  The  enthalpy  of  mixing  has  been  measured  throughout  the 
entire  composition  range  in  two  different  laboratories  [Mae72,  Feu88].  Both  of  the 
measurements  show  an  exothermic  enthalpy  change  with  a minimum  value  near  the 
composition  of  SbaTcs.  In  addition,  the  electrical  conductivity  and  viscosity  data  of 
Glazov  et  al.  [Gla69]  display  an  extremum  near  Nxe=  0.55.  Furthermore,  Abrikosov  et  al. 
[Abr59]  pointed  out  that  the  thermoelectric  power  shows  a sharp  increase  at  the  same 
composition.  These  results  suggest  that  the  melt  structure  has  short  range  order,  with  the 
extent  of  ordering  a maximum  at  a composition  close  to  the  stoichiometry  of  the  highest 
melting  compound,  Sb2Tc3. 

The  nearly  ideal  solution  behavior  for  Sb-rich  compositions  is  apparently  a result 
of  the  measurement  temperature  (-911  K)  being  sufficiently  greater  than  the  liquidus 
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Figure  2-12.  Calculated  entropy  of  mixing  values  for  Sb-Te  alloys  from  calorimetric 

measurements  and  smoothed  Gibbs  energy  of  formation  values. 
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temperature  (minimum  value  near  813K  at  Nje  = 0.31  [Abr59,  Bor86])  to  dissociate  the 
clusters.  As  the  measurement  temperature  becomes  closer  to  the  liquidus  value,  the 
antimony  activity  shows  negative  deviations  that  are  consistent  with  an  increasing  extent 
of  unlike  atom  clustering.  The  Sb-Te  system,  however,  does  not  show  as  high  a tendency 
towards  association  as  does  the  Sn-Te  and  Pb-Te  systems  [Kot88]. 

Figure  2-1 1 also  shows  the  antimony  activity  proposed  by  Ghosh  et  al.  [Gho89]  as 
a result  of  their  assessment  study.  In  this  work,  a molten  solution  model  with  partial 
ionization  developed  by  Hillert  et  al.  [Hil85]  was  used  to  describe  the  melt  behavior. 
Other  researchers  have  used  chemical  theories  in  which  molecular-like  species  are 
postulated  to  describe  the  apparent  short  range  ordering  between  unlike  atoms  [Jor70, 
Cla82  and  Som82].  The  assessed  antimony  activity  of  Ghosh  et  al.  [Gho89]  shows 
excellent  agreement  with  the  measured  values  of  this  study  for  0 < Nxe  < 0.4,  while  slightly 
less  negative  deviations  were  predicted  for  more  Te-rich  compositions.  In  this  early 
assessment,  these  authors  used  the  only  available  enthalpy  of  mixing  data  of  Maekawa  et 
al.  [Mae72].  These  data  were  measured  at  a higher  temperature  (935  K)  than  the  more 
recent  data  of  Feutelais  et  al.  [Feu88]  for  which  the  temperature  ranged  from  909  to  918 
K.  Since  the  enthalpy  of  mixing  is  temperature  dependent  and  more  exothermic  at  lower 
temperatures,  the  solution  behavior  predicted  from  the  assessment  is  expected  to  be 
slightly  more  ideal. 

There  are  three  reported  measurements  in  the  literature  of  the  Gibbs  energy  of 
formation  of  the  compound  Sb2Te3  with  a galvanic  cell  technique  [Ger61,  Sem69  and 
Mal86].  These  results  can  be  compared  for  consistency  with  the  smoothed  value  of 


Table  2-6.  Summary  of  experimental  results  and  calculated  thermodynamic  quantities  for  Sb-Te  system. 
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Table  2-7.  Summary  of  smoothed  thermodynamic  mixing  values  for  Sb-Te  alloys  at  91 1 K. 


Nxe 

AGmix 

( J mol'  ^ ) 
[this  work] 

AHnrix 

(J  mol'^) 
[Feu88] 

. TASn^x 
(J  moT^) 
T = 911  K 

0.000 

0 

0 

0 

0.050 

-1885 

-1057 

828 

0.100 

-3503 

-2120 

1328 

0.150 

-4839 

-3157 

1682 

0.200 

-5938 

-4194 

1744 

0.250 

-6864 

-5264 

1600 

0.300 

-7671 

-6376 

1295 

0.350 

-8387 

-7498 

889 

0.400 

-9015 

-8563 

452 

0.450 

-9532 

-9473 

59 

0.500 

-9894 

-10124 

-230 

0.550 

-10056 

-10423 

-367 

0.600 

-9975 

-10306 

-331 

0.650 

-9626 

-9759 

-133 

0.700 

-9004 

-8822 

182 

0.750 

-8129 

-7589 

540 

0.800 

-7030 

-6188 

842 

0.850 

-5734 

-4747 

987 

0.900 

-4224 

-3332 

892 

0.950 

-2399 

-1862 

537 

1.000 

0 

0 

0 

Table  2-8.  Values  of  the  standard  Gibbs  energy  of  formation  of  Sbo.4Teo,6  at  540  K and 

673  K (kJ/mole). 


Tn,  (K) 

[Mal86] 

[Sem69] 

[Ger61] 

This  Work 

540 

-13.63 

-12.90 

-12.28 

673 

-12.74 

— 

-12.34 

-12.36 
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AGmij^  of  this  study  at  the  Nje  = 0.6  composition  by  considering  the  following 
thermodynamic  sequence: 


0.4  Sb  (Tn,,  s)  ^ 0.4  Sb  (91 1 K,  1 ) AG,  (2.15) 

0.6  Te  (T„„  s)  ^ 0.6  Te  (9 1 1 K,  1 ) AGj  (2.16) 

0.4  Sb  (91 1 K,  1 ) + 0.6  Te  (91 1 K,  1 ) -»  Sbo.4Teo.6  (91 1 K,  1 ) AG^x  (2. 17) 


Sbo.4Teo.6(91 1 K,  1 ) — > Sbo.4Teo.6(Tni,  s)  AG3  (2.18) 

The  overall  reaction  in  this  sequence  is  simply  the  formation  reaction  at  the 

galvanic  cell  measurement  temperature,  Tm.  Values  of  A°Gj(T^)  were  calculated  by 
combining  the  selected  values  of  [Bar77]  for  AG, , AGj  and  AG3  with  the  value  AG,^^  = 
- 9.98  kJ/mole  suggested  by  this  study.  The  results  of  this  calculation  are  compared  to  the 
galvanic  cell  measurements  of  A“Gf  at  two  selected  temperatures  in  Table  2-6.  This 

calculation  indicates  good  agreement  between  the  results  of  the  solid  and  the  liquid 
solution  emf  measurements.  Using  this  sequence,  the  calculated  value  of 

A“Gf(298.15  K)  for  Sbo.4Teo.6is  - 12.05  kJ/mole. 


Experimental  Results  in  Ternary  Systems 

Data  in  ternary  system  Ga-In-As  are  scarce  due  to  high  vapor  pressures  associated 
with  these  alloys  owing  to  the  As  component  in  the  alloys.  Furthermore,  solid  phase 
diffusivities  are  slow  creating  problems  in  the  attainment  of  thermodynamic  equilibrium. 
Some  phase  diagram  data  have  been  determined,  primarily  in  the  Ga-,  In-  rich  composition 
ranges.  Very  scarce  data  however,  exists  on  the  thermodynamic  properties  of  these  alloys. 
Calorimetric  measurements  of  enthalpies  of  formation  and  mixing  were  reported  by  Rugg 
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et  al.  [Rug95].  In  this  study,  activity  measurements  have  been  performed  in  both  liquid  and 
solid  solutions  in  the  Ga-In-As  system. 

Activity  of  Ga  in  In-rich  Ternary  Liquid  GavInyAs i ...v  Solutions 

The  activity  of  gallium  in  the  ternary  liquid  solution  with  indium  and  arsenic  was 
determined  at  two  overall  In-rich  compositions. 

Most  of  the  data  exists  in  terms  of  phase  diagram  determinations  [Van63,  Hoc66, 
Woo58,  Shi56]  or  calculations  [CreSl,  Lit91,  Str69]  in  the  pseudobinary  phase  region.  It 
was  essential  to  perform  direct  measurements  for  activity  data  in  this  system  so  that  the 
validity  of  solution  models  used  to  describe  this  system  could  be  judged.  The  experiments 
were  designed  under  the  constraints  of  high  In  compositions  and  low  Ga  and  As  contents. 
This  was  because  at  Ga-rich  compositions  the  measurements  were  least  accurate  owing  to 
the  extremely  low  emf  values  generated;  while  at  high  arsenic  compositions,  the  high 
arsenic  vapor  pressures  caused  experimental  problems. 

The  following  cell  was  used  in  the  experiment: 

Ga203,  Ga(l,  pure)/A^SZ//GaxInyAS].x-y  (1),  Ga203  (VIII) 

Ga  (1,  pure)  + 3/2  O'^  = Ga203  -i-  3e‘  (2. 1 9) 

Ga203  -I-  3e’  = Ga  (in  GaxInyAsi.x.y  (1))  + 3/2  O'^  (2.20) 


Ga  (1,  pure)  = Ga  (in  GaxInyAsi.x.y  (1)) 


(2.21) 


The  activity  of  gallium  in  the  ternary  alloy  is  given  as: 


Soa  = exp 


-3FE 

RT 


/ 


(2.22) 
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The  working  electrode  was  designed  by  using  measured  amounts  of  InAs(s), 
GaAs(s)  and  pure  In  metal,  along  with  some  quantity  of  Ga203  powder.  The  amounts  of 
the  materials  were  chosen  to  yield  two  alloys  of  different  compositions  in  the  In-rich 
region  of  the  ternary  phase  diagram.  The  amount  of  Ga203  in  the  electrodes  was 
approximately  0.8g.  The  following  are  the  compositions  studied: 

Ino.8465Gao.o6788Aso.o8563  (Cell  I) 

Ino.8253  Gao.0756I  ASo.09911  (Cell  II) 

With  the  increase  in  cell  temperature,  the  indium  metal  melts  first  (m.p.  ~ 428  K). 
The  solids  InAs  and  GaAs  start  dissolving  in  the  liquid  indium  melt  to  form  an  In-rich 
ternary  liquid  phase  GaxIUyAsi.x.y  which  is  in  equilibrium  with  the  pseudobinary  solid  phase 
GaxIni.xAs  where  'X'  denotes  mole  fraction  of  GaAs  in  the  pseudobinary  solid  phase.  This 
was  supported  by  the  observation  that  emf  vs  time  stayed  constant  (within  +/-  0.05  mV  for 
2 hrs  or  more)  at  different  temperature  settings  indicating  achievement  of  cell  equilibrium. 
As  the  temperature  is  increased,  the  In-rich  ternary  liquid  phase  GaxIUyAsi.x.y  gets  leaner  in 
the  indium  component  while  moving  along  the  lines  of  constant  In/As  ratio  on  the  liquidus 
surface  of  the  ternary  system.  Once  the  terminal  point  is  reached  on  the  liquidus  surface 
corresponding  to  the  prepared  alloy  composition  in  the  ternary  phase  diagram,  the  activity 
measurements  correspond  to  the  single  ternary  liquid  phase  GUxIUyAsi.x.y.  The  transition 
from  the  two-phase  experimental  measurement  to  that  in  the  single  phase  ternary  liquid 
phase  gives  the  coordinates  of  the  liquidus  point. 
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Table  2-9.  Summary  of  emf  measurements  in  the  ternary  system  along  lines  of  constant 

In/As  ratio. 


Temperature 

(K) 

EMF  (mV) 
Cell  r 

EMF  (mV) 
Cell  II" 

^Ga 

Tea 

aL 

r" 

/ Ga 

867 

••MM 

200.0 

^ ^ 

... 

0.000325 

917 

MMM 

178.0 

MMM 

MMM 

0.001156 

963 

248.0 

151.0 

0.000127 

MMM 

0.00424 

998 

150.0 

140.0 

0.005334 

0.079 

0.00756 

1025 

140.0 

113.0 

0.00858 

0.126 

0.02150 

0.280 

1056 

136.5 

110.0 

0.01107 

0.163 

0.02650 

0.350 

1081 

134.0 

111.0 

0.01331 

0.196 

0.02796 

0.370 

1116 

122.9 

102.0 

0.0215 

0.317 

0.04136 

0.547 

1153 

116.0 

100.4 

V 

0.03013  ‘ 

0.443 

0.04825 

0.638 

1189 

107.0 

99.6 

0.04356 

0.642 

0.05410 

0.714 

1233 

94.5 

88.0 

0.06934 

=1.00 

0.08330 

=1.00 

[ (Cell  I)  xca=  0.068,  Xi„=  0.847 
**  (Cell  II)  xoa  = 0.076,  x^  = 0.825 
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The  emf  measurements  for  the  two  cells  in  the  ternary  Ga-In-As  system  are 
presented  in  Table  2-9  and  shown  in  Figure  2-13.  The  emf  versus  temperature  curves  for 
the  two  cells  show  a nearly  linear  dependence.  Both  the  cells,  however,  show  a 
discontinuity  in  the  emf  versus  temperature  dependence  at  certain  point  in  which  the  slope 
of  the  curve  changes  drastically  indicating  a transition  from  a solid-liquid  two-phase  region 
to  a single  liquid  phase.  The  temperature  at  which  this  occurs  is  the  liquidus  temperature 
of  the  ternary  alloy  for  that  particular  composition.  The  liquidus  temperature  obtained 
graphically  for  the  alloy  in  Cell  I (xoa=  0.068,  Xi„=  0.847)  is  1020  K,  and  that  for  alloy  in 
Cell  II  (xoa  = 0.076,  Xi„  = 0.825)  is  1040  K.  These  values  are  in  excellent  agreement 
(within  10  K)  extrapolated  estimates  at  these  compositions  from  experimental  phase 
diagram  data  of  [Van63,  Hoc66,  Woo58,'  and  Shi56]  and  the  calculations  of  [Str72].  The 
activity  values  as  seen  from  Table  2-9  indicate  negative  deviations  from  ideal  solution 
behavior  in  the  ternary  liquid  phase  GaxInyAsi.x.y,  the  degree  of  non-ideality  decreasing 
with  increasing  temperature. 

Activity  of  GaAs  (s)  in  the  Pseudobinarv  Section  GaAs-InAs 

Very  little  information  is  available  on  the  thermodynamics  of  the  solid  solution 
GaxIni-xAs.  There  is  a need  to  provide  insight  into  the  thermodynamic  behavior  of  the 
compounds  GaAs  and  InAs  in  the  pseudobinary  section  GaAs-InAs  of  the  ternary  system, 

. and  to  understand  the  interactions  between  the  Ga  and  In  atoms  on  the  metallic  sublattice. 
Such  information  makes  it  possible  to  examine  the  stability  of  the  solid  solution  at  different 


67 


temperatures  and  regions  of  the  phase  diagram,  thereby,  relating  to  the  difficulty  or  ease  of 
growing  alloys  from  the  melt  by  standard  epitaxial  growth  techniques.  The  calculation  of 
structural  properties  and  phase  diagram  aids  in  evaluating  these  alloys  for  use  in  long 
wavelength  infrared  material  (LWIR). 

Electrochemical  experiments  were  performed  by  Krestovnikov  et  al.  [Kre73] 
allowing  calculation  of  activities  of  GaAs  and  InAs  in  the  solid  solution,  which  they 
described  as  being  ideal.  Later  studies  on  low  temperature  heat  capacity  and  mass 
spectrometric  vapor  pressure  measurements  by  Tmar  et  al.  [Tma88]  suggest  immiscibility 
in  the  GaxIni.xAs  system  with  a critical  temperature  of  about  700  K.  Experiments 
conducted  by  Ding  et  al.  [Din86]  provided  direct  evidence  of  the  existance  of  a miscibility 
gap  by  high-resolution  electron  microscopy  and  electron  and  x-ray  diffraction  studies  on 
annealed  samples.  In  the  present  work  a direct  measurement  of  GaAs  activity  at  one 
composition  in  the  pseudobinary  section  as  a function  of  temperature  by  electrochemical 
technique  is  reported. 

The  galvanic  cell  used  in  the  experiment  is  shown  below  with  the  half  cell 
reactions: 

BzOafl),  GaaOsCs),  GaAs  (s).  As  (s)/ArSZ//GaxIni.xAs,  As  (s),  GajOsCs),  BjOjfl)  (IX) 

GaAs  (s,  pure)  -i-  3/2  O'^  = 1/2  GaaOs  + As  (s,  pure)  -t-  3e'  (2.23) 

1/2  Ga203  -I-  As  (s,  pure)  + 3e*  = GaAs  (ss)  + 3/2  0'^  (2.24) 

GaAs  (s,  pure)  = GaAs  (ss)  (2.25) 


where  ss  denotes  GaAs  in  the  solid  solution  GaxIn^xAs. 
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Table  2-10.  Summary  of  measured  activities  of  GaAs  in  the  pseudobinary  GaxIni.xAs  solid 

phase. 


Temperature 

EMF  (mV) 

^aAs 

AGoaAs 

y GaAs 

(K) 

(J/mol) 

(J/mol) 

873 

3.5 

0.961 

-1013.3 

20569 

2.13 

920 

5.6 

0.809 

-1621.2 

20569 

1.79 

928 

6.1 

0.802 

-1765.9 

20569 

1.78 

960 

8.75 

0.728 

-2533.1 

20569 

1.61 

968 

12.47 

0.638 

-3610.1 

20569 

1.41 

983 

13.30 

0.624 

-3850.4 

20569 

1.38 

1008 

15.00 

0.596 

-4342.5 

^ 

1.32 

1023 

13.5 

0.632 

-3474.0 

-25885 

1.47 

1056 

11.5 

0.684 

-3329.3 

-25885 

1.52 

1088 

9.5 

0.738 

-2750.3 

-25885 

1.64 

1113 

8.6 

0.764 

-2489.7 

-25885 

1.69 
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Experiments  were  conducted  to  verify  the  non-interference  of  B2O3  powder  with 
the  contents  of  the  electrode  material  in  both  the  compartments  as  discussed  earlier.  Once 
this  was  established  a ‘compression  design’  electrochemical  cell  was  assembled.  The 
working  electrode  was  designed  by  depositing  a 3 to  5p,  thick  film  of  Gao,45 1 Ino.549 As  on 
one  surface  of  the  solid  electrolyte  by  either  organometallic  chemical  vapor  deposition 
(MOCVD)  process  or  by  vapor  phase  epitaxy  (VPE)  done  at  Epitaxx,  Inc.  The  unwanted 
deposition  on  the  edges  of  the  YSZ  disk  were  etched  chemically  and  subsequently  cleaned 
by  rinsing  in  deionized  water.  Small  quantities  of  Ga203,  B2O3  and  As  powder  were 
sprinkled  uniformly  on  the  deposited  film.  The  preparation  of  the  reference  electrode  was 
fairly  simple.  The  constituent  powders  were  mechanically  mixed  and  placed  in  the 
reference  electrode  crucible.  The  quantity  of  B2O3  was  approximately  0.6  to  0.8g. 
Electrical  contact  was  made  by  means  of  tungsten  foils  sandwiched  between  the 
electrolyte  and  the  electrode  crucibles. 

The  open  circuit  emf  of  the  galvanic  cell  is  related  to  the  standard  chemical 
potential  difference  across  the  electrolyte  or  the  partial  Gibbs  energy  of  GaAs  in  the  solid 
solution.  At  equilibrium,  the  standard  Gibbs  energy  of  the  reaction  III  is  given  as: 

A°Gf  =-3FE  = RTln(ao3As,s)  (2.19) 

The  activities  obtained  by  these  electrochemical  measurements  are  presented  in  Figure  2- 
14  and  Table  2-10  as  a function  of  temperature. 

Summary  of  data  analysis.  The  activity  of  GaAs  in  the  pseudobinary  GaxIni.xAs 
solid  phase  is  a decreasing  function  of  temperature  in  the  experimental  range  of  873  to 
1113  K.  Beyond  this  temperature  range,  however,  the  activity  increases  monotonically 
with  temperature  as  shown  in  Figure  2-14.  It  appears  that  the  pseudobinary  section 


70 


GaAs-InAs  exhibits  a miscibility  gap  in  the  solid  phase.  The  temperature  at  which  the  emf 
vs  temperature  curve  shows  a discontinuity  in  temperature  dependence  (transition  from  a 
monotonically  increasing  to  monotonically  decreasing  function)  is  the  coordinate  of  the 
miscibility  gap  at  GaAs  = 0.451  composition  of  the  alloy  GaxIni.xAs.  The  activity  data  for 
temperatures  above  1008  K for  the  alloy  Gao.451Ino.549 As  in  the  liquid  single  phase  of  the 
pseudobinary  phase  diagram,  exhibit  positive  deviations  from  ideal  solution  behavior. 

The  present  results  are  in  disagreement  with  the  previous  emf  activity 
measurements  [Kre73]  that  reported  ideal  behavior  for  solutions  rich  in  both  GaAs  and 
InAs.  The  results  of  this  study,  however,  support  the  findings  of  Tmar  et  al.  [Tma88]  who 
identified  GaxIni-xAs  alloys  as  exhibiting  a miscibility  gap,  the  range  being  0.27  < XoaAs  < 
0.75  at  950  K.  They  estimated  the  critical  temperature  from  Knudson  cell  mass 
spectrometric  studies  of  InAs  activities  to  be  about  880  K.  This  is  significantly  higher  than 
the  value  of  700  K estimated  from  previous  calculations  [Tma87]  from  the  pseudobinary 
phase  diagram  data  and  estimates  of  the  interaction  parameter  assuming  regular  solution 
model  and  random  distribution  of  Ga  and  In  atoms  in  the  metal  sublattice.  The  present 
data  suggests  higher  critical  temperature,  an  estimate  being  between  950  to  1050  K.  This 
is  in  good  agreement  with  the  critical  temperature  estimate  by  Ding  et  al.  [Din86]  as  being 
in  the  range  of  848  to  923  K.  Shen  et  al  [She95]  performed  an  optimization  of  the  ternary 
system  in  the  pseudobinary  GaAs-In-As  phase  diagram  and  calculated  the  miscibility  gap 
in  the  solid  phase  region.  Their  calculation  showed  a critical  temperature  of  8 14.9  K. 

The  positive  deviations  from  ideal  solution  behavior  are  suggestive  of  a disordered 
state  and  strong  interaction  between  like  atoms  in  the  metal  sublattice  in  accordance  with 
the  findings  of  [Tma87]  suggesting  positive  excess  entropy  of  mixing.  The  positive 
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deviations  are  in  agreement  with  the  calculations  of  Ban  and  Ettenberg  [Ban73],  though 
the  measured  values  of  this  study  suggest  a higher  degree  of  non-ideality. 

In  summary,  the  present  work  illustrates  electrochemical  measurements  of  GaAs 
activity  in  the  pseudobinary  solid  region  of  the  GaAs-InAs  phase  diagram  confirming  the 
existence  of  a miscibility  gap  and  the  presence  of  strong  disorder  in  the  solid  phase. 
Earlier  optimizations  [Fos72,  Pan72,  Ant70,  CreSl,  Ona82a,  Ona82b,  Str69]  and 
estimates  from  physical  models  [Str82,  Str83,  Str72]  lead  to  values  for  the  interaction 
parameter  ranging  from  3 to  14.650  kJ/mol.  More  experiments  are  required  to  obtain  a 
complete  understanding  of  thermodynamic  behavior  in  the  pseudobinary  solid  region  and 
to  determine  precisely  the  range  of  immiscibility  in  the  solid  phase  as  well  as  the  critical 
temperature.  It  is  clear  from  the  present  work,  however,  that  the  interaction  parameter  has 
significant  temperature  dependence.  Similar  work  conducted  at  different  alloy 
compositions  will  yield  more  data  points  for  the  characterizing  the  solid  phase  in  the 
pseudobinary  cut  and  will  aid  in  providing  the  necessary  thermodynamic  information  for 


this  section  in  the  phase  diagram. 
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CHAPTER  3 

ASSESSMENT  STUDIES 
Assessment  ProcedurelLukas  Program! 

An  assessment  is  an  optimization  routine  in  which  optimal  values  for  the 
parameters  of  physical  or  chemical  models  describing  the  Gibbs  energies  of  phases  are 
estimated  from  the  available  thermodynamic  properties  of  a system  along  with  the  phase 
equilibria  data. 

The  values  of  the  estimated  parameters  permit  prediction  and  extrapolation  for 
unknown  properties  of  the  system.  Such  model  based  predictions  are  of  specific  interest 
to  systems  in  which  thermodynamic  property  determinations  are  limited  to  certain  regions 
of  the  phase  diagram  due  to  experimental  difficulties  (e.g.,  high  vapor  pressure  of  Group 
V elements)  or  to  systems  that  are  characterized  by  an  excessively  large  number  of 
experiments  that  contain  inconsistent  data.  Experiments  often  require  long  equilibration 
times  because  of  slow  diffusion  rates  in  the  solid  state.  It  is  thus,  economical  and  practical 
to  utilize  an  appropriate  thermodynamic  model  describing  the  thermochemical  behavior  of 
the  system.  The  model  parameters  can  be  used  to  make  predictions  in  unexplored  regions 
of  the  phase  domain,  and  in  the  case  of  ternary  and  higher  order  systems,  perform 
calculations  based  upon  the  determined  binary  limits. 
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The  first  step  in  an  assessment  is  to  choose  a suitable  model  for  all  phases, 
including  the  pure  component,  compound  and  the  solution  phases,  that  accurately 
describes  the  Gibbs  energy  of  each  phase.  For  solutions,  this  implies  the  formulation  of  a 

model  for  the  excess  Gibbs  energy  (Ag^^)  of  the  mixture  in  terms  of  the  temperature  and 
concentration  dependence.  It  is  desirable  to  use  as  few  terms  as  possible  in  the  expression 

for  Ag®*  in  order  to  avoid  extrapolation  difficulties.  It  is  frequent  that  very  different 
thermodynamic  formulations  result  in  very  similar  proposed  phase  diagrams  [Wag58].  If 
experimental  data  can  be  reasonably  well  represented  by  several  different  thermodynamic 
models,  a simpler  model  usually  provides  the  clearest  insight  into  the  basic  properties  of 
the  system  and  behaves  better  when  extrapolating  data.  Thus,  it  is  important  to  check  the 
adopted  thermodynamic  model  for  physically  unrealistic  parameters. 

V arious  types  of  solution  models  have  been  used  for  semiconductor  solutions.  In 
particular,  the  Redlich-Kister  solution  model  has  been  found  to  adequately  fit 
thermodynamic  data  for  most  of  III-V  semiconductor  systems.  Owing  to  it’s  widely 
accepted  use  for  representation  of  thermodynamic  properties  of  III-V  semiconductor 
systems  and  ease  of  making  reliable  extrapolations,  this  model  was  chosen  to  represent  the 
reported  enthalpy,  component  activities  and  phase  diagram  data  in  this  study.  The 
expression  is  as  follows: 


where  the  coefficients  Ly 
[Luk77]: 


V) 


(3.1) 


are  functions  of  temperature  according  to  Lukas  description 


- Bij^">T  + C>^(l-lnT)  - Dy^">  T'/2  - Eij<"V2T  - Fij<"'rV6 


(3.2) 
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where  the  coefficients  By^"*  ....  Fy^"^  are  the  model  parameters  fitted  to  the 

experimental  data  for  the  liquid  phase  in  the  optimization  routine  of  the  Lukas  software. 
The  coefficients  Ay  and  By^  * correspond  to  the  temperature  independent  contributions  to 
the  excess  enthalpy  and  entropy  of  mixing  in  the  liquid  phase.  Whereas,  the  terms  Cy*"^ ... 

Fy  are  the  specific  heat  terms  in  the  formation  equation  for  the  liquid  phase  from  the 
elements  in  their  standard  states. 

The  even  values  of  index  n yield  symmetric  ftinctions  in  composition  and  tend  to 

flatten  the  parabolic  nature  of  AG" . The  odd  values  of  n,  on  the  other  hand,  tend  to 
skew  the  parabolic  function.  Similar  coefficients  ay‘"* , by‘"’ ....  fy‘"'  are  obtained  during  the 
optimization  of  model  parameters  for  the  solid  compound  phase,  in  which  case  the 
expression  for  the  coefficient  for  the  solid  phase  is  as  follows: 

= aij‘">  - by<"'T  + Cy<">T(  1 -InT)  - dy<"^  TV2  - ey‘"V2T  - flj‘">TV6  (3.3) 

The  parameters  in  this  expression  have  the  same  connotation  in  the  solid  phase  as 
do  the  parameters  described  above  for  the  liquid  phase.  The  solid  III-V  compound  (e.  g., 
GaAs)  is  regarded  as  a line  compound,  because  of  the  typically  narrow  region  of  solid 
solution  about  its  ideal  stoichiometry  of  x^=  0.5.  The  pure  liquid  phase  was  chosen  as  the 
reference  state  for  the  liquid  solution.  For  the  stoichiometric  compound,  the  solid 
elements  in  their  standard  states  were  chosen  as  the  references  for  the  compound;  the 
choice  between  pure  solid  and  liquid  reference  states  is  entirely  arbitrary. 

Negligible  mutual  solubility  between  the  pure  solid  Group  III  and  V elements  was 
assumed.  Pure  component  SGTE  (Scientific  Group  Thermodata  Europe)  lattice  stabilities 
were  employed  to  express  the  Gibbs  energies  of  the  elements  in  the  stable  phase  relative  to 


76 


the  enthalpy  of  the  "Standard  Element  Reference"  { i.e,  the  enthalpies  of  the  pure  elements 
in  their  defined  phases  at  10^  Pa  and  298.15  K (denoted  as  G-Hser)  [Din89]}. 

The  Gauss  method  is  adopted  in  the  Lukas  software  for  least  square  fitting  the 
data.  The  strategy  of  the  calculations  was  based  upon  the  Newton-Raphson  method. 
Marquardt  s algorithm,  which  is  a combination  of  Newton-Raphson  method  and  the 
steepest  descent  method,  was  best  suited  to  calculate  complete  phase  diagrams  and  their 
sections.  Flexibility  is  assigned  in  the  optimization  by  the  ability  to  choose  different 

weighting  factors  for  the  source  literature,  and  is  taken  as  the  reciprocal  of  the  estimated 
accuracy  of  the  measured  values. 

Two  different  versions  of  error  equations  were  utilized,  one  being  an  explicit 
expression  while  the  other'  an  implicitely  defined  function.  The  choice  of  the  error 
equation  is  governed  by  the  type  of  data  in  hand,  the  expression  chosen  for  the  adjustable 
coefficients  and  partly  by  the  stage  of  optimization.  Especially  in  the  early  stages  of 
optimizations,  where  the  coefficients  are  far  from  the  best  values,  the  explicit  version  more 
probably  directs  the  calculation  to  convergence  than  the  implicit  expression  of  the  error 
equation. 

The  mean  square  of  error  was  the  criterion  used  to  obtain  the  best  fit,  however,  the 
lowest  value  did  not  necessarily  imply  the  best  thermodynamic  description  of  the  system. 
Examining  the  fit  of  each  experimental  data  point  to  the  values  calculated  from  the 
optimization  together  with  the  mean  square  of  error  for  the  entire  dataset  was  used  as  the 
criterion  for  a good  fit. 
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Assessments  of  Binary  Systems 

Ga-As  Assessment 

Ga-As  system  has  been  investigated  for  it’s  thermodynamic  properties  and  phase 
diagram  by  many  investigators  as  discussed  in  section  2.2.1.  The  selected  data  set  for 
phase  diagram,  enthalpy  and  Gibbs  energy  for  Ga-As  system  are  shown  in  Tables  3-1 
through  3-3  together  with  the  experimental  method  used,  the  range  of  measurements,  and 
the  estimated  errors  in  the  techniques  employed. 

Redlich-Kister  polynomial  was  used  to  model  the  liquid  phase,  while  a line 
compound  description  was  used  for  the  essentially  stoichiometric  compound  GaosAsos. 
The  pure  liquid  was  chosen  as  the  reference  phase  for  the  Ga-As  system.  However,  for 
the  stoichiometric  compound,  the  solid  elements  in  their  standard  states  were  chosen  as 
the  reference  phases.  Negligible  solubility  of  Ga  in  As  and  of  As  in  Ga  was  assumed,  this 
assumption  being  consistent  with  other  studies  [Hal63  and  PerSO]. 

Four  terms,  two  for  the  liquid  phase  and  two  for  soild  GaAs,  Ay^"* , ay*"’  and 

by  , were  sufficient  to  model  the  sytem.  The  terms  Cy^"^  and  Cy^"^  did  not  significantly 
improve  the  fit,  and  therefore  were  not  included  in  the  final  optimization.  This  is 
consistent  with  the  low  values  available  for  the  experimental  ACp  values  for  GaAs  (1)  and 
GaAs  (s).  The  selected  specific  heat  data  and  their  references  are  listed  in  Table  3-2.  The 
phase  diagram  is  fairly  well  characterized  for  the  system,  however,  the  thermodynamic 

information  in  the  Xa$  > 0.5  region  and  the  high  temperature  region  of  the  phase  diagram 
is  lacking. 
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Standard  steps  were  followed  to  perform  the  optimization  as  outlined  in  the 
previous  section.  The  salient  features  are: 

1.  The  coefficients  a3i  and  bai  were  set  to  values  determined  by  the  temperature 
dependent  expression  for  A“Gj  obtained  from  the  present  experiments. 

2.  Next,  using  only  the  Gibbs  energy  data  of  Lee  et  al.  [Lee92]  and  this  work  the  four 
coefficients  Ay , By,  ay  and  by  were  optimized  to  obtain  better  start  estimates  of  the 
parameter  and  also  to  ensure  physical  realism  to  the  parameters. 

3.  Using  these  coefficients  and  only  the  Gibbs  energy  data,  a reasonably  good  prediction 

was  made  for  the  phase  diagram  and  the  liquid  phase  enthalpy.  However,  the  asymmetry 

reported  in  the  experimental  phase  diagram  data  [Sol74,  Kos55,  Hsi77,  Hal63,  Rub66, 

Rak70,  Vig73,  PerSO,  DeW85]  and  the  enthalpy  of  mixing  data  of  Yamaguchi  [Yam88] 
remained  unfitted. 

4.  The  asymmetry  was  introduced  into  the  model  by  including  another  term  A51  in  the 

optimization.  This  parameter  skews  the  parabolic  nature  of  the  liquid  phase  excess  Gibbs 
energy  function. 

5.  Using  various  combinations  of  only  the  enthalpy  and  Gibbs  energy  data,  optimized 
coefficients  were  obtained  and  mean  square  error  compared.  Consistency  was  checked  by 
the  ability  of  the  coefficients  to  predict  the  experimental  phase.  Inconsistent  enthalpy  and 
Gibbs  energy  data  were  then  eliminated  if  an  acceptable  agreement  was  not  achieved. 

6.  The  coefficients  were  re-optimized  using  the  combined  thermochemical  and  phase 
diagram  data.  The  coefficients  yielding  an  optimal  fit  were  retained. 

7.  Lastly,  the  optimization  was  repeated  by  including  the  coefficients  C43  and  C33  to  check 
for  significant  Cp  contribution.  Since  the  fit  to  the  data  upon  their  inclusion  was  not 
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significant,  consistent  with  the  experimental  Cp  values,  these  parameters  were  excluded 
from  the  model  description. 

As  many  as  3 1 runs  were  performed  for  the  assessment  of  the  Ga-As  system.  The 
optimized  coefficients  for  the  Redlich-Kister  model  and  lattice  stability  of  GaAs  are  shown 
3-4.  The  comparision  of  experimental  data  with  the  optimized  curves  is  shown  in 
Figure  3-1  through  Figure  3-3.  The  general  agreement  of  the  curves  appears  to  be  good. 

When  compared  to  Chatillons  assessment  [Cha90],  the  calculated  liquidus  appears 
to  have  shifted  to  the  As-rich  side.  This  trend  was  also  seen  in  Lee's  experimental 
measurements  [Lee92]  on  Ga  activity  along  the  liquidus.  Chatillon's  assessment  [Cha90] 
predicts  higher  Ga  activities  especially  at  lower  temperatures,  while  results  of  this  work 
and  Lee's  measurements  [Lee92]  predict  lower  values  for  the  gallium  activity  and  stronger 
deviations  from  ideal  solution  behavior.  In  the  low  temperature  regions,  the  phase 
diagram  data  by  Sol  et  al.  [Sol74],  Hsieh  [Hsi77]  and  Hall  [Hal63]  were  inconsistent, 

i 

mostly  due  to  the  experimental  difficulties  encountered  in  accurate  observation  of 
dissolution  and  precipitation  of  crystals  in  gallium  rich  melts  during  the  heating  and 
cooling  cycles.  Subcooling  could  also  have  been  an  issue  during  the  cooling  cycle 
resulting  in  untrue  measurements.  The  As-rich  eutectic  data  of  Ceolin  et  al.  [Ceo89]  was 
determined  more  consistent  than  that  reported  by  Koster  and  Thoma  [Kos55]. 

The  standard  enthalpy  of  formation  of  GaAs  (s)  at  298.15  K agreed  well  with  that 
determined  calorimetrically  by  Martosudirjo  and  Pratt  [Mar74].  The  value  suggested  by 
Gutbier  [Gut61]  from  mass  spectrometry  is  too  low  (numerically).  The  enthalpy  of  fusion 
data  by  Lichter  [Lic69]  and  Yamaguchi  [Yam88]  were  retained  in  the  selected  data  set. 
The  liquid  phase  enthalpy  of  mixing  data  by  Yamaguchi  [Yam88]  was  well 
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Figure  3-2.  Optimized  Ga-As  liquid  phase  enthalpy. 
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represented.  However,  the  Cp  (T)  data  of  Lichter  [Lic69]  was  found  inconsistent  with 
other  data.  It  was  concluded  that  there  exists  marked  discrepancy  in  the  reported  Cp  (T) 
data,  and  this  was  evidenced  by  a third  law  analysis  discussed  below.  For  this  reason,  the 
expression  for  optimized  Cp  (T)  data  from  Chatillon's  least  square  fit  [Cha90]  was  adopted 
for  third  law  calculations  in  this  work. 

The  emf  studies  by  Krestovnikov  et  al.  [Kre71],  Sirota  [Sir68]  and  Abbasov  et  al. 
[Abb66]  for  the  measurement  of  Gibbs  energy  of  formation  of  GaAs  were  found 
inconsistent  with  the  phase  diagram  data  and  rest  of  the  thermochemical  properties.  The 
measurements  for  activity  of  gallium  in  the  liquid  phase  by  Wypartowicz  and  Fitzner 
[Wya88],  Hall  [Hal63]  and  Katayama  et  al.  [Kat89]  were  rejected  for  the  same  reason. 

To  summarize,  the  calculated  values  of  Gibbs  energy  values  for  Ga  activity  along 
the  liquidus  compare  well  with  the  experimental  data  of  Lee  et  al.  [Lee92]  as  is  shown  in 
Table  3-5  and  Figure  3-3.  The  partial  pressures  associated  with  As2  and  As4  vapor  phases 
in  the  Ga  -rich  segment  of  the  phase  diagram  are  shown  in  Figure  3-4.  The  Gibbs  energies 
of  formation  obtained  from  the  optimization  is  also  in  good  agreement  with  the 
experimental  data  of  the  present  work  as  can  be  seen  in  Table  3-6. 

In- As  Assessment 

Redlich-Kister  polynomial  was  well  suited  to  model  the  liquid  phase,  while  a line 
compound  description  was  used  for  the  essentially  stoichiometric  compound  Ino,5Aso.5. 
Either  solid  or  liquid  phase  could  be  used  as  the  reference  for  the  system.  As  with  the 
case  of  Ga-As  system,  pure  liquid  was  chosen  as  the  reference  phase  while  solid  elements 
in  their  standard  states  were  chosen  as  the  reference  for  the  stoichiometric  compound. 
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Figure  3-4.  Vapor  pressure  of  As  in  Ga-rich  region  of  Ga-As  phase  diagram. 
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Table  3-1 . Summary  of  selected  experimental  phase  diagram  data  for  Ga-As  system. 


Author  [Reference] 

Method  of 

Measurement 

Estimated 

Measurement 

Range 

Uncertainity 

T(K) 

+/-  AT  (K) 

Xas 

+/-  Axas 

Koster  [Kos55] 

Differential  Thermal 

1273-1511 

20 

Analysis 

0.15-0.80 

0.01 

Rubenstein  [Rub66] 

Filtration  Technique 

723-1511 

20 

— 

0.05 

Rakov  et  al. 

Dew  Point  and 

1447-1505 

5 

[Rak70] 

Continuous 

Weighing 

0.31-0.54 

0.05 

Vigdorowich 

Bourdon  Gauge  and 

1253-1353 

20 

[Vig73] 

Vapor  Density 

0.10-0.22 

0.05 

Perea  [PerSO] 

Heterogeneous 

853-943 

5 

Equilibrium  and 
Weight  Loss 

— 

0.05 

DeWinter  [DeW85] 

Visual  Observation 

, 623-923 

20 

and  Source 

0.05 

Dissolution 

Y amaguchi 

Drop  Calorimetry 

800-1500 

5 

[Yam88] 

— 

0.03 

Ceolin  et  al. 

Neutron  Diffraction 

1083-1443 

15 

[Ceo89] 

— 

0.02 
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Table  3-2.  Summary  of  selected  calorimetric  and  enthalpy  data  for  Ga-As  system. 


Author 

[Reference] 

Method  of 
Measurement 

Data  Type 

Temperature 
Range  or  Value 

Estimated 

Uncertainity 

Sirota  [Sir68] 

Oxygen 

Calorimetric 

Bomb 

Enthalpy  of 
Formation 

-43848  (J/mol) 
298  K 

■¥/-  5000  J/mol 

Lichter  [Lic69] 

Drop 

Calorimetry 

Enthalpy  of 
Fusion 

52702  (J/mol) 
1513K 

+/-  2500  J/mol 

Lichter  [Lic69] 

Drop 

Calorimetry 

Heat  Capacity 

421-1513  K 

+!-  10  K 

Martosudirjo 

[Mar74] 

Calorimetric 

Precipitation 

« 

Enthalpy  of 
Formation 

-41148  J/mol 
298  K 

+!-  5000  J/mol 

Yamaguchi 

[Yam88] 

Drop 

Calorimetry 

Enthalpy  of 

Formation 

1 

800-1550  K 

+1-  10  K 
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Table  3-3.  Summary  of  selected  vapor  pressure  and  Gibbs  energy  data  for  Ga-As  system. 


Author  [Reference] 

Method  of 

Temperature  Range 

Estimated 

Measurement 

(K) 

Uncertamity 

+!-  T (K) 

/^Asj 

Drowart  [Dro58] 

Knudsen-cell  Mass 

1030-1136 

10 

Spectrometry 

30 

Richman  [Ric63] 

Standard  Bourdon 

1323-1527 

3 

Gauge 

20 

Pupp  [Pup64] 

Knudsen-cell  Mass 

1035-1239 

10 

Spectrometry 

20 

Arthur  [Art67] 

Knudsen-cell  Mass 

890-1190 

10 

Spectrometry 

20 

Rakov  [Rak70] 

Continuous 

1448-1511 

3 

Weighing 

20 

Foxon  [Fox73] 

Knudsen-cell  Mass 

850-1100 

10 

Spectrometry 

20 

Khukhryanskii 

Dynamic  Method 

973-1273 

10 

[Khu74] 

30 

Panish  [Pan74] 

H2  + PH3  Flow 

1255-1435 

10 

Method 

20 

Lee  [Lee92] 

EMF 

923-1373 

5 

This  Work 

EMF 

813-1013 

5 

5 
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Table  3-4.  Optimized  coefficients  for  Redlich-Kister  model  and  lattice  stability  of  GaAs. 


Ajj  (J/mol) 

Bij  (J/mol  K) 

Suffix  (i,j) 

j=I 

j=2 

i=4 

-26966.58 

4.06734 

i=5 

12947.10 

10.11249 
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Table  3-5.  Calculation  of  partial  Gibbs  energies  of  Ga  and  As  along  liquidus  (in  J/mol). 


T 

(K) 

Xas 

Hca 

< 

II 

llGa'' 

InPls, 

inp;. 

YCa 

YAs 

RTInyca/xAs^ 

923 

0.005 

-39 

-67652 

-32 

-22.9 

-34.28 

1.000 

0.030 

— 

948 

0.007 

-57 

-66203 

— 

-21.5 

-32.24 

1.000 

0.031 

— 

973 

0.010 

-83 

-64746 

-67 

-20.13 

-30.22 

1.000 

0.033 

— 

998 

0.013 

-117 

-63281 

— - 

-18.82 

-28.28 

0.999 

0.036 

-49121.4 

1023 

0.018 

-163 

-61804 

-198 

-17.58 

-26.45 

0.999 

0.038 

-26263.8 

1048 

0.024 

-224 

-60312 

-265 

-16.40 

-24.71 

0.998 

0.041 

-30284.0 

1073 

0.031 

-302 

-58803 

-349 

-15.26 

-23.02 

0.997 

0.044 

-27890.7 

1098 

0.039 

-401 

-57274 

-453 

-14.17 

-21.40 

0.996 

0.048 

-24055.4 

1123 

0.048 

-524 

-55719 

-512 

-13.12 

-19.84 

0.993 

0.053 

-28466.2 

1148 

0.059 

-677 

-54135 

-723 

-12.12 

-18.35 

0.989 

0.058 

-30327.8 

1173 

0.071 

-864 

-52518 

-988 

-11.15 

-16.90 

0.985 

0.064 

-29238.9 

1198 

0.084 

-1089 

-50861 

-1228 

-10.21 

-15.49 

0.980 

0.072 

-28517.9 

1223 

0.099 

-1359 

-49160 

-1511 

, -9.30 

-14.13 

0.971 

0.080 

-30530.8 

1248 

0.115 

-1680 

-47408 

-1901 

-8.42 

-12.80 

0.961 

0.090 

-31210.7 

1273 

0.132 

-2060 

-45598 

-2286 

-7.56 

-11.50 

0.949 

0.102 

-31796.4 

1298 

0.151 

-2509 

-43718 

-2746 

-6.73 

-10.24 

0.933 

0.115 

-32822.0 

1323 

0.171 

-3038 

-41758 

-3189 

-5.90 

-8.97 

0.915 

0.131 

-33415.1 

1348 

0.193 

-3663 

-39702 

-4009 

-5.10 

-7.74 

0.891 

0.150 

-34796.3 

1373 

0.216 

-4406 

-37528 

-4996 

-4.30 

-6.50 

0.870 

0.173 

-34072.6 

^ Denotes  Ga  activity  from  Lee  [Lee92] . 

* Denotes  calculated  values  using  optimized  (0. 1 MPa)  and  Kp  (T)  data  from  Chatillon 
et  al.  [Cha90]. 
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Table  3-6.  Comparision  of  calculated  Gibbs  energy  of  formation  of  GaAs  (s)  referred  to 

Ga  (1)  and  As  (s)  with  experimental  data. 


Temperature  (K) 

A G®* , Optimized  (J/mol) 

A G®^ , Present  Work  (J/mol) 

813 

-33889.3 

— 

833 

-33541.2 

-33650 

853 

-33193.9 

-32420 

873 

-32844.9 

-31850 

893 

-32496.8 

-30970 

913 

-32148.7 

-30040 

933 

-31800.5 

-29460 

953 

-31452.4 

-29100 

973 

' ‘-31104.3 

-28740 

993 

-30756.2 

( 

— 

1013 

-30408.0 

— - 
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Negligible  solubility  of  In  in  As  and  of  As  in  In  was  assumed,  this  assumption  being 
consistent  with  phase  diagram  studies  [Hal63  and  PerSO]. 

Four  terms  Ay  , By,  ay  and  by  were  used  to  model  the  sytem.  The  terms  Cy  and  Cy 
did  not  significantly  improve  the  fit,  and  therefore  were  not  included  in  the  final 
optimization.  This  is  consistent  with  the  low  values  available  for  the  experimental  ACp 
values  for  InAs  (1)  and  InAs  (s).  The  selected  experimental  data  along  with  their 
methods  of  determination  and  estimated  uncertainities  for  phase  diagram,  enthalpy  and 
Gibbs  energy  for  In-As  system  are  listed  in  Table  3-7  through  3-9.  The  phase  diagram  is 
fairly  well  characterized  for  the  system,  however,  the  thermodynamic  information  in  the 
^As  > 0.5  region  and  the  high  temperature  region  of  the  phase  diagram  is  lacking.  Owing 
to  the  availability  of  enthalpy  of  mixing  values  only  for  the  Xas  < 0.5  region,  only  one 
enthalpy  term  A41  was  used  for  optimization. 

Standard  steps  of  optimization  similar  to  the  ones  utilized  on  the  Ga-As  system 
were  used  for  this  system  as  well.  First,  the  coefficient  A41  was  estimated  by  forcing  a fit 
to  the  liquid  phase  enthalpy  of  mixing  data  of  Yamaguchi  et  al.  [Yam88].  Next,  the 
coefficients  ay  and  by  were  estimated  based  upon  the  Gibbs  energy  of  formation  data  of 
InAs  as  determined  in  the  present  study.  The  coefficients  were  then  set  to  vary  in  a 
scheme  of  subsequent  optimizations  to  least  square  fit  the  entire  experimental  dataset. 

Each  data  point  was  checked  for  consistency  with  other  data  points  by  comparing 
the  degree  of  fit  to  the  model,  as  well  as  by  the  ability  to  back  calculate  other  quantities. 
Optimizations  were  performed  to  calculate  thermod)oiamic  data  from  pure  phase  diagram 
data.  In  a similar  set  of  reverse  calculations,  the  consistency  of  the  Gibbs  energy,  enthalpy 
and  pure  component  data  was  checked  by  calculating  the  system  phase  diagram  and 
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comparing  with  experimental  data.  This  was  helpful  particularly  because  the  phase 
diagram  is  the  most  well  characterized  information  available  on  the  system  as  opposed  to 
the  thermochemical  properties.  When  such  asssessments  were  completed  and  consistent 
data  separated  from  inconsistent  ones,  model  parameters  were  determined  and  are  shown 
in  Table  3-10. 

The  optimization  compared  fairly  well  with  the  recent  assessment  by  ChatUlon  et 
al.  [Cha90]  however,  at  temperatures  above  893  K there  is  a shift  of  the  phase  diagram 
towards  the  In-rich  region  as  shown  in  Figure  3-5.  This  indicated  that  the  calculated  In 
activities  were  higher  than  those  determined  by  pure  phase  diagram  data  or  vapor  pressure 
measurements  of  As.  Vapor  pressure  measurements  of  As  showed  great  discrepancy 
across  the  phase  diagram,  and  thus  were  not  reliable  tools  for  performing  accurate 
calulations.  Gibbs  energy  data,  of  formation  Of.  InAs  and  activity  of  In  was  used  to 
perform  predictions  of  the  phase  data.  The  Gibbs  energy  data  skewed  the  phase  diagram 
data  towards  In-rich  region  and  flattened  the  curve  at  mid  regions  of  the  phase  diagram. 
Lack  of  accurate  thermodynamic  data  in  the  As-rich  regions  of  the  phase  diagrams,  made 
it  difficult  to  provide  a better  fit  to  the  As-rich  segment  of  the  system. 

The  data  by  Itoh  et  al.  [Ito82]  compared  well  with  the  activity  data  of  the  present 
study  however,  exhibiting  a lower  degree  of  non-ideality.  Calorimetric  measurements  for 
enthalpies  of  fusion  data  by  Cox  and  Pool  [Cox67]  and  Lichter  and  Sommelet  [Lic69] 
were  used  in  the  optimization.  The  liquid  phase  enthalpy  data  reported  by  Yamaguchi  et 
al.  [Yam88]  was  best  represented  by  the  calculated  curve  as  shown  in  Figure  3-6. 
Invariant  equilibria  by  Liu  and  Peretti  [Liu53]  compared  well  with  the  calculated  curve. 
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The  optimized  curves  agree  well  with  experimental  data  as  is  evident  through  Figures  3-5 
through  3-7. 

Vapor  pressure  measurements  for  As2  and/or  As4  species  have  been  performed  by 
Knudsen  cell  mass  spectrometry  [Pup74,  Gut59  and  Gol59].  Total  arsenic  pressure 
measurements  have  been  performed  by  dew  point  techniques  [Van57,  Fol58,  Kar73], 
Bourdon  gauges  [Vig73,  Ric63]  and  gas  flow  methods  [Khu74,  Khu77].  In  order  to 
compare  all  the  published  results,  it  is  necessary  to  transform  the  data  by  recalculating  the 
partial  pressures  of  As2  and  As4  and  the  total  pressures.  Tmar  et  al.  [Tma84]  performed 
such  calculations  and  compared  the  corrected  experimental  results  with  each  other.  They 
found  significant  discrepancy  in  the  reported  data.  The  results  of  Khukhryanskii  et  al. 
[Khu74]  did  not  agree  with  our  activity  data  and  the  calculated  curve.  This  is  consistent 
with  the  analysis  of  their  data  performed  by  Tmar  et  al.  [Tma84]  indicating  that  these 
measurements  were  erroneous  because  they  treated  As2  as  the  main  species  while  ignoring 
the  arsenic  component  in  AS4,  AsHa,  ASH2  and  AsH  species.  The  situation  was 
aggravated  because  these  variations  were  temperature  dependent  and  uncertainity  in  the 
calculated  enthalpies  caused  the  activity  of  As  to  change  seriously.  Parasitic 
reevaporation,  dissociative  ionization  of  AS4  into  As2‘^,  the  reversibility  of  evaporation  and 
the  calibration  of  the  mass  spectrometer  are  critical  to  the  authenticity  of  the  vapor 
pressure  data  reported.  Tmar  et  al.  [Tma84]  corrected  the  values  of  Goldfinger  and 
Jeunehomme  [Gol59]  by  calculating  the  equilibrium  dissociation  constant  of  the  reaction 
(2.1)  and  the  calibration  constant  of  their  apparatus.  This  brought  their  values  closer  to 
the  values  by  Pupp  et  al.  [Pup74],  though  they  were  still  lower  and  were  thus  discarded. 
High  vapor  pressure  measurements  of  Van  den  Boomgaard  and  Schol 
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Figure  3-5.  Optimized  In-As  phase  diagram. 
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[Van57]  and  Karataev  et  al.  [Kar73]  especially  on  the  In-rich  side  of  the  In-As  system  do 
not  agree  with  the  results  of  the  present  work  and  the  data  reported  by  Pupp  et  al. 
[Pup74].  This  is  probably  due  to  the  fact  that  InAs  forms  a layer  at  the  surface  of  the 
indium  melt  creating  a heat  transfer  barrier,  thereby  causing  difficulties  in  accurate  dew 
point  temperature  measurement.  However,  their  data  agreed  well  at  high  arsenic  contents. 

The  emf  data  by  Abbasov  et  al.  [Abb64]  and  Krestovnikov  et  al.  [Kre71]  were 
retained  during  the  earlier  stages  of  optimization,  but  later  were  discarded  because  were 
found  inconsistent  with  the  low  temperature  data.  The  data  by  Krestovnikov  et  al. 
[Kre71]  however,  showed  better  agreement  with  the  calculated  curve  and  the 
measurements  of  the  present  study  when  extrapolated  to  the  high  temperature  regions  of 
the  phase  diagram  of  the  In-As  system.  The  discrepancy  is  probably  due  to  poor  stability 
of  cells  in  their  experiments  due  to  poor  diffusion  or  indium  surface  migration  on  the  two 
phase  mixtures  (InAs  -i-As). 

The  calculated  values  of  Gibbs  energy  values  for  In  activity  along  the  liquidus 
compare  well  with  the  experimental  data  of  this  work  and  are  summarized  in  Table  3- 1 1 
and  Figure  3-7.  The  partial  pressures  associated  with  As2  and  As4  vapor  phases  in  the  In- 
rich segment  of  the  phase  diagram  are  shown  in  Figure  3-8.  The  Gibbs  energies  of 
formation  obtained  from  the  optimization  is  also  in  good  agreement  with  the  experimental 
data  of  the  present  work  as  can  be  seen  in  Table  3-12. 
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Figure  3-8.  Vapor  pressure  of  As  in  In-rich  region  of  the  In- As  phase  diagram. 
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Table  3-7.  Summary  of  selected  experimental  phase  diagram  data  for  In- As  system. 


Author  [Reference] 

Method  of 
Measurement 

Measurement 

Range 

T(K) 

Xas 

Estimated 

Uncertainity 

+/-  AT  (K) 

+/-  AXas 

Liu  [Liu53] 

Differential  Thermal 

429-1215 

10 

Analysis  and  X-ray 

Analysis 

0-1 

0.03 

Hall  [Hal63] 

Heterogeneous 

673-1024 

10 

Equilibrium  and 

Weight  Loss 

0-0.16 

0.02 

Karataev  [Kar73] 

Dew  Point  and 

1149-1220 

10 

Continuous 

Weighing 

0.35-0.74 

0.02 

Perea  [PerSO] 

Heterogeneous 

849-929 

10 

e 

? 

Equilibrium  and 

Weight  Loss 

0-0.08 

0.02 

Knobloch  [Kno85] 

Dissolution  in  LPE 

724-906 

5 

0.007-0.07 

0.01 

DeWinter  [DeW86] 

Visual  Observation 

648-914 

5 

and  Source 

Dissolution 

0.002-0.08 

0.02 
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Table  3-8.  Summary  of  selected  calorimetric  and  enthalpy  data  for  In- As  system. 


Author 

Method  of 

Data  Type 

Temperature 

Estimated 

[Reference] 

Measurement 

Range  or  Value 

Uncertainity 

Schottky 

Dissolution 

Enthalpy  of 

-30962  J/mol 

+/-  3000  J/mol 

[Sch58] 

Calorimetry 

Formation 

298  K 

Sharifov 

Calorimetric 

Enthalpy  of 

-28870  J/mol 

+/-  3000  J/mol 

[Sha63] 

Bomb 

Formation 

298  K 

Cox  [Cox67] 

Drop 

Enthalpy  of 

38511  J/mol 

+/-  3000  J/mol 

Calorimetry 

Fusion 

1213  K 

+/-5K 

Cox  [Cox67] 

Drop 

Calorimetry 

Heat  Capacity 

298-1200  K 

+/-  10  K 

Lichter  [Lic69] 

Drop 

Heat  Capacity 

298-1200  K 

+/-  10  K 

' 

Calorimetry 

Glazov  [Gla75] ' 

Differential 

Enthalpy  of 

39330  J/mol 

+/-  3000  J/mol 

Thermal 

Analysis 

Fusion 

1215  K 

+/-5K 

Y amaguchi 

Drop 

Enthalpy  of 

38300  J/mol 

+/-  3000  J/mol 

[Yam88] 

Calorimetry 

Fusion 

1221  K 

+/-  10  K 

Yamaguchi 

Drop 

Enthalpy  of 

700-1500  K 

+/-  10  K 

[Yam88] 

Calorimetry 

Formation 
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Table  3-9.  Summary  of  selected  vapor  pressure  and  Gibbs  energy  data  for  In- As  system. 


Author  [Reference] 

Method  of 
Measurement 

Temperature  Range 
(K) 

Estimated 
Uncertainity 
+/-  T (K) 

AP  As2/P  As2 

Van  den  [Van57] 

Dew  Point  Method 

1004-1215  K 

10 

20 

Pupp  [Pup74] 

Knudsen-cell  Mass 

918-1084 

10 

Spectrometry 

50 

Itoh  [Ito82] 

EMF  Measurement 

920-1260 

5 

5 

Tmar  [Tma87] 

Knudsen-cell  Mass 

846-1086 

5 

Spectrometry 

10 

This  Work 

EMF  Measurement 

833-1153 

5 

1 

5 
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Table  3-10.  Optimized  coefficients  for  Redlich-Kister  model  and  lattice  stability  of  InAs. 


Aij  (J/mol) 

Bij  (J/mol  K) 

Suffix  (i,j) 

j=l 

j=2 

i=4 

-21600 

0.0 

i=5 

0.0 

0.0 

°G^as,s-0.5‘*G.„..-0.5°G 


In,bct 


As.rho 


-25113.3-I-7.780T 


Table  3-11.  Calculation  of  partial  Gibbs  energy  of  In  along  liquidus. 


cs 

< 

X 

c 

H 

-38191.5 

-35416.9 

-33472.1 

-31827.5 

-29171.3 

-26775.4 

-24098.4 

1009 13- 

-19938.0 

o 

oo 

1 

< 

8110 

0.146 

0.170 

96  TO 

0.236 

0.282 

0.341 

o 

• 

o 

O 

• 

o 

0.511 

Yin 

OOOT 

9660 

0.952 

00 

00 

00 

• 

O 

006  0 

O 

00 

• 

o 

0.756 

ouo 

0.711 

0390 

» ' ft- 

c ^ 
o 

-24.32 

-19.20 

-16.72 

-14.58 

-11.93 

-9.70 

-7.46 

1 

-5.54 

-4.50 

-3.13 

* Oh 

k*'  S 

c ^ 

o 

-19.47 

-16.01 

-14.20 

-12.60 

-10.70 

-9.08 

-7.48 

-6.13 

-5.40 

-4.46 

t "o 

< E 

=i-  ^ 

-38191.5 

-35416.9 

-33472.1 

-31827.5 

-29171.3 

-26775.4 

-24098.4 

-21600.1 

-19938.0 

-18073.3 

* 

£ 

03 

0.962 

0.936 

0.912 

00 

00 

• 

o 

0.838 

0.786 

0.720 

0.652 

0.605 

0.550 

* r-  "O 

s a 

H < 

-265 

-482.2 

-704.2 

-964. 1 

-1463.4 

-2073.5 

-2937.5 

-3943.6 

-4714.4 

-5728.5 

cd 

0.971 

0.941 

(N 

00 

00 

• 

o 

908  0 

00 

00 

• 

o 

0.733 

0.607 

0.542 

0.523 

9ZV0 

Uln 

(J/mol) 

-202.7 

-442.9 

-955.4 

-1708.1 

-1968.6 

-2663.4 

-4458.3 

-5645.3 

-6079.5 

-7961.3 

< 

X 

o 

• 

o 

0.0549 

0.0733 

0.0921 

OPZVO 

O 

• 

o 

0.1970 

0.2370 

0.2645 

0.2970 

T 

(K) 

833 

i 

00 

00 

918 

953 

994 

1033 

1073 

1108 

1128 

1153 

**  Denotes  calculated  values  (this  study). 


104 


Table  3-12.  Comparision  of  calculated  Gibbs  energy  of  formation  of  InAs  (s)  referred  to 

In  (1)  and  As  (s)  with  experimental  data. 


Temperature  (K) 

Ag®^,  Optimized  (J/mol) 

A G^^ , Present  Work  (J/mol) 

803 

-34880 

-34880 

833 

-34184 

-34300 

853 

-33720 

-33880 

867 

-33395 

-33300 

873 

-33256 

-33100 

891 

-32838 

-33000 

913 

-32328 

-32680 

921 

-32142 

-32600 

953 

' -31400 

-31400 
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In-Sb  Assessment 

Group  Ill-antimonides  are  important  constituents  to  the  solid  solutions  used  for 
long  wavelength  photodetector  systems,  owing  to  the  direct  nature  of  their  bandgap 
structures  and  high  electron  mobilities.  In  particular,  InSb  has  the  lowest  band  gap 
energyand  highest  room  temperature  electron  mobility  of  any  other  III-V  compound 
semiconductor.  To  understand  the  electronic  and  structural  characteristics  of  solid 
solutions  comprising  InSb,  it  is  important  to  study  the  thermodynamic  properties  in  the  In- 
Sb binary  in  both  the  liquid  and  solid  phase  so  as  to  predict  multicomponent  solution 
behavior,  stability  of  solid  solution  and  the  nature  of  interactions  within  the  alloy.  The  In- 
Sb system  is,  however,  over-specified  experimentaUy  to  the  extent  that  it  can  be  confusing 
and,  at  the  same  time  erroneous  to  select  data  without  checking  for  consistency.  The 

t * 

( 

purpose  of  this  work  is  to  systematically  and  criticaUy  analyze  the  data  on  In-Sb  system  in 

determining  a consistent  dataset  that  can  explain  the  solution  thermodynamics  in  the  In-Sb 
system. 

Phase  diagram  data.  The  phase  diagram  for  the  In-Sb  system  has  been  reported  by 
several  investigators  [Pog49,  Liu52,  Hal63,  Gor83]  and  a simple  eutectic  type  of  phase 
diagram  was  reported  on  either  side  of  the  compound  [Pog49,  Liu52].  The  Sb-InSb 
eutectic  occurs  at  a composition  of  Xsb  = 0.68  and  a temperature  of  767  K.  The  In-InSb 
eutectic  is  practically  degenerate  at  Xsb  = 0.007  and  a temperature  of  427  K,  approximately 
1 K below  the  pure  In  melting  temperature.  InSb  is  the  only  III-V  compound  having  a 
melting  temperature,  Tf,  below  that  of  a component  element  (i.e.,  Sb).  Liu  and  Peretti 
[Liu52]  investigated  the  In-Sb  system  by  thermal,  x-ray  and  metallographic  methods  of 
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analysis.  The  solid  solubility  of  antimony  in  indium  was  less  than  xsb  ~ 0.0047  as  was  the 
solid  solubility  of  indium  in  antimony.  The  data  of  Pogodin  et  al.  [Pog49]  agrees  quite 
well  with  the  results  of  Liu  and  Peretti;  their  data  are  only  slightly  higher  in  temperature 
for  Sb  compositions  Xsb  = 0.40.  Hall  [Hal63]  performed  liquidus  measurements  on  the  In- 
rich side  of  In-Sb  binary.  His  liquidus  measurements  are  significantly  higher  than  those 
made  by  Pogodin  and  Liu  and  Peretti.  Sharma  et  al.  [Sha87],  in  their  assessment  work, 
have  pointed  out  that  the  dissolution  technique  used  by  Hall  was  prone  to  error,  especially 
at  low  solubilities  of  InSb  in  liquid  In.  However,  the  more  recent  work  by  Goryacheva 
[Gor83]  is  in  good  agreement  with  Liu  and  Pogodin’s  reported  values.  A summary  of  the 
phase  diagram  data  is  presented  in  Table  3-13. 

Enthalpy  data.  Isothermal  data  for  the  enthalpy  of  mixing  in  the  liquid  phase  have 
been  reported  by  Yazawa  et  al.  [Yaz68],  itagaki  and  Yazawa  [Ita75,  Ita77],  Predel  and 
Oehme  [Pre76],  Oehme  and  Predel  [Oeh78],  Wittig  and  Gehring  [Wit67],  and  Vecher  et 
al.  [Vec74].  Rosa  et  al.  [Ros80]  and  Oehme  and  Predel  [Oeh78]  measured  the  enthalpy  of 
mixing  at  several  temperatures  in  the  temperature  range  910  K to  1200  K.  The  literature 
values  are  in  fair  agreement  with  each  other,  however,  the  range  for  the  maximum 
negative  enthalpy  of  mixing,  AHm,  is  2345  to  3725  J/mol  in  this  temperature  range.  The 
enthalpy  of  mixing  is  not  symmetric  with  re.spect  to  composition;  the  maximum 
exothermic  AHm  being  located  in  the  range  of  Xsb  = 0.40  - 0.45. 

Predel  and  Oehme  [Pre76]  performed  a calorimetric  investigation  of  the  enthalpy 
of  mixing  of  liquid  In-Sb  alloys  over  the  entire  concentration  range.  In  their  work,  Rosa  et 
al.  [Ros80]  determined  the  mixing  enthalpies  at  a function  of  temperature  with  the  aid  of  a 
high  temperature  calorimeter.  The  authors  attributed  the  deviations  from  the  regular 
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solution  behavior  to  the  existence  of  In3Sb  and  InSb  associates.  Furthermore,  their  results 
indicate  that  the  maximum  AHmix  shifts  to  a lower  magnitude  at  increased  Sb  content  as 
the  temperature  is  increased.  Wittig  and  Gehring  [Wit67]  obtained  mixing  enthalpies  at 
973  K.  Vecher  et  al.  [Vec74]  measured  the  enthalpy  of  mixing  at  910  K by  quantitative 
differential  thermal  analysis  while  Yazawa  et  al.  [Yaz68]  reported  AH^x  values  using  the 
same  technique.  The  enthalpy  of  fusion  has  been  determined  by  a number  of  investigators 
[Ita77,  Mec71,  Nac58,  Sch58,  Lic69,  Hul73,  Ric65,  Gar77,  Bla70]  and  a summary  is 
presented  in  Table  3-15.  The  results  of  Sirota  and  Yushkevich  [Sir67],  Nikol’skaya 
[Nik59],  Abbasov  and  Mamedov  [Abb70],  and  Terpilowski  and  Trzebiatowski  [Ter60]  for 
the  standard  enthalpy  of  formation  are  shown  in  Table  3-16.  Also  tabulated  are  the 
experimental  techniques  adopted  by  the  several  investigators  [Sch58,  Gad61,  Nik59, 
Sch61,  Jen67a,  Jen67b,  Kle55,  Abb70,  Lun63,  Sir67,  Ter60]. 

Gibbs  energy  data.  There  are  five  reported  experimental  determinations  of 
component  activities  in  liquid  In-Sb  solutions.  Terpilowski  [Ter59]  and  Terpilowski  and 
Trzebiatowski  [Ter60]  measured  the  In  activity  in  an  electrochemical  cell  using  a double 
bromide  salt  and  found  negative  deviations  throughout  the  entire  concentration  range  in 
the  temperature  range  643  to  763  K.  Hoshino  et  al.  [Hos65],  using  a double  iodide  salt 
concentration  cell,  measured  the  In  activity  at  900  K in  the  liquid  phase.  Their 
measurements  compared  well  with  the  results  of  Terpilowski.  Both  of  the  measurements 
assumed  that  the  In  ion  that  was  transported  in  the  fused  salt  electrolytes  was  monovalent. 
This  prompted  Chatterji  and  Smith  [Cha73]  to  investigate  the  In  activity  in  a galvanic  cell 
using  a solid  electrolyte  in  which  the  oxidation  state  of  In  was  trivalent.  Four 
compositions  were  studied  in  the  In-rich  portion  of  the  phase  field. 
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Anderson  et  al.  [And84]  have  pointed  out  the  uncertainties  involved  in  Chatterji  et 
al.’s  experiment  and,  after  performing  a third-law  analysis,  indicated  their  data  to  be 
inconsistent  with  other  thermodynamic  data.  Anderson  et  al.  set  up  two  solid  state 
electrochemical  cells  and  measured  the  In  activity  at  973  K.  Their  work  shows  less 
negative  deviation  when  compared  to  other  investigators. 

Four  emf  studies  of  the  standard  Gibbs  energy  of  formation  of  InSb(s)  are  available 
in  literature  [Nik59,  Abb70,  Sir67  and  Ter60]  as  summarized  in  Table  3.16. 

Analysis  of  data  on  In-Sb  system.  As  can  be  seen  from  Table  3-18,  six  parameters 
were  required  to  adequately  describe  the  In-Sb  liquid  phase  by  the  Redlich-Kister 
polynomial.  Three  parameters  were  determined  to  be  sufficient  to  model  the  solid  InSb 
phase  as  a stoichiometric  line  compound. 

The  calculated  phase  diagram  is  shown  in  Figure  3-9  along  with  the  experimental 
data  from  the  literature  including  those  omitted  from  the  selected  data  set.  The  reported 
phase  diagram  data  for  In-Sb  are  fairly  consistent  and  were  thus  used  as  the  basis  to  check 
the  consistency  of  the  other  thermodynamic  quantities.  Moderate,  if  not  severe, 
discrepancy  exists  in  the  component  activity  and  enthalpy  of  mixing  data.  The  idea  was  to 
begin  with  less  parameters  for  the  liquid  phase  and  as  soon  as  the  optimization  indicated  a 
decreasing  (converging)  mean  square  of  error,  increase  the  number  of  floating  parameters. 
The  procedure  was  as  follows: 

(1).  The  specific  heat,  Cp  (12.97040  J/mol  K)  was  estimated  from  the  data  of  Lichter  and 
Sommelet  [Lic69],  and  Hultgren  et  al.  [Hul73]. 
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Table  3-13.  Summary  of  In-Sb  phase  diagram  data. 


Reference 

Method 

Tm  (K) 

Ranges 

T(K) 

Xsb 

Accuracy 
AT  (K) 
A Xsb 

[Pog49] 

476  - 879 
0.01  - 0.91 
0.001-0.003 

5 

[Liu52] 

Thermal,  X-ray 
Metallography 

798 

556  - 904 
0.006-1 
0.001-0.003 

3-5 

[Hal63] 

Source 

Dissolution 

Technique 

• • V : 

438  - 538 
0.001  - 0.003 

5 

0.001 

[Gor83] 

mm 

73 1 - 797 
0.25  - 0.45 
0.001-0.003 

5 
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Table  3-14.  Summary  of  the  enthalpy  of  fusion  of  In-Sb. 


Reference 

Method 

AHf  (InSb)  (kJ/g  atom) 

[Ita75] 

Quantitative  Thermal 
Analysis 

24.35  ± 0.50 

[Ita77] 

Quantitative  Differential 
Thermal  Analysis 

25.5  ± 1.0 

[Nac58] 

Drop  Calorimetry 

23.47  ± 0.80 

[Sch58] 

Calorimetry 

25.52 

[Lic69] 

Drop  Calorimetry 

23.88  ± 0.40 

[Hul73] 

Evaluated 

24.04  ± 0.80 

[Ric65] 

Differential  Thermal 
Analysis 

18.8  ±6 

[Gar77] 

Quantitative  Differential 
Thermal  Analysis 

24.27  ± 1.5 

[Bla70] 

.Drop  Calorimetry 

24.69  ± 0.4 
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Table  3-15.  Summary  of  the  standard  enthalpy  of  formation  of  InSb. 


Reference 

Method 

T(K) 

AHf  (InSb) 
(kJ/g  atom) 

[Gad61] 

Bomb  Calorimetry 

298 

16.28  ±0.17 

[Nik59] 

EMF 

298 

15.36  ±0.04 

[Sch61] 

Drop  Calorimetry 

298 

14.23  ± 1.3 

[Jen67a] 

Solution 

Calorimetry 

278 

14.90  ±0.21 

[Jen67a] 

Solution 

Calorimetry 

273 

15.27 

[Jen67b] 

Solution 

Calorimetry 

298 

15.31 

[Kle55] 

Tin  Solution 
Calorimetry 

723 

18.07  ±0.08 

[Sch58] 

Tin  Solution 
Calorim.etry 

273 

14.52  ±0.46 

[Abb70] 

EMF 

298 

17.15  ±0.84 

[Lun63] 

Tin  Solution 
Calorimetry 

273 

15.32  ±0.42 

[Lun63] 

Tin  Solution 
Calorimetry 

750 

17.71  ±0.42 

[Sir67] 

EMF 

703 

17.20 

[Ter60] 

EMF 

703 

18.83  ±0.88 
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Table  3-16.  Summary  of  Gibbs  energy  of  formation  of  InSb  for  the  reaction: 

In  (1)  -I-  Sb(s)  = InSb(s). 


Reference 

Method 

AGf  @ 643  K 

AGf  @ 763  K 

(kJ/mol) 

(kJ/mol) 

[And84] 

Derived 

18.80 

15.73 

[Ter59] 

Derived 

20.87 

18.19 

[Hos65] 

Derived 

18.43 

15.30 

[Abb70] 

EMF 

17.78 

14.73 

[Nik59] 

EMF 

18.34 

15.56 

[Sir67] 

EMF 

18.47 

15.50 

[Ter60] 

EMF 

19.20 

15.76 
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Table  3-17.  Summary  of  selected  dataset  for  In-Sb  system. 


Data  T ype 

References 

Phase  Diagram 

[Pog49,  Liu52,  Gor83] 

Enthalpy 

[Wit67,  Lic69,  Bla70,  Hul73,  Ita77, 
Oeh78,  Ros80] 

Gibbs  Energy 

[Hos65,  Cha73,  And84] 

Table  3-18.  Optimized  coefficients  for  Redlich-Kister  model  and  lattice  stability  of  InSb. 


n 

Ain.sb  (J/mol) 

Bin-sb  (J/mol  K) 

Cin-sb  (J/mol  K) 

0 

-25631.2 

-89.47429 

13.45816 

1 

-2115.4 

1.31907 

— 

2 

2908.9 

— 

— 

Table  3-19.  Invariant  equilibria  in  the  In-Sb  system. 


Reaction 

Eutectic 

Composition 

Eutectic 

Temperature 

Reference 

At  % Sb 

(K) 

L = In(s)  -1-  InSb(s) 

0.66 

427.8 

[Liu52] 

— 

428 

[Pog49] 

0.47 

427.4 

This  study 

L = InSb(s)  + Sb(s) 

68.3 

767  ±3 

[Liu52] 

70.4 

779  ±2 

[Pog49] 

68.8 

767.35 

This  study 
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Table  3-20.  Summary  of  Ga-Sb  phase  diagram  data. 


Reference 

Method 

Tn,(K) 

Ranges 

T(K) 

Xsb 

Accuracy 
AT  (K) 
Axsb 

[DeW86] 

Source 

Dissolution 

Technique 

648  - 923 
0.003  - 0.2 

3-5 

0.001-0.003 

[Mag68] 

Differential 

Thermal 

Analysis 

982.2  ± 0.3 

823  - 982 
0.08  - 1 

3-5 

0.001-0.003 

[Gre55] 

Thermal,  X-ray 
Metallography 

978.9  ± 1.0 

820  - 980 
0.06  - 1 

3-5 

0.001-0.003 

[Kos55] 

Thermal 

Analysis 

976 

850  - 980 
0.05  - 1 

5 

0.003 

[Hal63] 

Source 

Dissolution 

Technique 

520  - 960 
0 - 0.30 

5 

0.001-0.003 

[Cha77] 

Source 

Dissolution 

Technique 

670  - 820 
0.006  - 0.07 

3-5 

0.001 

[Gla58] 

— 

920  - 960 
0.4  - 0.6 

5-  10 

0.001-0.003 

[Ded78] 

— — — 

— 

725  - 870 
0.015-0.10 

3-5 

0.001-0.003 
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Table  3-21 . Summary  of  the  enthalpy  of  fusion  of  GaSb. 


Reference 

Method 

AHf  (GaSb)  (kJ/g  atom) 

[Sch58] 

d.t.a. 

25.1  ±2 

[Mec71] 

d.t.a. 

32.4  ± 1.3 

[Gar77] 

d.t.a. 

29.1  ± 1.7 

[Gam75] 

d.t.a. 

27.6  ± 1.7 

[Lun63] 

Drop  Calorimetry 

31.0 

[Cox67] 

Drop  Calorimetry 

33.5  ±3.2 

[Lic69] 

Drop  Calorimetry 

32.55  ±0.4 

[Bla70] 

Drop  Calorimetry 

30.54  ± 1.3 

[Gam75] 

Drop  Calorimetry 

30.10  ± 1.3 

[Hul73] 

(Evaluated) 

33.12  ±0.6 

Table  3-22.  Summary  of  the  standard  enthalpy  of  formation  of  GaSb. 


Reference 

Method 

T(K) 

AHf  (GaSb)  (kJ/g  atom) 

[Sch58] 

Calorimetry 

298 

20.79  ± 0.9 

[Sch59a] 

Calorimetry 

298 

20.79  ± 0.9 

[Abb64] 

EMF 

298 

19.66  ±2.1 

[Lun63] 

Calorimetry 

633  - 833 

23.43  ±2.1 
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Table  3-23.  Summary  of  Gibbs  energy  of  formation  of  GaSb  for  the  reaction: 

Ga(i)  + Sb(s)  = GaSb(s). 


Reference 

Method 

T(K) 

AHf  (GaSb) 

(kJ/g  atom) 

[Sch58] 

Calorimetry 

298 

19.08  ± 1.5 

[Abb64] 

EMF 

298 

18.83  ± 1.3 

[Sir67] 

EMF 

298 

20.54 

[Sha87] 

Assessment 

298 

19.35* 

(*  denotes  recalculated  value  using  SGTE  transformation) 


Table  3-24.  Summary  of  selected  dataset  for  Ga-Sb  system. 


Data  Type 

References 

Phase  Diagram 

[Gre55,  Mag68,  Cha77,  DeW86] 

Enthalpy 

[Lic69,  Hul73,  Gam75] 

Gibbs  Energy 

[Ger65,  Ber74,  And83,  And75,  Kat87] 

Table  3-25.  Optimized  coefficients  for  Redlich-Kister  model  and  lattice  stability  of 

GaSb(s). 


n 

Aca-sb  (J/mol) 

Bca-sb  (J/mol  K) 

Cca-sb  (J/mol  K) 

0 

-13953.8 

-61.45546 

9.62320 

1 

1722.9 

1.92588 

— 

2 

2128.3 

— 

— 
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Table  3-26.  Invariant  equilibria  in  Ga-Sb  system. 


Reaction 

Eutectic 

Composition 

At.  % Sb 

Eutectic  Temperature 
(K) 

Reference 

L = Ga(s)  + GaSb(s) 

— 

302.5 

[Kos55] 

— 

302.8 

[Gre55] 

L = GaSb(s)  + Sb(s) 

87 

856 

[Kos55] 

88.7 

862.8  ± 0.5 

[Gre55] 

— 

861.5  ±0.3 

[Mag68] 

88.5 

862.7 

This  study 
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(2) .  Using  experimental  enthalpy  and  Gibbs  energy  data,  optimized  coefficients  were 
used  to  predict  the  phase  diagram.  Inconsistent  thermodynamicdata  were  discarded. 

(3) .  The  cornbined  thermochemical  and  phase  diagram  data  set  was  taken  together  and 
the  coefficients  allowed  to  float.  The  coefficients  obtained  hereafter  were  tested  for 
variation  with  experimental  data.  The  best  set  of  coefficients  was  selected  which  provided 
the  optimal  fit  to  the  combined  data  set. 

The  optimized  phase  diagram  (Figure  3-9)  fits  the  experimental  data  of  Liu  and 
Peretti,  Pogodin  et  al.,  and  Goryacheva  quite  well.  The  data  by  Pogodin  et  al.  was 
consistently  higher  than  those  by  other  authors  and  the  optimized  values,  but  were  still 
retained  as  the  deviations  were  within  the  experimental  error  limits.  Hall’s  data,  on  the 
contrary,  showed  marked  deviation  with  respect  to  the  other  data  values  at  indium-rich 
compositions  in  the  phase  diagram.  Eliminating  Hall’s  data  resulted  in  a significant 
improvement  in  the  phase  diagram  fit  and  mean  square  error  value.  The  calculated  melting 
temperature  of  InSb  (Tm  = 800  K)  and  the  two  eutectic  temperatures  (427.1  K,  767.5  K) 
are  in  good  agreement  with  the  reported  values  as  shown  in  Table  3-19. 

The  enthalpy  of  mixing  experimental  data  compares  quite  weU  with  the  optimized 
plots.  As  can  be  seen  from  Figure  3-10,  the  AHm  curve  is  assymmetric  in  composition. 
The  degree  of  assymmetry  is  fairly  moderate  and  is  taken  into  consideration  in  the 
calculation  by  the  use  of  parameters  in  the  phase  description  of  the  In-Sb  liquid  phase. 
There  is  a strong  temperature  dependence  in  the  liquid  phase  enthalpy  as  indicated  by 
Oehme  and  Predel  [Oeh78].  This  temperature  dependence  is  reproduced  well  by  the 
coefficients  of  the  liquid  phase  as  is  clearly  evident  in  Figure  3-10.  The  data  by  Vecher  et 
al.  [Vec74]  are  significantly  higher  in  value  when  compared  to  Itagaki  [Ita75,  Ita77] 
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Figure  3-11.  Optimized  enthalpy  changes  for  the  reaction:  In(s,l)  + Sb(s,l)  = InSb(s,l). 
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and  Predel’s  data  [Pre76].  Their  data  was  thus  not  selected  for  the  optimization.  Wittig’s 
data  [Wit67]  at  973  K was  found  to  be  consistent.  The  data  reported  by  Itagaki  for  923  K 
temperature  were  determined  to  be  low,  but  were  still  retained  as  their  inclusion  did  not 
affect  the  mean  square  error  significantly.  Furthermore,  their  data  were  found  to  be 
consistent  with  the  phase  diagram  data.  The  fit  of  the  calculated  curves  to  the  excess  heat 
content  data  of  Lichter  and  Sommelet  is  excellent,  as  shown  in  Figure  3-11.  The  8(AH)  of 
the  solid  and  liquid  phase  data  in  the  plot  was  accomplished  by  subtracting  out  the  heat 
content  contributions  from  the  pure  elements.  The  component  Cp  data  required  in  the 
calculation  was  taken  from  Hultgren  [Hul73]  . The  relations  used  were: 

(a)  Below  In  Melting  Point 

8(AH)  = H,„-0.5H^-0.5H,,  ^3  5^ 

(b)  Above  In  Melting  Point  and  Below  InSb  Melting  Point 

^In  melting,In ) Hullgrcn  (3.6) 

8(AH) , as  above 

(c)  Above  InSb  Melting  Point 

Sm)  = H„(T)  - H„, (T„ ) - (T - T,„ ) ('■‘’‘"Cp  / 2-h‘“’^'’Cp  / 2)  (3.7) 

The  enthalpy  of  fusion  of  InSb  is  determined  to  be  25.32  ± 0.9  kJ/g  atom  and  it  compares 
well  with  the  values  listed  in  Table  3-14.  Similarly,  the  Gibbs  energy  of  formation  of 
InSb(s)  from  In(s)  and  Sb(s)  at  643  K and  763  K are  calculated  as  -20.56  and  -18.15 
kJ/mol  respectively.  For  the  reaction  involving  In(I)  and  Sb(s)  the  respective  values  are  - 
18.93  and  -15.60  kJ/mol  respectively  and  agree  well  with  reported  values  in  Table  3-15. 
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Figure  3-12  shows  the  calculated  partial  Gibbs  energy,  min,  as  a function  of 
composition  at  973  K as  compared  to  the  data  of  Anderson  et  al.  [And84]  and  Chatterji 
and  Smith  [Cha73].  The  calculated  curve  agrees  well  with  these  data  and  also  gives  an 
excellent  fit  to  the  data  of  Hoshino  et  al.  [Hos65],  though  their  reported  values  are  at  900 
K.  The  data  by  Terpilowski  [Ter59]  showed  significant  discrepancy. 

Ga-Sb  Assessment 

Similar  to  the  In-Sb  system,  the  Ga-Sb  sytem  has  been  well  characterized  in  terms 
of  phase  diagram  information.  The  enthalpy  of  mixing  and  Gibbs  energy  data  for  the  Ga- 
Sb  system  show  some  disagreement  in  the  available  literature  values,  partly  because  of  the 
differences  in  experimental  methods  used  and,  partly  because  of  the  experimental  problems 
encountered  in  measurements  of  these  properties.  Thus,  the  motivation  towards 
performing  an  assessment  of  the  Ga-Sb  system  is  to  determine  a complete  and  a consistent 
data  set  which  can  be  used  to  make  accurate  predictions  for  multicomponent  systems 
containing  the  binary  Ga-Sb  system. 

Phase  diagram  data.  The  phase  diagram  for  the  Ga-Sb  system  has  been  reported  by 
several  investigators  [Hal63,  DeW86,  Mag68,  Gre55,  Kos55,  Cha77,  Gla58,  Deg78]  and 
have  indicated  a simple  eutectic  type  of  phase  domain  for  the  system  [Mag68,  Gre55, 
Kos55]  with  negligible  solubility  of  either  element  in  the  other.  A brief  summary  of  the 
methods  employed  by  the  various  investigators  for  the  liquidus  measurements  is  given  in 
Table  3-20.  Also  shown  in  Table  3-20  are  the  temperature  and  composition  ranges 
investigated  by  these  authors  and  the  reported  (or  estimated)  accuracies  in  their  respective 
measurements.  The  Sb-GaSb  eutectic  occurs  at  a composition  of  Xsb  = 0.88  and  a 
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temperature  of  863  K,  while  the  Ga-GaSb  eutectic  is  nearly  degenerate.  The  melting 
temperature  values  as  measured  by  various  investigators  [Sha87,  Mag68,  Gre55,  Kos55, 
And83]  range  form  975  to  998  K;  and  the  generally  accepted  value  is  985  K [Hul73], 
Greenfield  et  al.  have  reported  the  data  of  Maglione  and  Potier  to  be  about  5 K below 
their  data,  while  the  data  of  Hall  is  as  much  as  5-10  K higher  than  their  data.  However,  as 
was  also  mentioned  for  the  case  of  In-Sb  system,  Sharma  et  al.  suspect  the  data  by  Hall 
stating  that  the  source  dissolution  technique  used  by  Hall,  was  prone  to  error  at  low  Sb 
compositions.  Likewise,  Hall  s data  is  too  high  when  compared  with  the  liquidus 
measurements  of  Chang  and  Pearson  [Cha77]  and  Dedegkaev  et  al.  [Ded78].  Glazov  and 
Petrov  [Gla58]  made  measurements  in  the  range  0.4  < Xsb  <0.6.  Their  values  are 
significantly  lower  than  the  other  reported  values.  Koster  and  Thoma  report  values  which 

are  consistently  lower  than  other  measurements  throughout  the  composition  range  of  their 
investigation. 

Enthalpy  data.  The  exothermic  enthalpy  of  mixing  in  the  liquid  phase  has  been 
reported  by  Yazawa  et  al.  [Yaz68],  Predel  and  Stein  [Pre71]  and  Gambino  and  Bros 
[Gam75]  at  a temperature  near  1000  K.  The  data  by  Predel  and  Stein  and  Gambino  and 
Bros  agree  quite  well,  however,  the  data  by  Yazawa  et  al.  show  significant  deviations 
from  the  previous  studies.  Predel  and  Stein  reported  the  maximum  exothermic  enthalpy  of 
mixing  at  Xsb  = 0.5.  Yazawa  et  al.  used  an  adiabatic  calorimeter  to  obtain  mixing 

enthalpies  and  fitted  the  data  to  a fourth  degree  polynomial  expansion  in  terms  of  Sb 
content. 

The  enthalpy  of  fusion  has  been  reported  by  several  investigators  [Mec71,  Sch58, 
Lic69,  Hul73,  Gar77,  Bla70,  Lun63,  Gam75  and  Cox67]  and  is  tabulated  in  Table  3-21. 
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Schottky  and  Bever  [Sch58]  and  Lundin  et  al.  [Lun63]  determined  the  enthalpy  of 
formation  of  GaSb.  These  values  are  listed  in  Table  3-22. 

Gibbs  energy  data.  The  Ga  activity  in  the  liquid  phase  has  been  measured  by 
electrochemical  techniques  at  around  1000  K [And83,  Pon75,  Dan70,  And75,  Ger65, 
Kat87].  The  activity  data  of  Pong  [Pon75],  Gerasimenko  et  al.  [Ger65],  and  Danilin  et  al. 
[Dan70]  show  relatively  large  negative  deviations  from  ideality.  On  the  other  hand,  Hsi- 
Hsiung  et  al.  [Hsi66]  and  Bergman  et  al.  [Ber74]  derived  the  activities  from  vapor 
pressure  measurements.  Bergman  et  al.  directly  determined  the  activity  of  Ga  by  means  of 
a Knudsen  cell  coupled  to  a mass  spectrometer.  Their  data  as  well  as  the  values 
determined  by  Anderson  [And75],  Katayama  et  al.  [Kat87]  and  Anderson  et  al.  [And83] 
show  moderate  negative  deviation  from  ideality  and  are  in  close  agreement.  In  their  work, 
however,  Katayama  et  al.  reported  moderate  deviation  of  the  data  by  Anderson  et  al.  from 
their  measured  values  in  the  high  and  low  composition  ranges,  i.e.,  xsb  = 0.2  and  xsb  = 0.8. 

The  discrepancy  in  the  degree  of  negative  deviation  of  the  Ga  activity  from  ideality 
for  a few  authors  [Hsi66,  Ber74]  compared  to  the  others  [And83,  Pon75,  Dan70,  And75, 
Ger65,  Kat87]  might  be  attributed  to  the  uncertainty  in  the  Ga  valency  for  the  chloride 
molten-salt  galvanic  cells  used.  The  possible  existence  of  GaCl  and  GaCb  in  addition  to 
GaCls  would  result  in  an  effective  charge  of  less  than  the  value  of  3 assumed. 

Anderson  et  al.  [And83]  found  the  Non-Random  Two  Liquid  (NRTL)  equation  to 
best  represent  the  activities  reported  in  their  work  while  also  performing  a consistency  test 
on  the  other  activity  data.  A summary  of  the  Gibbs  energy  of  formation  of  GaSb  is 
presented  in  Table  3-23. 
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Analysis  of  data  on  Ga-Sb  system.  The  selected  data  set  for  Ga-Sb  is  giyen  in 
Table  3-24.  Six  parameters  were  required  to  adequately  describe  the  liquid  phase  by  the 
same  Redlich-Kister  formalism  and  are  shown  in  Table  3-25,  while  four  parameters  were 
sufficient  for  the  description  of  the  solid  GaSb  phase  as  a stoichiometric  line  compound 

The  calculated  phase  diagram  is  shown  in  Figure  3-13  and  appears  to  be  in 
excellent  agreement  with  experimental  determinations  of  the  liquidus.  Howeyer,  there 
exists  marked  differences  in  the  melting  temperature  of  GaSb  amongst  the  yarious 
inyestigators  [Sha87,  Mag68,  Gre55,  Kos55,  And83].  Similar  to  the  In-Sb  system,  the 
Ga-Sb  system  exhibits  a fairly  consistent  phase  diagram  data  set  and  thus  was  used  as  a 
basis  for  a consistency  check  of  the  other  thermodynamic  data  of  the  system.  A similar 
process  scheme  was  adopted  for  the  Ga-Sb  system  as  for  the  In-Sb  system  mentioned  in 
the  previous  section  and  the  same  convergence  criteria  used.  An  initial  estimate  of  9.623 
J/mol  K was  employed  utilizing  the  specific  heat  data  of  Lichter  and  Sommelet  [Lic69]  for 
liquid  GaSb  and  that  of  the  liquid  elements  fi'om  Hultgren  et  al.  [Hul73]. 

As  is  evident  from  Figure  3-13,  the  optimized  phase  diagram  fits  the  data  of 
Greenfield  et  al.,  Maglione  and  Potier,  Cheng  and  Pearson  and  DeWinter  et  al.  quite  well. 
The  data  by  Koster  and  Thoma  and  Glazov  and  Petrov  are  at  least  5- 1 0 K lower  than  the 
rest  of  the  liquidus  data.  As  in  the  case  of  In-Sb  system,  the  data  by  Hall  showed 
significant  deviation  in  the  gallium-rich  region  of  the  phase  diagram,  and  for  a similar 
reason  was  eliminated  from  the  selected  dataset.  Table  3-26  indicates  good  agreement 
between  the  calculated  melting  temperature  of  GaSb  (Tm  = 984.8  K)  or  the  eutectic 
temperature  (862.7  K)  with  the  reported  values. 
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The  liquid  phase  enthalpy  of  mixing  data  [PreVl,  Gam75]  compare  well  with  the 
optimized  data  as  shown  in  Figure  3-14.  The  parameters  give  an  acceptable  assymmetry 
and  temperature  dependence  to  the  liquid  phase  enthalpy  to  compare  well  with  the 
experimental  measurements  of  Predel  and  Stein  [Pre71]  and  Gambino  and  Bros  [Gam75]. 

The  data  of  Yazawa  et  al.  [Yaz68]  show  discrepancy  with  the  other  two  authors. 
Their  data  provide  an  incorrect  temperature  dependency  of  the  mixing  enthalpy  and  were 
considered  to  be  inconsistent. 

A similar  set  of  calculations  was  performed  for  determining  the  excess  heat  content 
data  of  the  solid  and  liquid  phases  as  was  done  for  the  case  of  In-Sb  system;  an  excellent 
fit  is  obtained  as  shown  in  Figure  3-15. 

The  enthalpy  of  fusion  of  GaSb  is  calculated  to  be  34.5  ± 0.8  kJ/g  atom  which 
compares  well  with  literature  values  in  Table  3-21.  Furthermore,  the  Gibbs  energy  of 
formation  of  GaSb  from  Ga(s)  and  Sb(s)  at  298  K is  determined  to  be  -15.49  kJ/g  atom. 
SGTE  transformation  yields  a value  of  -15.58  kJ/g  atom  for  the  pure  components,  Ga(l) 
and  Sb(s),  which  agrees  fairly  well  with  the  data  in  Table  3-22. 

Figure  3-16  shows  good  agreement  of  the  partial  molar  Gibbs  energy  (mca)  data 
with  the  reported  values  [And83,  And75,  Ger65,  Kat87].  However,  the  data  of  Danilin 
und  Yatsenki  [Dan70]  do  not  match  well  with  the  optimized  curve  and  were  determined  to 
be  inconsistent  based  on  their  ability  to  generate  the  phase  diagram  and  fit  the  liquid  phase 

enthalpy.  As  mentioned  previously,  their  data  exhibited  significantly  negative  deviations 
from  ideality. 
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Figure  3-13.  Optimized  Ga-Sb  phase  diagram. 
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Figure  3-16.  Optimized  liquid  phase  Ga  activity  in  the  Ga-Sb  system. 
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Assessments  of  Ternary  Systems 


Ga-Sb-As  Assessment 

The  Ga-Sb-As  system  exhibits  a eutectic  valley  extending  from  the  Sb-rich  region 
side  of  the  binary  Ga-Sb  phase  diagram  to  the  As-rich  region  of  the  GaAs  phase  diagram. 
Given  the  high  vapor  pressure  of  As,  the  phase  diagram  determinations  have  been  limited 
to  As-dilute  compositions  [Ant70,  Woo58,  Ino70,  Fos72,  Gra73,  Nah77,  Sug77,  GraSO, 
Pes83,  Man86  and  Mik65]. 

The  pseudobinary  section  GaAs-GaSb  in  the  ternary  phase  diagram  is  a region  of 
extended  solid  solution  with  a large  region  of  immiscibility  in  the  solid  phase.  There  exists 
significant  disagreement,  however,  regarding  the  solid  phase  behavior  in  this  pseudobinary 
section.  A eutectic-type  phase  diagram  has  been  proposed  [Ino71]  in  the  GaAs-GaSb 
section,  though  more  recent  measurements  have  indicated  a miscibility  gap  in  the  solid 
solution  [Gra73,  Pes83,  Man86  and  Ish88]. 

Measurements  of  the  liquidus  have  generally  been  made  using  differential  thermal 
analysis  (D.T.A.),  while  those  on  the  solidus  have  been  made  using  microprobe  analysis  on 
solid  solutions  of  GaAs-GaSb  deposited  on  GaAs  substrates  by  liquid  phase  epitaxy  (LPE) 


or  by  X-ray  analysis. 
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Binary  systems.  The  binary  systems  Ga-As,  Ga-Sb  and  As-Sb  [Ans93,  Nga88] 
haye  been  assessed  preyiously,  and  the  optimized  coefficients  for  the  Gibbs  energy 
expressions  for  the  contained  phases  have  been  used  as  the  limiting  values  for  assessment 
of  the  ternary  system  Ga-Sb-As.  The  binary  systems  Ga-Sb  and  Ga-As  were  assessed 
[Ish88,  Sha89  and  Lee92]  as  also  by  this  work.  For  the  As-Sb  system,  the  thermodynamic 
parameters  were  taken  from  Ishida  et  al.  [Ish89a]. 

Solid  phases.  There  has  been  no  evidence  for  the  formation  of  a ternary  compound 
in  this  system.  For  the  pseudobinary  solid  solution,  Straumanis  and  Kim  [Str65]  showed 
that  Vegard's  law  was  obeyed  to  34  at.  % GaAs  Ifom  their  measurements  on  samples  of 
alloys  obtained  by  annealing  quenched  melts.  In  contrast,  Muller  and  Richards  [Mul64] 
observed  complete  solid  solution  over  the  entire  composition  range  in  flash  evaporated 
samples  and  lattice  parameter  measurements  indicated  that  Vegard's  law  was  obeyed  to  50 
at.  % GaAs  composition.  Later,  Gratton  and  Woolley  [Gra73]  determined  the  variation  of 
lattice  parameter  with  composition.  They  confirmed  the  validity  of  Vegard's  law  to  38  at. 
% GaAs  compositions.  At  compositions  greater  than  61  at.  % GaAs,  they  observed  slight 
negative  deviations  from  Vegard's  Law,  which  was  consistent  with  the  predicted  results  of 
Straumanis  and  Kim  [Stra65].  A summary  of  lattice  parameters  is  given  in  Table  3-27. 

Pseudobinarv  section.  The  liquidus  in  the  pseudobinary  section  has  been  fairly  well 
characterized  [Ino71,  Fos72,  Ish88,  Mik65]  and  the  data  are  in  agreement.  The  solid 
solutions  in  the  pseudobinary  section,  however,  have  been  the  subject  of  controversy. 
Inoue  et  al.  [Ino71]  reported  a eutectic  type  behavior  at  a composition  of  3.0  mole  % 
GaAs  composition  and  temperature  of  981  K,  while  other  investigators  [Gra73,  Pes83, 


135 


Man86,  Ish88]  found  that  a solid  solution  occured  over  the  entire  composition  range.  The 
critical  temperature  for  the  miscibility  region  lay  near  the  solidus  curve.  Given  the 
difficulty  in  achieving  equilibrium  in  the  solid  solution,  it  was  not  surprising  that  the  range 
of  immiscibility  varied  significantly  in  the  authors’  works.  Pessetto  et  al.  [Pes83], 
Stringfellow  [Str83],  and  Cooper  et.  al.  [Coo82]  determined  a symmetric  miscibility  gap, 
while  Gratton  and  Woolley  [Gra80],  Mani  et  al.  [Man86]  and  Ishida  et  al.  [Ish88] 
suggested  an  asymmetric  miscibility  gap. 

Invariant  equilibria.  The  invariant  reaction  from  liquid  to  solid  in  the  As-dilute 
pseudobinary  section  has  been  reported  to  be  peritectic  by  Stringfellow  [Str72],  Gratton 
and  Woolley  [Gra80]  and  Ishida  et  al.  [Ish83];  but  a eutectic  reaction  has  been  proposed 
by  Inoue  et  al.  [Ino71].  A peritectic  behavior  is  accepted  in  this  assessment. 

The  four  phase  reaction  L <=>  (Ga)  + Ga(As,Sb)'  + Ga(As,Sb)“  is  degenerate.  The 
calculated  mo  no  variant  equilibria  is  shown  in  Figure  3-17  and  the  reaction  scheme  is 
presented  in  Figure  3-18.  The  temperatures  and  compositions  of  the  various  binary  or 
pseudo-ternary  reactions  are  given  in  Table  3-28. 

Liquidus  surface.  Liquidus  data  have  been  reported  only  for  the  As-dilute 
pseudobinary  section  and  the  Ga-Sb  rich  region  in  the  ternary  system.  A number  of 
calculated  phase  diagrams  using  either  semi-empirical  models  or  optimized  coefficients 
have  appeared  in  the  literature.  Antypas  [Ant70]  calculated  the  ternary  phase  diagram 
using  Darken's  quadratic  formalism  for  the  ternary  liquid  and  a regular  solution  model  for 
the  solid  solution  and  compared  the  results  with  the  data  of  Woolley  [Woo62].  The 
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Figure  3-17.  Calculated  mono  variant  equilibria  in  the  Ga-As-Sb  system. 
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Figure  3-19.  Calculated  liquidus  lines  in  the  Ga-As-Sb  system. 
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analysis  of  Foster  and  Woods  [Fos72]  indicated  nearly  ideal  liquid  phase  behavior  with 
strong  positive  deviations  in  the  solid  solution.  Sugiyama  et  al.  [Sug77]  obtained  As 
solubility  data  from  weight  loss  measurements.  Gratton  and  Woolley  [GraSO]  found  large, 
three-phase,  liquid-solid-solid  fields  corresponding  to  the  two-phase,  solid-solid  fields  of 
the  pseudobinary  section.  They  also  determined  the  locus  of  the  cusp  in  the  liquidus  sheet. 
Calculated  liquidus  lines  are  shown  in  Figure  3-19. 

In-Sb-As  Assessment 

The  phase  diagram  of  the  In-Sb-As  system  is  dominated  by  the  pseudobinary  solid 
solution  formed  between  the  only  contained  binary  compounds  InAs  and  InSb.  These 
alloys  are  attractive  semiconducting  materials  due  to  their  narrow  and  direct  bandgap 
energy  (infrared  optoelectronic  devices)  and  high  electron  mobility  (high  speed  electronic 
devices).  Several  consistent  measurements  of  the  liquidus  surface  have  been  reported, 
however,  the  variation  of  the  solidus  has  been  controversial,  primarily  due  to  the  strain 
energy  associated  with  the  size  difference  between  As  and  Sb  atoms,  the  solid  solution 
exhibits  immiscibility.  The  measured  extent  of  this  immiscibility  has  varied  considerably, 
with  both  eutectic  and  isomorphous  type  phase  diagrams  reported.  Experimental 
difficulties  beyond  those  normally  encountered  with  temperature  and  composition 
measurement  include  low  diffusivities  in  the  solid.  As  vapor  loss,  a large  As  segregation 
coefficient,  and  a tendency  for  supercooling  in  the  liquid. 
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Table  3-27.  Summary  of  solid  phases  in  the  Ga-As-Sb  system. 


Phase 

Pearson 

Symbol/Prototype 

Lattice  Parameter 
(pm) 

Comments 

Ga(AsxSbi-x) 

cF8 

ZnS 

Sphalerite 

565.26  to  609.5 

Vegard’s  Law  valid 
for  GaSb-rich 
compositions  and 
slight  negative 
deviations  at  GaAs- 
rich  compositions 

GaAs 

same 

565.26 

X=1 

GaSb 

same 

609.5 

x=0 
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Table  3-28.  Summary  of  invariant  equilibria  in  the  Ga-As-Sb  system 


T(K) 

Type 

Phases 

Composition  At  % 

Reaction 

As  Ga  Sb 

L = (As)  -t-  GaAs 

Eutectic 

Liquid 

0.0  0.0 

As 

GaAs 

1001.6 

Liquid 

0.0 

L = (Sb)  -1-  GaSb 

Eutectic 

Sb 

GaSb 

L-t-GaAs^Sb,.,' 

Peritectic 

Liquid 

0.052  0.5  0.321 

= GaAs^Sb,./ 

GaAs^Sb,., 

302.8 

Liquid 

O 

• 

o 

1 

o 

• 

1 

o 

• 

o 

1 

L = (Ga)  + 

Degenerate 

(Ga) 

GaAs^Sb,_J  + 

GaAs,Sb,_, 

GaAs^Sb,./ 
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Binary  systems.  The  phase  diagrams  of  the  In-As  and  In-Sb  binary  systems  are 
centered  about  the  compound  InAs  or  InSb.  Simple  eutectic  type  phase  diagrams  are 
exhibited  on  both  sides  of  the  compounds.  Though  little  information  is  available  for  As- 
rich  solutions,  considerable  experimental  data  exists  for  the  compound  InAs  and  In-rich 
compositions  in  the  In-As  system.  The  In-Sb  system  has  been  extensiyely  studied, 
including  seyeral  consistent  measurements  of  the  liquidus  and  compound  properties  as 
well  as  the  liquid  phase  enthalpy  of  mixing  and  component  activities.  The  As-Sb  system 
shows  a continuous  solid  solution  and  is  less  well  studied  than  the  other  two  binary 
systems.  The  As-Sb  system  shows  a continuous  solid  solution  and  is  less  well  studied  than 
the  two  other  binary  systems.  The  binary  systems  In-Sb  and  In-As  were  assessed  by 
Ishida  et  al.  [Ish88].  For  the  As-Sb  system,  the  thermodynamic  parameters  are  taken  from 
[Ans94]. 

Solid  phases.  No  ternary  compounds  have  been  reported  in  the  In-Sb-As  ternary 
system.  The  structural  characteristics  of  the  pseudobinary  InAsxSbi.*  have  been  studied  by 
a number  of  investigators.  The  solution  crystallizes  in  the  zinc  blende  structure  and  shows 
a slight  deviation  from  Vegard’s  Law.  Woolley  and  Warner  [Woo64]  improved  the 
accuracy  of  their  original  measurements  [Woo58]  of  the  lattice  parameter  by  annealing 
compressed  powders  for  6 months  at  a temperature  10  to  15°C  below  the  relevant  solidus 
temperature.  Slight  negative  deviations  from  Vegard’s  Law  were  observed  for  InAs 
compositions  greater  than  70  at.  % InAs,  while  somewhat  larger  positive  deviations  were 
reported  for  InSb-rich  compositions.  The  cross  over  composition  (~  30  at.  % InAs) 
corresponds  to  a minimum  in  the  bandgap  energy,  suggesting  that  considerable  ordering 
exists  on  the  group  III  sublattice  [Sem74].  The  room  temperature  lattice  parameter 
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measurement  of  Muller  and  Richards  [Mul64]  on  a flash  evaporated  equimolar  film  is  in 
good  agreement  with  the  results  of  Woolley  and  Warner  [Woo64].  A summary  of  the 
lattice  parameters  is  given  in  Table  3-29. 

Pseudobinarv  system.  The  pseudobinary  InAs-InSb  section  of  this  ternary  system 
has  been  examined  by  a number  of  investigators.  Thermal  analysis,  X-ray,  and 
metallographic  studies  [Shi53,  Shi56]  indicated  a eutectic-type  phase  diagram  that  was 
nearly  degenerate  with  respect  to  InSb.  The  extrapolated  InAs-rich  solidus  gave  an 
estimated  2 at.  % InSb  maximum  solubility  at  the  eutectic  temperature.  In  contrast,  the 
solidus  determination  by  X-ray  measurements  [Woo58]  indicated  an  isomorphous-type 
phase  diagram  with  a broad  liquidus-solidus  gap.  This  finding  was  later  confirmed  by  near 
equilibrium  liquid  phase  epitaxial  growth  studies  at  constant  temperature  [StrVlJ.  With 
the  exception  of  the  solidus  data  of  Shih  and  Peretti  [Shi53,  Shi56],  the  measured  values 
of  the  more  recent  solidus  data  are  in  good  agreement. 

Liquidus  measurements  of  the  pseudobinary  phase  field  have  been  reported  by 
several  investigators  [Woo58,  Str71,  Shi53  and  Gro83].  Although  Gromova  et  al. 
[Gro83]  extrapolated  their  DTA  measurements  to  zero  cooling  rate  to  account  for 
supercooling  of  the  melt,  their  four  liquidus  temperatures  for  Sb-rich  melts  were 
somewhat  below  the  values  reported  by  other  investigators.  Further  evidence  for  the 
formation  of  a continuous  solid  solution  is  found  in  the  ability  to  grow  the  full 
compositional  range  of  solutions  by  conventional  film  growth  techniques  such  as  organo- 
metallic  chemical  vapor  deposition  (MOCVD)  [Nat81]. 
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Table  3-29.  Summary  of  solid  phases  in  the  In-As-Sb  system. 


Phase 

Pearson 

Symbol/Prototype 

Lattice  Parameter 
(pm) 

Comments 

In(AsxSbi-x) 

cF8 

ZnS 

Sphalerite 

605.84  to  647.94 

Slight  posittive 
(InSb-rich)  and 
slight  negative 
(InAs-rich) 
deviations  from 
Vegard's  Law 

InAs 

same 

605.84 

X=1 

InSb 

same 

647.94 

x=0 

Table  3-30.  Summary  of  invariant  equilibria  in  the  In-As-Sb  system. 


T(K) 

Reaction 

Type 

Phases 

Composition  At  % 
As  In  Sb 

998.6 

L = (As)  + InAs 

Eutectic 

Liquid 

As 

InAs 

0.85  0.15  0.0 

763.0 

L = (Sb)  + InSb 

Eutectic 

Liquid 

Sb 

InSb 

0.0  0.288  0.702 

778.0 

Critical 

0.334  0.5  0.166 

InAs,Sb,.J  = 
InAs^Sb,,," 

Point 

InAs.Sb,., 
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Part  of  the  reason  for  the  discrepancy  in  the  solidus  is  related  to  the 
difficulty  in  achieving  equilibrium  at  the  relatively  low  temperatures  of  solid  solution 
existence.  Although  no  solid  solution  thermodynamic  data  are  available,  semi-empirical 
and  assessment  procedures  have  been  used  to  estimate  the  interaction  energy  in  the  solid 
solution.  Using  these  procedures,  solid  solution  immiscibility  was  predicted  above  room 
temperature  with  a calculated  critical  temperature  ranging  from  well  below  to  just  above 
the  solidus  curve  [Str69,  Str72,  Kau81,  Str82,  Ona82a,  Ona82b,  Pat87,  Ich89].  The 
results  of  EXFAS  measurements  show  that  the  In-V  bond  lengths  are  close  to  their  binary 
compound  values  as  the  composition  varies.  On  this  basis,  the  stored  elastic  strain  energy 
has  been  suggested  to  be  the  major  contribution  to  the  thermodynamic  enthalpy  of  mixing. 
The  existence  of  solid  solution  immiscibility  has  recently  been  experimentally  confirmed  by 
the  electron  probe  microanalysis  of  [Ish89].  Their  data  indicate  an  asymmetric  miscibility 
gap  with  the  critical  temperature  located  just  below  the  solidus  curve  at  an  InAs-rich 
composition  [Ish89]. 

Invariant  equilibria.  The  four  phase  reaction  L<=>  (In)-Hln(As,Sb)‘  +In(As,Sb)“ 

is  degenerate.  The  calculated  monovariant  equilibria  is  shown  in  Figure  3-20  and  the 
reaction  scheme  is  presented  in  Figure  3-2 1 . The  temperatures  and  compositions  of  the 
various  binary  or  pseudo-ternary  reactions  are  given  in  Table  3-30. 

Liouidus  surface.  In  addition  to  the  pseudobinary  liquidus  measurements  cited  in 
the  previous  section,  several  other  determinations  of  the  ternary  liquidus  have  been 
reported  in  the  literature.  Extensive  measurements  of  the  constitution  of  the  In-Sb-As 
system  were  performed  by  Shih  and  Peretti  [Shi54,  Shi56]  using  thermal  analysis  and 
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Figure  3-20.  Calculated  monovariant  equilibria  in  the  In-Sb-As  system. 
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Figure  3-22.  Calculated  liquidus  lines  in  the  In-Sb-As  system 
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metallographic  techniques.  In  these  studies  [Shi54],  the  InAs-Sb  section  as  well  as  three 
sections  in  the  In-InAs-InSb,  four  sections  in  the  Sb-InAs-InSb,  and  two  sections  in  InAs- 
Sb-As  subsystems  [Shi56]  were  examined.  Dissolution  experiments  were  performed  by 
Stringfellow  and  Greene  [Str71]  for  solutions  dilute  in  As  at  temperatures  typically 
encountered  in  liquid  phase  epitaxy  (LPE).  Liquidus  data  were  also  obtained  at  773  and 
793  K by  direct  visual  observation  of  the  dissolution  of  a ternary  crust  [AbrSl].  In 
addition  to  the  pseudobinary  liquidus  measurements  [Gro83],  these  authors  also 
investigated  six  compositions  dilute  in  As.  It  is  well  established  in  the  literature  that  a 
significant  amount  of  supercooling  (as  great  as  50°C)  can  exist,  even  when  using  low 
cooling  rates  or  a solid  substrate  as  in  the  LPE  experiments.  In  view  of  this  difficulty,  the 
As-dilute  liquidus  data  are  in  surprisingly  good  agreement.  Calculated  liquidus  lines  are 
shown  in  Figure  3-22. 

Solidus  surface.  Given  the  relatively  low  melting  temperature  of  the  elements  (In: 
429.75  K,  As:  1090  K,  Sb:  903.78  K)  and  the  compounds  (InAs:  1211  K,  InSb:  800  K) 
and  the  high  vapor  pressure  of  As,  solidus  measurements  have  been  difficult.  This  is 
apparent  from  the  discussion  previously  given  for  the  extent  of  the  miscibility  gap  in  the 
pseudobinary  section.  A considerable  amount  of  solidus  data  was  collected  [Shi54, 
Shi56],  though  those  data  involving  the  InAsxSbi.*  solid  solution  might  be  questioned  in 
view  of  the  subsequent  pseudobinary  solidus  measurements.  Strinfellow  and  Greene 
[Str71]  reported  solidus  data  in  the  temperature  range  793  to  1002  K by  measuring  the 
composition  of  films  grown  by  an  isothermal  LPE  process.  Similarly,  Abrikosov  et  al. 
[Abr81]  measured  the  composition  of  LPE  grown  films;  the  data  being  in  excellent 


agreement  with  [Str71]. 
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Solution  thermodynamics.  There  have  been  no  enthalpy  measurements  in  either  the 
liquid  or  solid  solution  phases.  Sorokina  et  al.  [Sor84]  has  reported  Knudsen  effusion 
vapor  pressure  measurements  for  a wide  range  of  compositions.  Rather,  the  solution 
thermodynamics  of  this  ternary  system  has  largely  relied  on  the  results  of  assessment  work 
or  semi-empirical  models  [Str69,  Ste70,  Str72,  BatSl,  Kau81,  Ona82a,  Ona82b,  Str82, 
Pat87,  Ish89a  and  Ish89b]. 


CHAPTER  4 
CONCLUSIONS 


Compound  semiconductor  systems  have  been  studied  experimentally.  The  binary 
arsenides  and  antimonides  were  the  subjects  of  specific  interest  to  this  study.  The  reason 
for  this  was  that  these  systems  determine  the  properties  of  higher  order  semiconductor 
systems  to  which  they  form  important  constituents.  Accurate  predictions  of 
multicomponent  properties  are  possible  only  through  good  understanding  of  properties  in 
these  binary  limits. 

The  present  work  was  successful  in  conducting  direct  measurement  of 
thermodynamic  properties  through  the  use  of  solid  electrolytic  cells.  Containment  of 
arsenic  vapors  was  the  most  difficult  part  in  these  experiments.  Another  problem  was  due 
to  the  long  diffusion  times  required  for  attaining  equilibrium  in  the  solid  phases  of  these 
systems.  However,  carefully  designed  electrochemical  cells  and  the  use  of  B2O3  powder  in 
the  electrodes  were  able  to  achieve  this.  This  was  evident  through  the  stability  of  the 
electrochemical  cells  during  operation  and  the  return  of  the  cells  to  equilibrium  condition 
upon  applying  electrical  and  thermal  perturbation  to  the  experimental  cells. 

The  results  indicate  negative  deviations  for  Group  III  activity  in  the  binary  liquid 
melts,  which  are  consistent  with  theoretical  predictions.  For  the  binary  systems  In- As,  Ga- 
As  and  Sb-Te,  the  experimental  results  on  Gibbs  energy  data  are  first  of  their  kind. 
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Activity  of  In  along  the  liquidus  in  the  binary  In- As  phase  diagram  was  obtained  for  the 
first  time  up  to  a temperature  of  1 1 50  K. 

For  Ga-As  system,  the  Gibbs  energy  of  formation  of  GaAs  does  not  exist  above 
743  K.  In  this  study,  the  experiments  were  successfully  extended  to  temperature  of  up  to 
973  K.  The  values  of  the  Gibbs  energy  of  formation  of  GaAs  at  these  high  temperatures 
were  checked  for  consistency  using  the  second  and  third  laws.  The  calculated  value  of  the 

standard  enthalpy  of  formation  of  GaAs  at  298. 15K,  A°Hf(298.15  K),  was  in  good 

agreement  with  experimental  data  [Mar74  and  Yam94].  Taking  into  due  consideration  the 
uncertainities  involved  in  the  specific  heat  data,  Cp  and  its  thermal  dependence,  it  was 
concluded  that  the  measured  values  were  consistent  and  within  acceptable  limits  of 
experimental  errors. 

A similar  anaysis  was  done  on  the  Gibbs  energy  data  for  InAs  in  the  In- As  system. 
The  analysis  showed  excellent  agreement  with  experimental  data  for  [Sha63,  Sch58  and 
Yam94].  The  existing  discrepancy  between  the  experimental  and  calculated  values  of  the 
standard  enthalpy  of  formation  of  InAs  at  298. 15K,  A°Hf(298.15  K),  indicates  random 

errors  inherent  to  experimental  measurement,  and  systematic  errors  inherent  to  the  choice 
of  specific  heat,  Cp,  data  used  in  the  calculations  by  numerous  aouthors.  The  error 
analysis  as  presented  in  chapter  2 shows  that  the  dominant  effect  is  that  of  the  temperature 
dependence  of  the  specific  heat,  Cp,  data. 

The  results  obtained  from  experiment  were  analysed  for  consistency  with  other 
thermochemical  and  phase  diagram  data.  The  systems,  Ga-As,  In-As,  In-Sb  and  Ga-Sb 
were  thermodynamically  assessed  in  the  present  study.  Such  data  was  utilized  to  perform 
phase  equilibria  calculations  and  predict  data  in  the  relatively  difficult  region  of 
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experimentation,  namely  the  Group  V-rich  region  of  the  phase  diagram.  Optimizations  of 
this  nature  were  able  to  define  the  necessary  binary  limits  for  ternary  and  higher  order 
semiconductor  systems,  for  which  data  is  usually  scarce  or  non-existent.  Such  predictions 
were  successfully  demonstrated  in  two  ternary  systems,  namely  Ga-As-Sb  and  In-As-Sb. 

Another  aspect  of  this  study  was  to  understand  the  pseudobinary  section  GaAs- 
InAs.  Thermodynamic  study  of  the  GaAs-InAs  pseudobinary  phase  diagram  and  the  Ga- 
In-As  system  provides  reliable  data  for  use  in  different  processes  such  as  liquid  phase 
epitaxy,  molecular  beam  epitaxy,  and  chemical  vapor  deposition,  employed  for  production 
of  multilayer  devices  with  GaxIni.xAs  solid  solutions.  In  fact,  near  equilibrium  epitaxial 
growth  conditions  are  proven  to  yield  devices  of  optimum  electrical  properties.  A 
discussed  in  chapter  2,  the  existence  of  a miscibility  gap  in  the  GaAs-InAs  pseudobinary 
solid  phase  region  has  been  debated  in  the  past.  The  experiments  in  the  present  work 
support  the  existence  of  a miscibility  gap,  and  estimate  the  critical  temperature  to  be 
between  950  K and  1050  K.  More  experiments  in  the  solid  phase,  however,  are  required 
to  characterize  the  solid  solution  of  the  alloys  in  terms  of  the  miscibility  range  through 
activity  measurements  at  different  alloy  compositions. 
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